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Chapter  1 


Introduction 


With  high  resolution  traffic  measurements  from  a  variety  of  networking  applica¬ 
tions  becoming  widely  available,  there  has  been  renewed  interest  in  understand¬ 
ing  the  statistical  nature  of  emerging  network  traffic.  Large  data  sets,  obtained 
from  Ethernet  LANs,  VBR  video  sequences,  ftp,  telnet  and  WWW  applications 
in  WANs,  have  been  extensively  studied  and  results  point  to  the  conclusion  that 
real  traffic  is  very  bursty,  exhibiting  great  variability  over  extended  periods  of 
time  that  are  much  longer  than  previously  expected.  As  these  dependencies  may 
have  a  pronounced  effect  on  performance  [40],  they  should  be  taken  into  account 
when  modeling  network  traffic  for  buffer  and  link  provisioning,  or  for  evaluating 
scheduling  policies. 

In  the  recent  literature  on  traffic  modeling  such  persistent  correlations  are  often 
reported  to  be  best  captured  by  stochastic  processes  that  are  long-range  depen¬ 
dent  [4,  12,  22,  23,  36,  48].  Roughly  speaking,  this  happens  when  the  data  stream 
displays  correlations  which  span  multiple  time  scales,  and  which,  despite  being 
individually  small,  decay  in  such  a  slow  hyperbolic-like  manner  as  to  be  consid¬ 
ered  non-summable.  In  other  cases  (e.g.,  the  studies  of  VBR  video  traces  [30,  34]) 
sample  autocorrelation  functions  are  found  to  conform  with  a  more  general  subex- 
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ponential  decay.  In  any  case,  the  observed  dependencies  cannot  be  exclusively 
attributed  to  classical  Markovian  models,  with  bounded  exponential  moments.  A 
consensus  seems  to  be  emerging  to  the  fact  that  non-traditional  stochastic  models 
should  be  considered;  these  will  most  likely  play  an  increasing  role  in  capturing 
the  dynamics  of  traffic  that  networks  are  expected  to  carry  in  the  near  future. 

Going  beyond  the  statistical  findings  mentioned  above,  an  interesting  line  of 
current  research  focuses  on  quantifying  the  consequences  of  high  variability  and 
dependence  on  network  performance.  The  initial  experimental  work  in  [18]  indi¬ 
cates  that  the  impact  of  long-range  dependencies  is  adverse  and  significant,  yet 
evidence  to  the  contrary  also  exists  [24,  38,  52],  suggesting  that  in  many  practical 
situations  queueing  measures  are  not  seriously  affected.  In  addition,  the  traffic 
measurement  studies  have  generated  interest  in  queueing  systems  with  correlated 
arrival  processes.  Few  analytical  results  are  currently  available  for  queues  with 
long-range  or  subexponentially  dependent  arrivals.  These  include  the  fractional 
Brownian  motion  model  of  Norros  [43],  fractional  Gaussian  noise  [1],  the  popular 
independent  on/off  source  model  with  Pareto  activity  periods  [9,  28]  and,  more 
recently,  the  multiplexed  on/off  sources  [29].  In  all  these  cases  buffer  overffow 
probabilities  display  a  slow  non-exponential  decay;  this  is  in  sharp  contrast  with 
the  exponential  tails  that  typically  characterize  queues  with  short-range  dependent 
Markovian  inputs.  Moreover,  a  closer  examination  confirms  that  classification  in 
terms  of  the  short  vs.  long-range  dependent  nature  of  the  traffic  is  often  insuffi¬ 
cient:  Both  within  the  short  and  the  long-range  dependent  regime  further  details 
matter,  and  vastly  different  queueing  behaviors  arise. 
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1.1  From  on/off  to  M\G\oc  sources 

In  this  dissertation,  we  consider  the  class  of  discrete-time  M\G\oo  input  processes. 
An  M\G\oo  input  process  is  understood  as  the  busy  server  process  of  a  discrete- 
time  infinite  server  system  fed  by  a  discrete-time  Poisson  process  of  rate  A  (cus¬ 
tomers/slot)  and  with  generic  service  time  a.  Such  M\G\oo  processes  can  account 
both  for  short  and  long-range  dependent  behaviors,  with  the  correlation  patterns 
controlled  through  a  [Proposition  2.1.1],  Furthermore,  asymptotic  self-similarity 
arises  when  a  is  Pareto-like,  i.e. ,  has  a  regularly  varying  tail  of  the  form  (3.9). 
M\G\oo  processes  have  already  been  used  by  Paxson  and  Floyd  to  successfully 
model  WAN  traffic  [48].  However,  perhaps  the  most  convincing  justification  sup¬ 
porting  their  use  as  plausible  traffic  models  is  provided  by  the  following  limiting 
result  [37]:  Consider  M  identical  and  independent  sources,  with  alternating  inde¬ 
pendent  emission  and  silence  periods.  Assume  that  during  its  “on”  periods  each 
source  generates  information  at  a  constant  rate  of  one  unit  per  time  slot,  while 
during  the  “off”  periods  it  remains  inactive.  We  alternatively  view  these  “on” 
periods  as  corresponding  to  information  sessions.  Denote  by  a  the  duration  (in 
slots)  of  the  generic  information  session  and  allow  the  duration  of  the  generic  “off” 
period,  denoted  by  T0q,  to  depend  on  the  number  of  sources,  i.e.,  T0q  =  T0s(M). 
The  resulting  total  session  arrival  rate  A (M)  is  given  by 

M 

MM)  =  E  [a]  +  E  [Toff(M)] ' 

Let  the  number  of  sources  M  go  to  infinity,  while  simultaneously  reducing  the 
number  of  sessions  of  an  individual  source,  so  that  the  aggregate  session  arrival  rate 
remains  finite.  This  can  be  achieved  by  selecting  T0g(M )  such  that  E  \T0q(M)\  = 
M/A  for  some  A  >  0,  in  which  case  lim  A (M)  =  A.  Consider  now  the  process 

M—>  oo 
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which,  in  each  time  slot,  records  the  total  number  of  newly  initiated  “on”  periods. 
This  converges,  as  M  goes  to  infinity,  to  a  discrete-time  Poisson  process  with 
session  arrival  rate  A  (sessions/slot).  Since  a  is  the  generic  session  duration  random 
variable  (rv),  we  readily  identify  the  M\G\oo  busy  server  process  with  parameters 
(A,  a)  as  the  limiting  process  counting  the  number  of  active  sessions  at  any  given 
time  slot.  Hence,  the  class  of  M | G | oc  processes  is  one  that  naturally  arises  from  a 
Poisson  superposition  scheme  of  infinitely  many  simpler  on/off  sources. 

1.2  Summary  and  discussion 

As  shown  in  [14,  39,  46,  47],  M\G\oo  processes  induce  a  wide  variety  of  asymptotic 
behaviors  for  the  buffer  probabilities  at  a  multiplexer  with  constant  release  rate. 
In  particular,  when  a  has  a  regularly  varying  tail  -  the  M\G\oo  process  is  now 
asymptotically  self-similar  -  the  buffer  asymptotics  are  hyperbolic  in  nature,  in 
stark  contrast  with  the  Weibullian  tails  induced  by  fractional  Gaussian  noise  (or 
fractional  Brownian  motion)  [43] .  A  key  contribution  of  this  dissertation  is  to  elu¬ 
cidate  the  noted  difference  in  buffer  asymptotics  between  M\G\oo  and  fractional 
Gaussian  noise  inputs  by  further  exploring  this  discrepancy  in  the  heavy  traffic 
regime.  One  might  expect  that,  with  asymptotically  identical  correlation  patterns, 
both  models  necessarily  have  a  heavy  traffic  characterization  in  terms  of  fractional 
Brownian  motion,  in  very  much  the  same  manner  that  different  short-range  depen¬ 
dent  models  eventually  collapse  to  a  single  description  involving  Brownian  motion. 
However,  this  turns  out  not  to  be  the  case. 

In  Chapter  3  we  show  that,  under  short-range  dependence,  the  class  of  M\G\oo 
inputs  belongs  to  the  domain  of  attraction  of  the  standard  Brownian  motion,  as 
expected.  However,  under  long-range  dependence,  with  a  belonging  to  the  domain 
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of  attraction  of  a  non-normal  stable  law,  the  M\G\oo  process  is  not  attracted  to  a 
fractional  Brownian  motion,  but  instead  to  a  non-Gaussian,  a  stable'  Levy  motion 
which  is  1/a  self-similar.  As  a  consequence,  the  distribution  of  the  heavy  traffic 
queue  length  is  given  by  a  Mittag-Leffler  function,  thus  displaying  not  a  Weibullian, 
but  a  Pareto  tail,  with  power  1  —  a  [Theorem  3.4.3].  These  results  undescore  the 
fundamentally  different  nature  of  the  long-range  dependent  M\G\oo  process  (when 
compared  to  fractional  Gaussian  noise),  and  also  point  to  the  fact  that  fractional 
Brownian  motion  does  not  necessarily  play  for  long-range  dependence  the  same 
key  role  that  standard  Brownian  motion  assumes  under  short-range  dependence. 
Within  long-range  dependence,  there  seems  to  be  a  choice  for  distinct  modeling 
possibilities,  and  it  is  not  at  all  difficult  to  find  rather  simple,  potentially  useful 
traffic  models  that  are  attracted  to  non-Gaussian  limits. 

In  Chapter  4  we  shift  attention  to  the  light  traffic  regime.  This  refers  to  the 
limiting  situation  where  the  traffic  intensity  approaches  zero.  Noting  that  the 
M\G\oo  process  is  Poisson  driven,  we  apply  the  Reiman-Simon  theory  [49,  50,  51] 
to  obtain  information  in  the  form  of  derivatives  of  system  quantities  with  respect 
to  the  intensity  of  the  driving  Poisson  process,  when  this  intensity  tends  to  zero 
[Propositions  4.2.3,  4.2.4],  In  addition,  when  the  “on”-state  rate  of  each  con¬ 
stituent  on/off  source  exceeds  the  multiplexer  release  rate  explicit  expressions  for 
the  expected  queue  size  become  available.  These  results  quantify  the  differences 
between  the  gradual  M\G\oo  inputs  and  the  point  arrivals  of  a  classical  GI\GI\l 
queue,  and  suggest  a  classification  of  the  light  traffic  behavior  of  the  buffer  con¬ 
tent  distribution  in  terms  of  the  short-  vs  long-  range  dependent  property  of  the 
M\G\oo  process  [Corollaries  4.3.1,  4.3.2],  However,  in  light  traffic  further  subcases 
arise  depending  on  comparisons  of  the  on/off  source  “on” -state  rate  and  the  server 
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capacity.  Thus,  in  general,  the  correlation  function  and  the  short-  vs  long-range 
dependence  property  of  the  M\G\oo  inputs  are  not  the  sole  factors  that  impact 
performance. 

Work  on  queueing  analysis  under  long-range  dependence  appears  to  have  been 
initiated  by  Norros  [43],  where  the  presence  of  fractional  Brownian  motion  is  pos¬ 
tulated.  This  line  of  inquiry  is  further  pursued  in  [58],  while  in  [9]  Brichet  et  al. 
show  how  fractional  Brownian  motion  can  arise  from  a  Gaussian  superposition 
scheme  of  infinitely  many  on/off  sources  with  heavy  tailed  on/off  periods.  In  the 
limiting  setup  of  [9]  the  sources  are  “small”,  i.e.,  the  peak  rate  of  the  individual 
source  becomes  infinitely  small  in  comparison  with  the  multiplexer  capacity.  In 
view  of  the  fact  that  M\G\oo  processes  arise  from  a  different  superposition  scheme 
of  infinitely  many  on/off  sources  [37],  where  the  peak  on/off  source  rate  is  “large”, 
i.e.,  remains  comparable  to  the  link  capacity,  it  is  not  too  surprising  that  these 
lead  to  a  different  heavy  traffic  limit  involving  Levy  motions.  More  recently,  in  [8], 
a  heavy  traffic  limit  of  this  type,  giving  rise  to  a  Mittag-Leffler  function,  is  ob¬ 
tained  in  the  standard  GI\GI\1  queueing  setup,  for  the  case  where  the  service  time 
distribution  is  heavy  tailed.  Heavy  traffic  results  similar  and  related  to  the  ones 
given  here  have  also  been  reported  in  [33],  where  only  convergence  of  finite  dimen¬ 
sional  distributions  is  announced.  The  conclusions  discussed  here  were  obtained 
independently,  and  were  summarized  in  the  conference  paper  [61]. 

The  asymptotic  characterization  of  the  queue  size  distribution  in  the  heavy 
and  light  traffic  regimes  is  exploited  in  Chapter  5.  By  suitably  interpolating  be¬ 
tween  the  two  extremes  we  derive  approximations  to  the  queue  size  distribution, 
applicable  to  all  traffic  intensities.  For  some  common  choices  for  the  session  du¬ 
ration  distribution  G  the  approximants  assume  a  simple  final  form.  The  accuracy 
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of  the  proposed  expressions  as  well  as  pitfalls  of  this  technique  are  discussed  via 
several  numerical  examples  [Section  5.4],  In  Chapters  2  and  4  we  make  occasional 
detours  and  study  simpler  queueing  systems;  these  help  us  obtain  exact  results 
[Propositions  4.3.2,  4.3.3]  and  establish  stochastic  comparisons  that  can  provide 
bounds  whenever  exact  expressions  are  not  available.  In  Appendices  A  and  B  we 
have  summarized  several  needed  facts  concerning  functions  of  regular  variation  and 
stochastic  orderings. 

A  few  words  about  the  notation  adopted  here.  We  use  ==^r  to  denote  weak 

p 

convergence  [5],  and  — >r  to  denote  convergence  in  probability  (with  r  going  to 
infinity).  We  write  f(x )  ~  g(x )  ( x  — »  oo)  when  lim  f(x)/g(x)  =  1.  Equality 

X^-OO 

in  distribution  is  denoted  by  =st,  inequality  in  the  strong,  convex  and  increasing 
convex  stochastic  ordering  sense  are  denoted  by  <st,  <cx  and  <iCX,  respectively. 

1.3  Self— similarity,  stable  distributions  and  reg¬ 

ular  variation 

This  section  provides  a  quick  tour  into  some  background  concepts  which  recur 
throughout  the  dissertation.  The  material  presented  here,  and  much  more,  can  be 
found  in  [3,  7,  16,  20,  53].  We  start  with  a  definition  of  long-range  dependence: 

Definition  1.3.1  The  R -valued  wide  sense  stationary  process  {Yk,  k  =  0,  ±1, . . .} 
is  said  to  be  long-range  dependent  if 

OO 

^  |cov  [Yk,  Ip]  |  =  oo  (1.1) 

fc=i 

and  short-range  dependent  otherwise. 
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We  proceed  with  self-similarity.  Roughly  speaking,  a  structure  is  “self-similar” 
if  it  appears  the  same  on  any  scale,  large  or  small.  The  term  self-similar  and  the 
definition  below  are  due  to  Mandelbrot  [41]. 

Definition  1.3.2  The  1R -valued  process  {X(t),  t  G  1R}  is  (strictly)  self-similar 
with  index  (or  Hurst  parameter)  H  >  0  if  for  all  a  >  0  the  finite- dimensional 
distributions  of{X(at),  t  G  R}  are  identical  to  the  finite- dimensional  distributions 
of  {aH  X(t),  t  G  R},  i.e.,  if  for  any  n  =  1,2,...,  t1;  t2, . . .  ,tn  in  R,  and  a  >  0, 

(X(at1),X(at2),  •  •  .,X(atn))  =st  (aHX(ti),  aHX(t2), .  ..,aHX(tn)). 

From  this  definition  we  see  that,  in  the  context  of  stochastic  processes,  self¬ 
similarity  is  tantamount  to  scale  invariance  of  the  finite-dimensional  distributions, 
but  not  necessarily  of  the  sample  paths. 

Definition  1.3.3  We  say  that  the  R -valued  process  { X(t),  t  G  R}  is  H-sssi  if  it 
is  strictly  self-similar  with  index  H  >  0  and  has  stationary  increments. 

Among  the  II -sssi  processes  the  Gaussian  one  is  the  most  prominent;  this  is  in 
part  due  to  the  fact  that  it  has  been  widely  applied  in  the  context  of  long-range 
dependence. 

Definition  1.3.4  A  H  -sssi  Gaussian  process  with  index  0  <  H  <  1  is  called 
fractional  Brownian  motion  and  is  denoted  by  { Bn(t ),  t  G  R}.  It  is  called  standard 
fractional  Brownian  motion  if  yax[BH{l)\  =  1. 

Proposition  1.3.1  With  0  <  H  <  1,  the  fractional  Brownian  motion  {B nit),  t  G 
R}  has  the  following  properties: 

(a)  Bh( 0)  =  0  a.s. 
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(b)  Its  covariance  function  is  given  by 


cav[BH(ti),BH(t2)}  =  ^  {\ti\2H  +  \t2\2H  -  \ti  -  t2\2H}  var[£#(l)],  tx,t2  G  R. 

(c)  When  0  <  H  <  1,  we  have  E  [Bait)}  =  0  for  all  t  in  R. 

(d)  When  H  =  1,  we  have  B1(t)  =  tB i(l)  a.s.  for  all  t  in  R. 


In  the  case  H  =  1/2  Proposition  1.3.1(b)  reads 


cov[£>i/2(ti),  Bi/2(t2) 


var  [-Bi/2(1)]  min(£i,  t2)  if  txt2  >  0 
0  if  t\t2  <  0 


and  {£1/2  (t),  t  G  R}  is  classical  Brownian  motion. 

Consider  now  the  increment  process  (Z^(n),  n  =  0,  ±1, . . .}  associated  with 
{BH(t),  t  G  R}  and  defined  by 


ZH(n)  :=  BH(n  +  1)  -  BH(n),  n  =  0,±1,....  (1.2) 


This  Gaussian  sequence  is  stationary,  because  fractional  Brownian  motion  {£#(£),  t  G 
R}  has  stationary  increments. 


Definition  1.3.5  The  stationary  process  {Z H{n) ,  n  =  0,±1, . . .}  of  (1.2)  is  called 
fractional  Gaussian  noise.  It  is  called  standard  fractional  Gaussian  noise  if 
var[ZH(l)]  =  1. 


From  Definition  1.3.5  and  Proposition  1.3.1(b)  it  follows  that  the  covariance  func¬ 
tion  rH[n )  :=  cov  [ZH(n  +  1),  ZH(1)\  of  the  fractional  Gaussian  noise  process 
{Zn{n),  n  =  0,  ±1, . . .}  is  given  by 

ruin)  =  ^  [\n  +  1\2H  —  2\n\2H  +  \n  —  1\2H]  var  [ZH{1)\ ,  n  =  0,  ±1, _  (1.3) 
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If  H  =  1/2  then  r^(n)  =  0  for  n  ^  0,  in  which  case  {Zi/2(n),  n  =  0,±1,...} 
forms  a  sequence  of  i.i.d.  Gaussian  rvs  and  fractional  Gaussian  noise  reduces  to 
the  familiar  white  Gaussian  noise.  For  H  ^  1/2  it  follows  from  (1.3)  that 

rH(n)  ~  var  [ZH(1)\  H(2H  —  1)  n2H~2  (n  — »  oo),  (1.4) 

so  that  ruin)  behaves  like  a  power  function.  Note  that  lim  ruin)  =  0  for  all 

n— >-oo 

0  <  H  <  1,  however,  for  1/2  <  II  <  1  the  covariance  function  rjj  (n)  decays  so 
slowly  as  n  — »  oo  that  the  corresponding  sum  (1.1)  diverges.  Thus,  in  the  case 
1/2  <  H  <  1  the  fractional  Gaussian  noise  process  {Zn(n),  n  =  0,±1,...}  is 
long-range  dependent. 

We  next  present  another  class  of  processes  which  also  contains  the  classical 
Brownian  motion  as  a  special  case.  First,  recall  the  family  of  stable  distributions. 

Definition  1.3.6  A  non- degenerate  rv  X  is  said  to  have  a  stable  distribution  if 
for  any  positive  numbers  c\  and  c2  there  is  a  positive  number  a(ci,c2)  and  a  real 
number  b(ci ,  c2)  such  that 

C\X\  +  c2X2  =st  a(ci,  c2)X  +  b(c\,  c2), 

where  Xi  and  X2  are  i.i.d  copies  of  X. 

The  characteristic  functions  of  all  stable  distributions  were  discovered  by  Levy 
(1924): 

Proposition  1.3.2  A  rv  X  has  a  stable  distribution  if  and  only  if  its  characteristic 
function  is  of  the  form 

E  [exp(i0X)]  =  exp  {i/id  —  8a\9\a  (1  —  i/3sgn(9)  z(9,  a))}  ,  (1.5) 
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where  p  is  a  real  constant,  8  >  0,  a  in  (0,  2],  (3  in  [—1, 1]  and 

[  tariff)  if  o+l, 

z(0,a )  :=  <  2  ^  2 

- In  \6\  if  a  —  1. 

v  7 r 

Stable  distributions  are  the  only  possible  limit  distributions  for  normalized  and 
centered  sums  of  i.i.d.  rvs,  and  in  that  sense  they  generalize  the  Gaussian  distri¬ 
bution,  which  is  obtained  by  setting  a  =  2  in  (1.5).  When  0  <  a  <  2  methods  for 
generating  deviates  from  stable  laws  are  available  but,  with  a  few  exceptions,  closed 
forms  expressions  for  stable  densities  are  not  known.  However,  series  expansions 
and  the  tail  behavior  of  stable  distributions  are  known. 

Let  Sa(8,/3,fi)  denote  the  generic  stable  rv  distributed  according  to  (1.5).  The 
next  two  results  can  be  found  in  the  monograph  by  Samorodnitsky  and  Taqqu  [53, 

pp.  16,  18]. 

Proposition  1.3.3  If  X  —st  Sa(8,/3,n )  with  0  <  a  <  2,  then 


lim  xaP  [X>x]  =  K(^±l  8° 

rr->-+oo  2 


(1,6) 


and 


where  Kn  := 


lim  xQP  [X  <  - x \  =  Ka- — -  (5° 

£—>•+00  2 


POO  \  -1 

x~a  sin  x  dx 


As  a  consequence,  we  have 


Proposition  1.3.4  If  X  =st  Sa(8,/3,fi)  with  0  <  a  <  2,  then 


E  [|A+]  <  cx)  for  0  <  p  <  a, 


and 


(1,7) 


E  [|W|P]  =  cx)  for  p>  a. 
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In  particular,  for  1  <  a  <  2  it  holds  that 


var  [X]  =  oo  and  E  [|X|]  <  oo. 

Because  of  this  infinite  pth  moment  property  for  p  >  a,  stable  distributions  are 
candidates  for  modeling  phenomena  with  high  variability. 

Let  us  now  consider  the  following  stochastic  process: 

Definition  1.3.7  With  a  in  (0,2],  the  R -valued  process  { La(t ),  t  >  0}  is  called 
(standard)  a-stable  Levy  motion  if 

(a)  La(0)  =  0  a.s., 

(b)  {La(t),  t  >  0}  has  independent  increments,  and 

(c)  La(t)  —  La(s)  — st  Sa((t  —  s)l^a,(3,0)  for  all  0  <  s  <  t  <  oo,  for  some 
(3  in  [—1, 1] . 

Clearly,  the  2-stable  motion  t  >  0}  is  simply  the  Brownian  motion.  More¬ 

over,  using  (1.5)  it  can  be  verified  that  for  all  c  >  0  the  processes  {La(ct),  t  >  0} 
and  {c1/" La(t),  t  >  0}  have  the  same  finite-dimensional  distributions,  and  an 
a-stable  Levy  motion  { La(t ),  t  >  0}  is  self-similar  with  index  H  =  1/ a  (unless 
a  =  1  and  (3  ^  0). 

In  view  of  the  power-like  tail  behaviors  encountered  in  (1.4),  (1.6)  and  (1.7),  it 
is  appropriate  to  review  the  notion  of  regular  variation.  The  definition  below  can 
be  interpreted  as  introducing  a  class  of  generalized  power  functions: 


Definition  1.3.8  A  Lebesgue  measurable  function  f  :  R+  — *  R+  is  said  to  be 
regularly  varying  (at  infinity)  with  index  p  in  R  if 


f(xy) 


lim  =  yp, 


z-y+oo  f(x) 

If  p  =  0  in  ( 1.8 )  then  f  is  called  slowly  varying. 


y  >  o. 


(1.8) 
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From  (1.8)  it  follows  that  if  /  :  R+  — »  R+  is  regularly  varying  with  index  p  then 
it  can  be  written  as 

f(x)  =  xp  h(x),  x  >  0, 


where  the  function  h  :  R+  — >  R+  is  slowly  varying. 

Work  on  functions  of  regular  variation  was  initiated  by  Karamata  (1930).  Later, 
regular  variation  and  in  particular  its  relevance  in  probability  was  popularized  by 
Feller  [20];  an  authoritative  treatment  of  the  subject  is  found  in  the  monograph 
[7].  The  following  theorem  due  to  Lamperti  [35]  (see  also  [7,  p.  356],  [3,  p.  50]) 
provides  a  connection  between  regular  variation  and  self-similarity:  Self-similar 
processes  are  exactly  those  that  arise  from  probability  limit  theorems  where  a 
process  is  centered  and  rescaled. 


We  use 


f.d.d. 


to  denote  convergence  of  finite-dimensional  distributions. 


Theorem  1.3.1  Suppose  that  the  IR.-valued  process  {X(t),  t  G  1R.}  is  such  that, 
with  suitably  chosen  mappings  /,  g  :  R  — *  R 

[X{Ut)f{u)9{U\  t  e  b}  = #  {Y(t),  t  g  R}  (u  ^  oo)  (1.9) 

for  some  R -valued  process  (T(t),  t  G  R}  with  non- degenerate  T(l).  Then  {Y(t), 
t  G  R}  is  strictly  self-similar,  and  all  self-similar  processes  arise  in  this  way. 
Moreover,  f  is  regularly  varying  with  index  H ,  where  H  is  the  Hurst  parameter  of 
the  limiting  self-similar  process  {Y(t),  t  G  R}. 

The  fact  that  the  norming  functions  /  appearing  in  (1.9)  are  necessarily  regularly 
varying  shows  that  regular  variation  is  intrinsically  present  in  limit  theorems  of 
probability. 
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Chapter  2 


M\G\oo  and  related  models 


In  this  chapter  we  introduce  a  class  of  traffic  models  based  on  the  M\G\oo  busy 
server  process.  In  addition,  we  discuss  simpler  related  queueing  models,  with  which 
stochastic  comparisons  will  be  sought. 

2.1  The  M\G\  oo  arrival  processes 

We  start  by  presenting  the  M\G\oo  arrival  processes,  together  with  the  assumptions 
and  notation  that  will  be  used  throughout.  Several  key  properties  concerning  this 
class  of  processes  are  stated  here  without  proof;  additional  details  are  given  in 
[11,  47], 

2.1.1  Definitions  and  basic  properties 

Consider  a  population  of  infinitely  many  information  sources,  operating  in  discrete¬ 
time.  Sources  can  be  in  one  of  two  states,  active  or  idle.  During  time  slot  [n,  n  + 
1),  n  =  0,1,...,  f3n+ 1  new  sources  become  active.  Source  j,  j  =  1, . . . ,  (3n+i, 
begins  generating  information  by  the  start  of  slot  [n  +  1,  n  +  2),  its  activity  period 
has  duration  crn+ij  (in  number  of  slots).  While  active,  each  source  generates 
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information  at  a  constant  rate  of  one  information  unit  (packet)  per  time  slot.  After 
its  activity  period  expires,  each  source  switches  off  permanently,  never  to  generate 
packets  again.  Let  bn  denote  the  number  of  active  sources,  or  equivalently,  the 
number  of  packets  generated  by  the  active  sources  at  the  beginning  of  time  slot 
[n,n  +  1).  If  initially  (i.e. ,  at  time  n  =  0)  there  were  already  b  active  sources, 
we  denote  by  a0j  the  residual  activity  duration  (in  time  slots)  for  the  jth  active 
source,  j  =  1, ...  ,b. 

For  each  n  —  0, 1, ...  we  have  the  decomposition 

bn  =  6l0)  +  <£>  (2.1) 

where  the  rvs  b\^  and  bffl  describe  the  contributions  to  the  number  of  active  sources 
at  the  beginning  of  slot  [n,  n  +  1)  from  the  sources  already  active  at  n  =  0  and 
from  subsequent  activations  in  slots  [k,  k  +  1),  k  =  1,2  respectively.  We 

readily  check  that 

b  n  (3k 

bn]  =  1  Kj  >  ^  aIld  1  ^  >n~k}-  (2-2) 

j= 1  k= 1  j= 1 

The  rv  b^  can  also  be  interpreted  as  the  number  of  active  sources  at  the  begin¬ 
ning  of  slot  [n,  n  +  1)  given  that  all  sources  were  silent  at  time  n  =  0.  On  the 
other  hand,  to  obtain  a  stationary  process  {bn,  n  —  0, 1, . . .},  the  rv  b^}  should 
be  specified  as  the  number  of  active  sources  at  time  n  =  0  given  that  the  sources 
have  been  operating  since  time  n  =  —  oo.  This  requirement  dictates  the  appropri¬ 
ate  distributional  assumptions  on  the  rvs  b  and  {a 0j-,  j  =  1,2, . . .}.  With  these 
considerations  in  mind,  we  now  record  a  set  of  assumptions  enforced  throughout. 

Assumption  (A)  The  IN -valued  rvs  b,  {/3n+ 1,  n  =  0, 1, . . {<Jn,j,  n  =  1,2,...; 
j  =  1,  2, . . .}  and  {a0j,  j  =  1,2,.. .}  are  defined  on  a  common  probability  space 
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{Cl,  IF.  P)  and  satisfy  the  following: 

(i)  These  rvs  are  mutually  independent. 

(ii)  The  rv  h  is  a  Poisson  rv  with  parameter  AE  [a], 

(in)  The  rvs  {/3n+i,  n  —  0, 1, . . .}  are  i.i.d.  Poisson  rvs  with  parameter  A  >  0. 

(iv)  The  rvs  {a n  =  1, . . . ;  j  =  1,2,...}  are  i.i.d.  with  distribution  function 
G  on  {1,  2, . . Let  a  denote  the  generic  IN -valued  rv  distributed  according  to  G. 
We  assume  that  E  [a]  <  oo. 

(v)  The  rvs  {<7o,j,  j  =  1,2,...}  are  i.i.d.  IN -valued  rvs  distributed  according 
to  the  forward  recurrence  time  distribution  G  associated  with  G,  i.e.,  if  a  denotes 
a  generic  IN-valued  rv  distributed  according  to  G,  then 

P  \(j  79] 

gn:=P[a  =  n}  =  ^~nJ,  n=  1,2,....  (2.3) 

E  [<j\ 

The  proposition  below  summarizes  the  properties  of  the  resulting  process  {bn,  n  = 
0, 1, . . .}  and  is  a  consequence  of  Assumption  (A),  (2.1)  and  (2.2)  [47]. 

Proposition  2.1.1  The  process  {bn,  n  =  0, 1, . . .}  is  a  (strictly)  stationary  ergodic 
process  with  the  following  properties: 

(a)  For  each  n  —  0, 1, . . .,  the  rv  bn  is  a  Poisson  rv  with  parameter  AE  [a]; 

(b)  Its  covariance  function  is  given  by 

cov  [bn+j,  bn]  =  AE  [((7  -  j)+]  =  AE  [a]  P  [a  >  j] ,  n,j  =  0, 1, ... ; 

(c)  Its  index  of  dispersion  of  counts  (IDC)  is  given  by 

OO  OO  . 

IDC  :=  J]  cov  [bn+j,  bn]  =  AE  [a]  ^  P  [a  >  j\  =  ^E  [a(a  +  1)]  (2.4) 

3= 0  3= 0 

and  the  process  is  short-range  dependent  (i.e.,  IDC  finite)  if  and  only  if  E  [a2]  is 
finite. 
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From  part  (b)  above  it  is  clear  that  the  sequence  { bn ,  n  =  0, 1, . . .}  exhibits  some 
form  of  positive  dependence.  In  fact,  as  mentioned  in  [45],  the  rvs  {bn,  n  —  0, 1, . . .} 
are  strongly  positively  correlated,  in  a  sense  that  can  be  made  precise  by  using  the 
notion  of  association  [19]: 

Proposition  2.1.2  The  rvs  { bn ,  n  =  0, 1, . . .}  are  associated,  in  that  for  any 
n  =  0, 1, . . .  and  any  pair  of  non-decreasing  mappings  /,  g  :  lNn+1  — *  R  we  have 

E  [/(6o,  •  •  • ,  bn)g(b0 , . .  • ,  bn)}  >  E  [f(b0, . . . ,  bn)]  E  [g(b0, . . . ,  bn)}  (2.5) 
provided  the  expectations  exist  and  are  finite. 

In  summary,  the  process  {bn,  n  =  0, 1, . . .}  results  from  discrete-time  Pois¬ 
son  arrivals  of  information  sessions,  where  the  generic  session  duration  rv  a  is 
distributed  according  to  the  pmf  G  and  the  packet  generation  rate  of  an  on¬ 
going  session  is  one  packet  per  time  slot.  It  is  fully  characterized  by  a  pair 
(A ,G),  with  A  the  Poisson  arrival  rate  (per  slot).  Under  Assumption  (A)  the 
sequence  { bn ,  n  =  0, 1, . . .}  can  be  identified  as  the  stationary  busy  server  pro¬ 
cess  of  a  discrete-time  M\G\oo  queue;  for  this  reason  the  packet  arrival  process 
{bn,  n  —  0, 1, . . .}  is  henceforth  referred  to  as  the  M\G\oo  input  process.  From 
Propositions  2.1.2  and  2.1.1(b)  we  see  that  { bn ,  n  =  0,1,...}  is  an  associated 
process,  whose  positive  correlation  structure  is  completely  determined  by  the  dis¬ 
tribution  of  a  (and  thus  of  cr).  In  many  cases  the  inverse  is  also  true,  i.e. ,  it  is 
possible  to  extract  M\G\oo  model  parameters  to  match  a  given  autocorrelation 
function. 

Proposition  2.1.3  An  R  +-valued  sequence  (</>(n),  n  —  0, 1, . . .}  is  the  autocor¬ 
relation  function  of  an  M\G\oo  process  (A ,  a)  if  and  only  if  the  mapping  n  — *  f>(n) 
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is  decreasing  and  integer-convex  with  0(0)  =  1  >  0(1)  and  lim  0(n)  =  0,  in  which 

n— >•  oo 

case  the  corresponding  distribution  of  a  is  given  by 


P  [a  >  n 


0(n )  —  0(n  +  1) 

i  -  0(1) 


n  =  0, 1, . . . 


Based  on  this  property,  M\G\oo  processes  have  been  used  to  model  VBR  video 
traffic  in  [34], 


2.1.2  Second  order  self— similarity 


In  [11]  Cox  observed  that  when  G  is  a  Pareto  distribution  with  parameter  a,  1  < 
a  <  2,  the  M\G\oo  busy  server  process  has  the  so-called  second  order  asymptotic 
self-similarity  property.  That  is,  its  correlation  structure  is  asymptotically  that 
of  the  increments  of  a  strictly  self-similar  process.  The  covariance  function  of 
standard  fractional  Gaussian  noise  is  given  by  (1.3).  Using  (1.3)  and  Proposition 
2.1.3  we  find  that  if  the  activity  rv  a  is  distributed  according  to 


P  \<7  >  n]  = 


In  —  1 


2  H 


3\n\2H  +  3|n  +  1 


2  H 


I  n  +  2 


2  H 


4(1  -  22H 


-2\ 


n  =  0,l,...,  (2.6) 


with  1/2  <  H  <  1,  then  the  corresponding  M\G |oo  input  process  has  the  same 
correlation  function  (1.3)  as  a  fractional  Gaussian  noise  process  with  Hurst  param¬ 
eter  H.  This  already  provides  a  point  of  contact  between  the  M\G\oo  process  and 
the  increments  of  a  strictly  self-similar  process.  Clearly,  the  particular  distribu¬ 
tion  (2.6)  achieving  this  match  is  too  restrictive.  It  is  relaxed  as  follows:  For  each 
m  —  1,2,...,  we  introduce  the  process  n  =  0, 1, . . .}  defined  by 

1  m—  1 

i-ir*  -  - y; fw,  n=i,2,...,  (2.7) 
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so  that  { b ^  ,n  =  0, 1, . . .}  is  also  a  stationary  process,  obtained  from  {bn,  n  = 
0, 1, . . .}  by  averaging  over  blocks  of  size  m.  Denote  its  covariance  function  by 

r‘m>(fc):=cov[6iy^],  *  =  0,1,.... 


We  say  that  the  original  process  {bn,  n  —  0, 1, . . .}  is  asymptotically  second-order 
self-similar  if  the  correlation  function  of  {b^n\  n  —  0, 1, . . .}  tends,  as  the  block 
size  m  goes  to  infinity,  to  the  correlation  function  of  fractional  Gaussian  noise,  i.e., 
if  for  each  lag  k  —  1, 2, . . .,  we  have 


r(m)  (k)  rH(k ) 

lim  — — —  = - 

m.->oof(m)  (0)  rn(  0) 

where  rn(k)  is  given  by  (1.3).  Noting  that 


(2.8) 


1 


rM(0)  =  _  var[60]  +  2^(1 - )  cov[6n, b0] 

777,  ^ ^  \  m  J 


n= 1 


(2.9) 


it  was  shown  in  [37]  that  the  M\G\oo  process  is  asymptotically  second-order  self¬ 
similar  with  parameter  H,  1/2  <  H  <  1,  if  the  tail  of  a  is  regularly  varying  with 
index  —(3  —  2 H),  i.e., 


P  [a  >  n]  =  n  ^  2H^h(n),  n  =  1,2,...,  (2-10) 


for  some  slowly  varying  function  h  :  1R+  — >  H+.  The  specihc  distribution  (2.6)  is 
simply  one  instance  of  (2.10).  Whenever  1/2  <  H  <  1,  the  tail  behavior  (2.10) 
implies  E  [a2]  =  oo,  and,  by  Proposition  2.1.1(b),  the  corresponding  second-order 
asymptotically  self-similar  M\G\oo  process  is  also  long-range  dependent. 

The  study  of  queueing  systems  in  heavy  traffic  typically  involves  operations 
such  as  accumulating  over  time  and  rescaling.  For  this  reason  we  expect  that  a 
heavy  traffic  analysis  of  a  queue  with  M\G\oo  arrivals  will  provide  a  natural  way  to 
further  explore  the  connection  with  self-similar  processes  discussed  in  this  section, 
and  address  this  problem  in  detail  in  Chapter  3. 
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2.2  The  queueing  system 


We  now  feed  the  M\G\oo  arrival  stream  {bn,  n  =  0, 1, . . .}  presented  in  Section 
2.1.1  into  a  discrete-time  single  server  queue  with  infinite  buffer  capacity.  Such 
a  queueing  system  routinely  serves  as  a  model  for  a  network  multiplexer:  If  qn 
denotes  the  number  of  packets  remaining  in  the  multiplexer  buffer  by  the  end  of 
slot  [n  —  1  ,n),  and  the  multiplexer  output  link  can  transmit  c  packets/slot,  then 
the  buffer  content  sequence  {qn,  n  =  0, 1, . . .}  evolves  according  to  the  Lindley 
recursion 

q0  =  q ;  qn+l  =  [qn  +  bn+1  -  c]  +  ,  n  =  0, 1, . . .  (2.11) 

for  some  initial  buffer  content  q  >  0.  To  identify  conditions  ensuring  existence  of 
a  finite  stationary  version  of  {qn,  n  =  0, 1, . . .}  and  to  determine  its  properties,  we 
appeal  to  established  results  on  recursions  of  the  form  (2.11).  We  introduce  the 
partial  sums  {sn,  n  =  0, 1, . . .}  defined  by 

n 

So  :=  0;  sn:=^2bj,  n  =  1,2,...  (2.12) 

3= 1 

and  specialize  the  results  from  [2,  54]  to  the  present  setup. 

Proposition  2.2.1  If  AE  [a]  <  c,  then  the  Lindley  recursion  (2.11)  is  termed 
stable  and  the  following  statements  hold: 

(a)  There  are  infinitely  many  n  such  that  qn  =  0. 

(b)  The  sequence  {qn,  n  =  0,1,...}  of  (2.11)  and  the  sequence  {g°,  n  = 
0, 1, . . .}  constructed  from  (2.11)  with  q0  =  0  strongly  couple,  in  the  sense  that 

Qn  =  Qn ,  n>m(q),  (2.13) 

where  m(q)  :=  min{n  —  0, 1, . . .  :  qn  =  0}  <  oo  (by  part  (a)). 

(c)  The  convergence  qn  =^n  qoo  takes  place,  where  the  stationary  lR+-valued 
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rv  q0 o  is  a.s.  finite  and  given  by 


<?oo  =st  sup{sn  -  nc]  n  =  0, 1, . . (2.14) 

Part  (a)  above  is  (1.2.5)  in  [2,  p.  71]  (or  Lemma  6.1(3)  in  [54])  and  (b)  is  implied 
by  Lemma  6.1(4)  in  [54],  Weak  convergence  in  (c)  is  a  consequence  of  (4.2.6)  and 
Remark  4.1.1  in  [2],  Expression  (2.14)  for  the  stationary  rv  follows  from  (2.2.3) 
of  [2], 

From  parts  (b)  and  (c)  of  Proposition  2.2.1  we  see  that  if  AE  [a]  <  c,  then  weak 
convergence  to  the  stationary  version  q^  takes  place  for  any  initial  condition  g; 
furthermore  the  distribution  of  q^  does  not  depend  on  q.  It  thus  suffices  to  restrict 
attention  to  the  choice  q  =  0  and  we  implicitly  do  so  from  now  on  whenever  we 
refer  to  (2.11).  In  that  case  the  system  is  initially  empty  and  the  output  to  the 
Lindley  recursion  admits  an  equivalent  representation  given  by 

qo  =  0;  qn  =  sn  -  nc  -  inf  (sy  -  jc,  j  =  0, 1, . . . ,  n) ,  n  =  1,2,...  (2.15) 

where  the  partial  sums  {sn,  n  =  1,2,...}  are  defined  by  (2.12).  This  is  useful  for 
establishing  heavy  and  light  traffic  limit  theorems. 

2.3  Instantaneous  inputs 

By  “instantaneous”  inputs  we  refer  to  the  situation  where  each  arriving  session 
brings  all  of  its  workload  to  the  system  in  one  time  slot,  immediately  upon  arrival. 
These  inputs  are  to  be  contrasted  with  the  gradual  M\G\oo  inputs,  where  arriving 
work  is  spread  over  the  entire  duration  of  a  session.  Such  instantaneous  arrivals 
are  represented  by  the  LV-valued  sequence  of  i.i.d.  rvs  {un,  n  —  0, 1, . . .}  given  by 

fln+l 

Un+ 1  :=  y;  <Tn+i,i,  n  =  0,1,...,  (2.16) 

i= 1 
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where  the  families  of  i.i.d.  rvs  {/3n+i,  n  —  0, 1, . . .}  and  {crn+i,i,  n  =  0, 1, ,  i  = 
1,2,.. .}  are  as  in  Section  2.1.1.  These  arrivals  are  also  characterized  by  the  pair 
(A,  G )  and  we  use  u  to  denote  the  generic  rv  for  the  i.i.d.  sequence  {un,  n  = 

0,1,...}. 

We  offer  the  instantaneous  inputs  { un ,  n  =  0, 1, . . .}  to  the  same  multiplexer 
with  constant  release  rate  c.  Assuming  that  the  queue  is  initially  empty,  we  write 
the  corresponding  Lindley  recursion  for  the  queue  length  sequence  {q\f\  n  = 
0, 1, . . .}  as 

QoU)  =  0;  Qn+1  =  [Qn]  +  un+ 1  -  c]  +  ,  n  =  0, 1, ... .  (2.17) 

If  E 

M  =  AE  [a]  <  c  the  system  is  stable  and  the  convergence  q ^  q ^  takes 

place  for  some  1R+ -valued  rv  . 

Owing  to  the  independence  of  the  rvs  {un,  n  —  0, 1, . . .}  recursion  (2.17)  can, 
at  least  in  principle,  be  handled  by  standard  generating  function  techniques.  The 
details  of  this  approach  are  given  in  Section  2.3.1. 


2.3.1  A  Markov  chain  of  the  M|G|1  type 


Consider  a  Markov  chain  on  {0, 1, . . .},  whose  transition  probability  matrix  P  is 
of  the  M\G\1  type  and  is  given  by 


P  = 


fo  +  fl  $2  f% 

fo  fl  /2 

0  fo  fi 


(2.18) 


OO 

Assume  that  the  probability  vector  (/0,  f\, . . .)  satisfies  <  1,  in  which  case 

i= o 

the  Markov  chain  is  positive  recurrent.  Denote  by  (7T0,  tti,  . . .)  the  steady  state 
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probability  vector  associated  with  P,  and  set 

OO  OO 

II(z)  :=  'S^J'KiZl  and  F(z)  :  =  fjZ1,  z  G  D,  (2-19) 

i= 0  i=0 

where  I?  =  {sG(D:|s|<l}is  the  unit  disk  in  the  complex  plane  (D.  The  vector 
(7r0, 7r1; . . .)  satisfies 

TTo  =  (/o  +  /l)7To  +  /07Ti  (2.20) 


and 


k+ 1 

TTfc  =  fk-l+l^h  k=  1,2, ... . 
1=0 


We  multiply  (2.20)  by  2;  and,  for  each  fc  =  1,2,...,  the  kth  equation  in 
by  zk .  Adding  up,  invoking  definitions  (2.19)  and  manipulating  we  find 


(2.21) 

(2.21) 


zU(z)  =  {z-  l)/o7T0  +  (/o  +  fiz)U(z)  +  ( F(z )  -  fo  -  fiz)U(z), 


so  that 

nM  =  z  e  <2'22) 

We  compute  the  limit  as  z  — >  1  by  applying  F  Hospital’s  rule  on  the  right-hand 
side  in  relation  (2.22).  Since  n(l)  =  1,  we  get 


_  1  -  F'(  1) 


(2.23) 


and  inserting  (2.23)  back  in  (2.22)  yields 


n(z) 


(l-z)(l-F'(l)) 
F{z)  -  z 


z  e  D. 


(2.24) 


OO 

To  calculate  n'(l)  =  we  differentiate  (2.24)  and  apply  F  Hospital’s  rule 

i= 0 

twice,  so  as  to  conclude  that 


n'(i) 


i) 

2(1  -F'{1)Y 


(2.25) 
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2.3.2  Case  c  =  1:  A  solution  by  generating  function 

Set 

Q{u\z)  :  =  E  \zq~]  and  U(z)  :=  E  [zu]  ,  z  E  D.  (2.26) 

When  the  multiplexer  release  rate  is  c  =  1  the  sequence  {q^\  n  —  0, 1, . . .}  is 
a  Markov  chain  on  {0, 1, . . With  the  notation  of  Section  2.3.1,  its  transition 
matrix  is  of  the  form  (2.18)  with 

fr.=  P[u  =  i],  i  =  0,1,....  (2.27) 

Under  (2.27),  we  have  the  identification 

F(z)  =  U(z),  n(z)  =  Q(u\z),  zeD  (2.28) 

and 

f'(l)  =  E[«],  n'(l)  =  E  [<&>]  .  (2.29) 

As  we  plan  to  specialize  the  results  of  Section  2.3.1  to  the  inputs  {un,  n  =  0, 1, . . .}, 
given  by  (2.16),  we  note  that 

U(z)  =  exp  (A(E  [U7]  -  1))  and  U"{  1)  =  AE  [a(a  -  1)]  +  A2E  [a]2  .  (2.30) 

Thus,  because  of  (2.28)  and  (2.29),  relations  (2.23),  (2.24)  and  (2.25)  imply 


P[<A  =  0]  =  (l-AE[<r])eA 

(2.31) 

q(u\z)  =  (1 7 d)1  7 A1? (f1*  ■  2  6  D 

exp  (A(E  [za]  —  1 ))  —  z 

(2.32) 

prj.11  A(AEH2  +  EK«T-1)]) 

^°°J  2(1  —  AE  [cr]) 

(2.33) 
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2.3.3  Case  c  =  1:  An  equivalent  representation 

We  now  derive  a  representation  of  the  stationary  queue  size  q ^  in  terms  of  the 
forward  recurrence  times  of  the  input  sequence  { un ,  n  =  1,2,.. To  do  this  we 
introduce  the  sequence  of  i.i.d.  rvs  {un,  n  =  1,2, . . .}  with  generic  rv  u  whose 
distribution  is  given  by 

P  [u  =  0]  =  0;  P  [u  —  n]  —  \  P  [u  >  n] ,  77-  =  1,2, _  (2.34) 

E  [u\ 

The  corresponding  generating  function  is  given  by 

OO 

U(z)  :=  J^P  [u  =  n\zn 

n=  1 

00 

L  J  n—  1 

1  OO  OO 

L  1  n=  1  k=n 
1  _°°_  ~k+ 1  _ 

-,  z  <E  D. 


EM 


EM 


E 


fc=i 

^  1  -  U(z) 

EM  1  —  z 


(2.35) 


Under  the  stability  condition  E  [u]  <  1,  relations  (2.24)  (with  the  identification 
(2.28)  and  (2.29))  read 

(l-M(i-EM) 


Q{u\z)  = 


U(z)  -  z 


=  (1  -  E[u])  |  1  —  1  U{Z))  ,  z£D.  (2.36) 


1-z 


Therefore,  with  the  help  of  (2.35),  we  obtain 

1  -  E  M 


Q(u)(z )  = 


1  -  E  M  E  T 


t-11  ’  z  £  D, 


(2.37) 
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where  we  note  that  the  rv  u—  1  is  non-negative  because  of  (2.34).  Upon  rewriting 
(2.37)  as 

OO 

Q{u\z)  =  (1  -E[u])J^E[u]nE  [U2-1]”,  z  <E  D, 

71=0 


we  get  the  representation 

V 

gW=sf^(i-l),  (2.38) 

71=1 

where  the  rv  v  is  independent  of  {un,  n  =  1,2,...}  and  geometrically  distributed 
with  parameter  E  [it]  according  to 


P  [v  —  n]  =  E  [u]n  (1  —  E  [it]),  77.  =  0, 1, ,  (2.39) 


(with  the  convention  that  empty  sums  in  (2.38)  have  value  zero).  The  stationary 
queue  size  is  thus  expressed  as  a  geometric  sum  of  the  forward  recurrence  times 
{ un ,  n  =  1,2,...}  associated  with  the  input  sequence  {un,  n  =  1,2,.. .}. 


2.3.4  Case  c  =  1:  Idle  and  busy  periods 

We  now  take  one  more  look  at  recursion  (2.17);  this  will  prove  useful  in  the  de¬ 
velopments  of  Chapter  4.  We  view  the  queue  size  sequence  {q^\  n  —  0, 1, . . .}  as 
evolving  in  a  series  of  independent  regenerative  cycles,  alternating  between  zero 
and  positive  values.  If  the  queue  is  initially  empty  then,  for  each  n  =  1,2, . . ., 
the  nth  cycle  consists  of  an  idle  period  followed  by  a  busy  period,  with  respective 
lenghts  1^  and  B^\  (expressed  in  time  slots).  We  say  that  a  time  slot  is  part  of 
a  busy  period  if  the  queue  length  at  the  beginning  of  the  time  slot  is  positive.  If 
the  queue  length  at  the  left  slot  boundary  is  zero,  the  slot  is  considered  to  belong 
to  an  idle  period.  That  is,  for  each  n  —  1,2, . . .,  the  family  of  i.i.d.  pairs  of  rvs 
n  =  1,2,...}  associated  with  (2.17)  are  recursively  defined  by 

:=  mf{t  =  0,1,...  :  >  0}  (2.40) 
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and 


inf{t  =  0, 1, . . .  :  q 


(«) 


0}, 


(2.41) 


with  the  convention  that  empty  sums  are  zero. 

A  clarification  is  needed  as  the  terms  “idle”  and  “busy”  are  slightly  abused 
here.  For  example,  it  is  possible  that,  when  the  queue  length  is  zero  at  consecutive 
instants,  say  t  and  t+1,  this  occurs  because  of  a  single  arriving  packet  in  [t,  t  +  1) 
which  was  served  by  the  end  of  the  time  slot.  Such  a  slot  is  considered  to  belong 
to  an  “idle”  period,  despite  the  fact  that  the  server  was  busy  serving  the  arriving 
packet.  Thus  “idle”  and  “busy”  are  defined  here  in  reference  to  queue  content, 
and  not  to  server  activity.  For  lack  of  better  terminology,  we  shall  continue  to  use 
“idle”  and  “busy”  in  forthcoming  arguments,  referring  to  definitions  (2.40)  and 
(2.41)  to  resolve  any  confusion. 

Next,  denote  by  (I^u\  B^)  the  generic  idle  and  busy  period  pair  associated 
with  {(/<«>,  B<«>),  n  =  1,2,.. .}.  By  the  Renewal-Reward  Theorem  we  can  alter¬ 
natively  express  P  =  0]  (which  has  already  been  evaluated  in  (2.31))  as 

p  r  Ju)  _  o]  _ E  [/H] _  (2  42) 

[lo°  J  E  [JM]  +  E  [£(«)]  ’  1  j 

To  obtain  the  distribution  of  the  idle  period  1^  let 


77  :=  P  [ft  =  0]  T  P  \/3  =  1]  P  [a  =  1]  (2.43) 


denote  the  probability  that  at  most  one  unit  of  work  arrives  during  a  time  slot. 
When  c  =  1,  the  definitions  (2.40)  and  (2.41)  imply  that  I is  geometric  with 
parameter  77,  i.e., 

P  [/(u)  =  k]  =  P  [<21  =  ■  •  •  =  qk-i  =  0,  qk  >  0] 

=  k  =  1,2,...,  (2.44) 
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so  that 


E  [/(u)‘ 


1  —  T) 


(2.45) 


2.3.5  Case  c  =  1/m  (m  =  l,2, ...) 


The  arguments  presented  in  Section  2.3.2  can  be  easily  extended  to  address  the 
following  situation:  Fix  some  integer  m  =  1,2, .. .  and  consider  the  queue  length 
sequence  Ui”1,  n  =  0, 1, . . .}  resulting  from  the  recursion  (2.17)  with  multiplexer 
release  rate  c  =  1/m.  In  this  case  {q^\  n  =  0,1,...}  is  a  Markov  chain  on  the 
lattice  {0, 1/m,  2/m, . . .}.  The  transition  probability  matrix  is  again  of  the  form 
(2.18)  but  this  time,  instead  of  (2.27),  it  holds  that 


(2.46) 


IP  u  =  —  if  i  =  0  (mod  m) 

L  mj 

0  if  i  7^  0  (mod  m). 

In  place  of  (2.28)  we  now  have  the  identification 


F(z)  =  U{zm)  and  U(z)  =  Q{u) (zm) ,  z  E  D. 


(2.47) 


Consequently, 


F'(  1)  =  mE  [u] ,  F"(  1)  =  m2U"{  1)  +  m(m  -  1)E/'(1) 


(2.48) 


n'(l)  =  mE  [<&>]  . 


(2.49) 


As  the  analysis  leading  to  (2.23)  and  (2.25)  still  applies,  we  use  (2.48)  in  (2.23)  to 


obtain 


P  =  0]  =  (1  -  AmE  [u])eA, 


and  (2.30),  (2.48)  and  (2.49)  in  (2.25)  to  collect 

r  (u)]  =  A  (AmE  [a]2  +  mE  [a2]  —  E  [a]) 
Lg°°J  2(1  -  AmE  [cr]) 


(2.50) 


(2.51) 
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2.3.6  Poisson  inputs:  Stochastic  comparisons 

The  temporal  correlations  in  the  M |  G |  oc  arrival  process  are  expected  to  have  an 
adverse  effect  on  queueing  performance.  Insight  to  this  effect  can  be  obtained  from 
a  comparison  with  the  situation  where  these  correlations  are  altogether  eliminated, 
while  maintaining  the  same  Poisson  marginal  distribution.  To  that  end  we  consider 
the  companion  sequence  {q„\  n  =  0.1,...}  evolving  according  to 

Qo]  =  0;  Qnh  =  [qn  +  Cn+i  -  c]  +  ,  n  =  0, 1, . . .  (2.52) 


where  the  rvs  {£„,  n  =  1,2, . . .}  form  a  sequence  of  i.i.d.  Poisson  rvs  with  pa¬ 
rameter  AE  [cr] .  Clearly,  these  independent  Poisson  arrivals  fall  in  the  category  of 
the  instantaneous  inputs  of  (2.16),  where  the  pair  (A,  a)  is  replaced  by  (AE  [cr]  ,  1). 
Under  the  stability  condition  AE  [cr]  <  c,  convergence  to  the  stationary  ]R+-valued 
rv  takes  place.  Then,  for  c  =  1/m,  with  m  =  1,  2, . . .,  we  can  use  (2.50)  and 
(2.51)  to  obtain 

P  [g(0  =  0]  =  (1  -  AmE  [cr])eAE[crl  (2.53) 


and 


Noting  the  inequalities 


AE  [a]  (AmE  [cr]  +  m  —  1) 
2(1  —  AmE  [a]) 


(2.54) 


P  [«£’  =  0]  >  P  [<£>  =  0]  and  E  [««>]  <  E  [g«] 

we  suspect  that  q^  and  qj£  may,  at  least  in  some  circumstances,  act  as  stochastic 
lower  and  upper  bounds  respectively  to  the  stationary  queue  size  q^  induced  by 
M\G\oo  arrivals.  As  a  first  step  in  this  direction  we  now  establish  a  comparison 
between  q^  and  q^\  in  the  increasing  convex  stochastic  ordering  sense. 
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Proposition  2.3.1  If  AE  [<x]  <  c  and  E  [cr2]  <  oo,  then  the  stationary  queue 
lengths  and  qjff)  associated  with  recursions  (2.52)  and  (2.17)  satisfy 

q[S<icXq^.  (2.55) 

The  proof  of  Proposition  2.3.1  relies  on  the  following  fact: 

Lemma  2.3.1  For  Poisson  rvs  XlX  and  Xx  with  parameters  7A  and  X,  respec¬ 
tively,  it  holds  that 

XlX  <cx  ^Xx,  7  >  1.  (2.56) 


Proof.  We  first  prove  a  similar  comparison  result  for  certain  Bernoulli  rvs;  these 
are  subsequently  used  to  construct  the  Poisson  rvs  of  interest.  Let  denote  a 
generic  Bernoulli  rv  with  parameter  p,  0  <  p  <  1,  i.e., 

p  [W{p)  =  l]=p=l-P  [W{p)  =  0]  . 

Fix  some  integer  n  >  7A  and  consider  two  sequences  of  i.i.d.  Bernoulli  rvs 
|  w\ ^ ,  i  =  1, 2, . . . ,  nj  and  |  W-  n  \  i  =  1,  2, . . . ,  n|  with  generic  rvs  and 

W^\  respectively.  For  all  convex  mappings  0  :  H  — »  H  and  7  >  1,  it  holds  that 

0(7)  -  0(0)  -  7(0(1)  -  0(0))  >  0. 


Therefore,  upon  comparing 


with 


E 


0(W^?))|  =  ^  0(1)+  (  1- 


7A 


n 


7A 


n 


m 


-  9(7)  +  (1  -  1)  9(0), 

n  \  nj 
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we  get 


E 


cf)  (w^ 


<  E 


0  (VfP(> 


Thus,  by  Definition  B.4  (see  Appendix  B),  we  arrive  at 


<cx  7 W& 


(2.57) 


which  in  turn  [59,  Proposition  1.1.2]  implies  that 


£  (2.68) 

i= 1  i=  1 

By  Poisson’s  Convergence  Theorem  the  rvs  X a7  and  "jX\  can  be  obtained  as  the 

n  n 

Wi  "  and  y,  bbbj  "  ,  respectively,  by  letting  n  go  to 

i= 1  i=l 

infinity.  The  expectations  of  both  left  and  right  hand  side  of  (2.58)  are  finite  and 
equal  to  yA,  the  common  value  of  the  expectations  of  the  limiting  rvs.  Therefore 
Proposition  B.8  applies  and  (2.56)  follows  by  taking  the  limit  in  (2.58).  ■ 


Proof  of  Proposition  2.3.1.  Recall  that  the  generic  rvs  (3  and  £  are  Pois¬ 
son  rvs  with  parameters  A  and  AE[cr],  respectively,  with  E  [cr]  >  1.  Applying 
Lemma  2.3.1,  we  obtain 

£  <cz  p  E  [a]  (2.59) 

or,  equivalently,  by  Definition  B.4, 

E  [</>(£)]  <  E  [</>(/?  E  [cr])]  (2.60) 


for  all  convex  mappings  (j)  :  H  — »  1R  for  which  the  expectations  exist.  On  the  other 
hand,  Jensen’s  inequality  implies 


</> 


-  p 


Lfc=i 


1/3 
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whence 


E[#E[a])]<E[^)].  (2.61) 

Combining  (2.61)  and  (2.60)  we  collect  E  [</>(£)]  <  E  [<j){u) ],  which  (again  by  Defi¬ 
nition  B.4)  is  tantamount  to 

£  <cx  U.  (2.62) 

Finally,  we  use  (2.62),  together  with  the  fact  that  E  [q^]  and  E  [qj£]  exist  and 
are  both  finite  when  E  [ex2]  <  oo,  to  conclude  that  (2.55)  holds  true  by  appeal¬ 
ing  to  the  external  monotonicity  of  GI\GI\  1  recursions  given  by  Proposition  B.8.  ■ 

In  the  case  c  =  1  the  rv  q (f)  (and  qj£)  as  well)  admits  the  equivalent  repre¬ 
sentation  (2.38)  given  in  Section  2.3.3.  This  enables  us  to  sharpen  the  result  of 
Proposition  2.3.1  as  follows: 

Proposition  2.3.2  Let  c  =  1  in  the  recursions  (2.17)  and  (2.52).  If  XE[a]  <  1, 
then  the  stationary  queue  lengths  q and  qj£)  satisfy 

V{£  <st  (2-63) 


Proof.  When  c  =  1  relation  (2.38)  for  is  in  effect,  and  a  corresponding 
expression  holds  true  for  qj£) .  That  is, 

=.i  E(6.  -  1),  (2-64) 

n= 1 

where  the  generic  rv  £  for  the  sequence  of  i.i.d.  rvs  {£„  ,  n  =  1,  2, . . .}  is  distributed 
according  to 


£  =  0 


=  0;  P 


£  =  n 


E[£] 


P[£>n],  to  =  1,2,. 


(2.65) 
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and  the  rv  /j  is  independent  of  {£n,  n  =  1,2,.. with 


P  =  n\  =  E  [£]"'  (1  -  E  [£]),  n  =  0, 1, . . . .  (2.66) 


Noting  that  E  [£]  =  E  [it]  we  have  fi  =st  /y,  where  the  distribution  of  v  was  given 
in  (2.39)  .  Also,  the  convex  stochastic  comparison  (2.62)  implies 

E  [(£  -  ™)+]  <  E  [{u  -  n)+]  ,  n  =  0, 1, . . . .  (2.67) 


From  (2.34)  and  (2.65)  we  see  that 

1 


f  >  n 


EB] 


E[({-n)+]  r»  =  0,l, 


with  the  corresponding  relation  for  u.  Inequality  (2.67)  now  yields 


f  >  n 


<  P  \u  >  n]  ,  n  —  0,1, 


or,  equivalently, 


£  <st  u. 


(2.68) 


Thus,  the  convex  stochastic  comparison  between  £  and  u  translates  into  a  strong 
stochastic  comparison  between  £  and  u,  and  the  conclusion  (2.63)  follows  from 
li  =st  v,  (2.38),  (2.64),  (2.68)  and  Proposition  2.2.5  in  [59,  p.  45].  ■ 


33 


Chapter  3 


Heavy  traffic:  Levy  motion  limits 

3.1  Introduction 

In  this  chapter  we  derive  the  non-degenerate  limiting  distribution,  as  the  traf¬ 
fic  intensity  AE  [cr]  tends  to  the  multiplexer  release  rate  c,  of  the  appropriately 
normalized  queue  length  induced  by  an  M\G\oo  arrival  process,  for  a  generally 
distributed  activity  rv  cr.  The  arising  limits  are  classified  in  terms  of  the  short- 
vs.  long-range  dependent  property  of  the  M\G\oo  process,  as  determined  by  the 
tail  behavior  of  a.  In  the  short-range  dependent  regime  the  limiting  distribution 
distribution  is  exponential,  as  is  the  case  in  the  classical  GI\G\1  queue,  originally 
studied  by  Kingman  in  [31,  32],  However,  under  long-range  dependence  the  results 
do  not  involve  the  fractional  Brownian  motion  model  of  Norros  [43,  44],  Different 
self-similar  limits  arise  in  the  form  of  Levy  motion,  leading  to  a  buffer  content 
distribution  with  hyperbolic  decay. 

The  basic  idea  behind  the  proof  of  these  results  is  a  “convergence  together” 
argument  which  allows  us  to  identify  processes  with  well-known  heavy  traffic  be¬ 
havior,  under  both  short-  and  long-range  dependence.  This  is  accomplished  chiefly 
by  combining  standard  results  on  stable  rvs  and  their  domain  of  attraction  [20], 
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with  a  general  functional  convergence  result  for  processes  with  stationary  indepen¬ 
dent  increments  due  to  Skorokhod  [57].  We  point  out  that,  even  in  the  short-range 
dependent  case,  convergence  to  Brownian  motion  does  not  appear  to  follow  from 
standard  results  for  stationary  processes  [5,  Thm.  20.1,  p.  174],  as  it  is  not  ob¬ 
vious  that  the  M\G\oo  busy  server  process  satisfies  the  required  mixing  property. 
However,  as  mentioned  in  Proposition  2.1.2,  the  M\G\oo  busy  server  process  is 
strongly  positively  correlated  -  it  is  an  associated  process.  Because  of  this  prop¬ 
erty,  it  is  then  possible  under  short-range  dependence  to  develop  an  alternative 
approach  similar  to  that  used  by  Newman  and  Wright  in  [42]  in  establishing  the 
Invariance  Principle  for  sequences  of  associated  random  variables.  This  approach 
is  not  pursued  here. 

3.2  The  heavy  traffic  regime 

We  seek  to  understand  the  behavior  of  the  (stable)  queue  with  the  the  correlated 
M\G\oo  arrival  process,  under  the  assumption  that  it  is  almost  fully  utilized,  i.e., 
AE  [cr] ,  though  less  than  the  release  rate  c,  is  very  close  to  c.  This  typically  involves 
obtaining  limiting  expressions  of  properly  rescaled  quantities  of  interest,  as  the 
packet  arrival  rate  AE  [cr]  tends  towards  its  critical  value  c.  Here,  the  quantity 
of  interest  is  the  steady-state  queue  size  q^.  A  natural  setup  to  investigate  this 
problem  consists  of  embedding  the  discrete-time  queue  with  release  rate  c  driven  by 
an  M\G\oo  input  process  (A,  cr)  into  a  parametric  family  of  like  queueing  systems, 
indexed  by  an  integer  parameter,  say  r.  More  precisely,  for  each  r  =  1,2,... 
we  take  the  rth  system  to  be  a  discrete-time  queue  with  release  rate  c  driven 
by  an  M\G\oo  input  process  {brn,  n  —  0, 1, . . .}  characterized  by  the  pair  (Ar,cr). 
The  corresponding  queue  size  sequence  {qrn,  n  =  0.1,...}  also  obeys  the  Lindley 
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recursion  (2.11),  and  admits  a  representation  of  the  form  (2.15),  i.e., 


=  qrn=  srn-nc-  inf  (s'  -  jc,  j  =  0,1,...,  n)  ,  n  =  1,2,... 


(3.1) 


where  {s”,  n  =  1,2,...}  is  the  sequence  of  partial  sums  (2.12)  associated  with 
{brn,  n  =  1,2,.. .}.  We  take  ArE  [cr]  <  c  for  all  r  =  1,  2, ...  for  some  fixed  c  >  0, 
so  that 


lim  Xr  =  c/E  [cr]  .  (3.2) 

r— too 

Thus,  each  one  of  these  systems  is  stable  with  lim  ArE  [cr]  =  c,  thereby  capturing 

r— >co 

the  notion  that  “the  system  is  driven  to  heavy  traffic.”  We  seek  a  scaling  sequence 
{Cr,  r  =  1,2,...}  such  that  the  convergence  in  distribution 


Cr 


Q 


(3.3) 


takes  place  to  some  R-valued  rv  Q. 

Unfortunately,  this  heavy  traffic  program  cannot  be  carried  out  in  this  form  as 
exact  expressions  are  unavailable  for  the  distribution  of  q ^  owing  to  the  correla¬ 
tions  present  in  the  M\G\oo  input  process,  and  we  need  to  resort  to  the  following 
indirect  approach  where  the  buffer  content  is  rescaled  in  both  the  time  and  state 
space  variables:  For  each  r  =  1,2,...,  we  define  the  R.-valued  continuous-time 
processes  {Sr(t)i  t  >  0}  and  {Qr(t),  t  >  0}  by 

Sr(t)  :=  }  (4,1  -  E  [sfrt]])  and  Q'(t)  :=  t  >  0, 

Sr  Sr 

and  the  function  7”  :  R+  — *  R  by 

Y(t)  ■=  7-  ([rt]c  -  E  [sr[tr]\)  =  M  (c  -  ArE  [a]) ,  t>  0. 

Sr  Sr 

The  convergence  (3.3)  can  be  stated  informally  as 


lim  lim  Qr(t)  =  Q  (3.4) 

cot— >00 
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with  limits  understood  in  the  sense  of  weak  convergence.  The  approach  to  heavy 
traffic  followed  here  is  to  interchange  the  order  of  these  limits,  i.e.,  to  evaluate 

lim  lim  Qr{t )  (3.5) 

t—¥ oo  r— >oo 

which  corresponds  to  first  taking  r  to  infinity,  and  then  letting  t  go  to  infinity. 
Assuming  that  the  limits  can  be  taken  in  that  order,  we  are  then  left  with  the  task 
of  showing  that 

lim  lim  Qr(t)  =  Q  =  lim  lim  Qr{t).  (3.6) 

r—>oo  t—>oo  t — yoo  r—>oc 

In  this  chapter  we  concentrate  only  on  establishing  the  first  step  (3.5),  and  it  is 

well  known  [27,  65]  that  the  theory  of  weak  convergence  on  function  spaces  provides 

a  natural  framework  for  doing  so.  To  that  end,  we  pause  briefly  to  introduce  the 

needed  notation,  as  well  as  to  highlight  several  points  from  the  theory  of  weak 

convergence  of  processes;  this  material  is  drawn  mostly  from  [5,  pp.  150-153]  to 

which  the  reader  is  referred  for  additional  information: 

For  each  T  >  0,  let  D[0,T]  denote  the  space  of  mappings  [0,  T]  — *  R  which  are 

right-continuous  with  left  limits;  the  space  D[0,T]  can  be  equipped  with  either 

the  uniform  topology  or  the  standard  Skorokhod  topology  [5,  p.  Ill],  As  in  [5,  p. 

150],  a  concept  prefixed  with  U  (resp.  S )  refers  to  the  uniform  (resp.  Skorokhod) 

topology.  For  probability  measures  defined  on  the  collection  of  I/-Borel  (resp.  S- 

Borel)  sets  on  D[0,T],  we  refer  to  weak  convergence  in  the  sense  of  the  uniform 

(resp.  Skorokhod)  topology  by  [/-weak  (resp.  S-weak)  convergence,  and  we  write 
u  s 

=^r  (resp.  =^r)  (with  the  understanding  that  r  goes  to  infinity).  For  probability 
measures  defined  on  the  collection  of  U-  Borel  sets,  U  weak  convergence  implies  S'- 
weak  convergence  but  the  converse  is  false.  This  implication  will  be  used  repeatedly 
in  various  technical  arguments  [Sections  3.5  and  3.7]. 


37 


Finally,  let  _D[0,oo)  denote  the  space  of  mappings  R+  — >  R  which  are  right- 
continuous  with  left  limits.  In  this  chapter,  we  present  results  on  the  S  weak 
convergence  of  the  restrictions  to  finite  intervals  of  sequences  of  R-valued  processes 
with  sample  paths  in  D[ 0,  oo).  More  precisely,  consider  the  sequence  of  R-valued 
processes  { Xr(t ),  t  >  0},  r  =  1,  2, . . .,  with  sample  paths  in  R[0,oo).  Whenever 
for  each  T  >  0  we  have  the  S'-weak  convergence 

{Xr(t),  0  <  t  <  T}  =^r  {X(t),  0  <  t  <  T}  in  D[0,T\ 

for  some  R-valued  process  {A(f),  t  >  0}  with  sample  paths  in  D\Q.  oo),  we  sim¬ 
plify  the  notation  by  writing 

{Xr(t),  t  >  0}^>r{X(t),  t  >  0}. 

Now,  noting  that  (3.1)  can  be  rewritten  as 

Qr(t)  =  Sr(t)-Y(t)~  inf  (S'(x)-Y(x)),  t>  0.  (3.7) 

0<x<t 

and  recalling  the  continuous  mapping  theorem  [5,  Thm.  5.1,  p.  30],  we  conclude 
that  the  first  limit  in  (3.5)  requires  at  the  very  least  identifying  a  scaling  sequence 
{£r,  r  =  1,2,...}  that  ensures  the  convergence 

{Sr(t),  t  >  0}  =^r  {S{t),  t  >  0}  (3.8) 

for  some  non-trivial  limiting  process  {S(t):  t  >  0}. 

3.3  The  main  heavy  traffic  results 

As  will  become  apparent  shortly,  the  choice  of  the  scaling  sequence  { (r ,  r  = 
1,2,...}  and  the  characterization  of  the  limiting  process  {S(t),  t  >  0}  entering 
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(3.8)  both  depend  on  the  distribution  of  the  rv  cr  which  controls  the  correlations 
in  the  input  packet  stream.  It  is  nevertheless  easy  to  see  that  in  order  to  avoid 
collecting  only  a  law  of  large  numbers,  any  candidate  scaling  sequence  { (r ,  r  = 
1,2,.. .}  should  obey  the  following  necessary  condition: 

Condition  (A)  The  scaling  sequence  {Cr,  r  =  1,2,...}  satisfies 

Cr 

lim  (r  =  +oo  with  lim  —  =  0. 

r—¥ oo  r—¥ oo  T 

The  heavy  traffic  assumption  below  rehnes  (3.2),  and  guarantees  that,  as  r  goes 
to  infinity,  the  family  of  queueing  systems  described  by  (3.7)  gradually  approaches 
instability  at  the  appropriate  speed: 


Assumption  (B) 

lim  (ArE  [a]  — 

r—¥  oo 


The  scaling  sequence  {Cr,  r  = 


c)  —  =  —7  or  equivalently , 
Sr 


1,2,.. .}  satisfies 
ArE  [cr]  =  C  -  -^(7  +  0(1)) 


for  some  7  >  0. 


Condition  (A)  and  Assumption  (B)  are  enforced  throughout.  It  is  worth  point¬ 
ing  out  that  the  scaling  sequence  {Cr,  r  =  1, 2, . . .}  is  essentially  unique,  i.e.,  any 
other  scaling  sequence  {Cr,  r  —  1,2, . . .}  yielding  a  non-degenerate  limit  in  (3.8) 
must  satisfy  lim  C(/Cr  =  C  for  some  finite  constant  C  >  0. 

r—¥  oo 

We  begin  with  the  case  where  the  M\G\oo  process  is  short-range  dependent 
and  let  {B(t),  t  >  0}  denote  a  standard  Brownian  motion. 


Theorem  3.3.1  (Short-range  dependence)  If  E  [a2]  <  oo,  then  with  (r  = 
yfr,  r  =  1,2,...,  it  holds  that 


{ sr(t ),  t  >  0} 


'cE 


E  cr 


B(t),  t  >  0}. 
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The  remaining  results  are  obtained  under  the  additional  assumption  that  the 
tail  of  a  is  regularly  varying  of  order  a  (1  <  a  <  2),  i.e.,  of  the  form 

P  [a  >  n]  =  n~ah(n),  n=  1,2,...  (3.9) 


for  some  slowly  varying  function  h  :  H+  — »  H+,  in  which  case  the  mean  E  [a]  is 
finite.  Of  particular  interest  for  the  forthcoming  discussion  is  the  realization  that 
the  truncated  second  moment  of  a  is  (2  —  a)-regula.rly  varying.  Writing 

(  d 

I  -  h(x)  if  f  <  a  <  2 

2  —  a;  v  ' 


x  :=  < 


h(r) 


r= 1 


for  all  x  >  0,  we  can  show  via  Proposition  3 
is  slowly  varying  and  that  whenever  E  [cr2]  = 


(3.10) 

if  a  =  2 

>.l  that  the  function  la  :  R+  — »  1R+ 
oo,  it  holds 


E  1  [a  <  n]a2  ~  n2  ala(n )  (n  — >  oo).  (3-11) 


The  details  of  the  proof  of  this  asymptotic  equivalence  are  identical  to  those  of 
(3.43)  and  (3.44). 

The  next  proposition  handles  the  boundary  value  a  =  2,  which  represents  a 
hybrid  case  between  short-  and  long-range  dependence. 


Theorem  3.3.2  Assume  a  =  2  in  (3.9)  with  E  [a2]  =  oo.  Then,  with  {Cr,  r  = 
1,2,...}  satisfying 


lim  ^Z2(Cr)  =  Jim  y^E  [l  [cr  <  Cr]  o'2]  =  K 


c 


for  some  positive  constant  K,  it  holds  that 


{Sr(t),  t  >  0}  =^r  {) 


cK 

EM 


B(t),  t  >  0}. 


(3.12) 
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Finally,  we  turn  to  the  case  of  bona  fide  long-range  dependence,  i.e.,  1  <  a  <  2. 
We  recall  Definition  1.3.7  and  let  {La(t),  t  >  0}  denote  a  spectrally  positive,  ca¬ 
stable  Levy  motion,  i.e.,  an  a  stable  Levy  motion  such  that  for  all  t  >  0,  the  rv 
La(t )  is  a  stable  rv  Sa(t1/ai  1,0)  [53,  p.  9]  characterized  by 

E  [exp(i9La(t))}  =  exp  t\9\a  ^1  —  i  sgn(d)  tan(— ,  $GR.  (3.13) 


Theorem  3.3.3  (Long-range  dependence)  If  1  <  a  <  2  in  (3.9),  then  with 
{Cr ,r  =  1,2, . . .}  satisfying 


lim  — h((r)  =  lim  rP  [a  >  ([.r]  =  K 

r— >oo  Ca  r— >oo 

Sr 


(3.14) 


for  some  positive  constant  K ,  it  holds  that 


2  —  a. 


{Sr(t),  t  >  0}  =>r  {  Ck  cos(7t — - — )  La(t),  t>  0  } 


1  /a 


where 


Ck  — 


cAT(2  -  a) 


(3.15) 


(3.16) 


,K'  (a-l)BH' 

We  close  with  a  characterization  of  the  scaling  sequences  encountered  in  Theorems 
3.3.2  and  3.3.3;  its  proof  is  given  in  Proposition  A. 3  of  Appendix  A. 


Proposition  3.3.1  The  scaling  sequence  {Cr,  r-  =  1,2,...}  of  Theorems  3.3.2  and 
3.3.3  is  1  / a-regularly  varying,  1  <  a  <  2,  i.e.,  of  the  form  (r  =  r1/“h(r)  for  some 
slowly  varying  function  h  :  R+  — »  R+. 


3.4  Consequences  and  comments 

Several  interesting  inferences  follow  from  the  heavy  traffic  results  obtained  so  far. 
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3.4.1  Queue  size 

We  start  with  the  heavy  traffic  behavior  of  the  normalized  queue  length.  Whenever 
the  convergence  (3.8)  holds,  we  can  immediately  conclude  from  (3.7)  and  from  the 
continuity  of  the  reflection  mapping  (via  the  continuous  mapping  theorem  [5,  Thm. 
5.1,  p.  30])  that 

{Qr(t),  t  >  0}  =^r  {Q(t),  t  >  0}  (3.17) 

with 

Q(t)  :=  S(t)  —  yt  —  inf  (S(x)—'yx),  t  >  0.  (3.18) 

0<x<t 

The  form  of  the  limit  derives  from  (3.7)  and  the  fact  that  lini7r(f)  =  —yt  under 

r— >-o o 

Assumption  (B). 

This  observation  can  now  be  used  to  provide  a  characterization  of  Q( oo),  the 
steady-state  buffer  content  in  heavy  traffic,  under  the  assumptions  of  Theorems 

3.3.1  3.3.3. 

In  the  short-range  dependent  case,  Theorem  3.3.1  combines  with  a  classical 
result  on  the  reflection  functional  of  Brownian  motion  [27,  p.  15]  to  yield  the 
following. 

Theorem  3.4.1  Under  the  assumptions  of  Theorem  3.3.1,  the  resulting  stationary 
heavy-traffic  buffer  content  is  exponentially  distributed,  with 

P  [<3(oo)  >  x\  —  exp  x^j  >  ^  >  0. 

Theorem  3.3.2  leads  via  (3.17)-(3.18)  to  a  similar  result. 

Theorem  3.4.2  Under  the  assumptions  of  Theorem  3.3.2,  the  resulting  stationary 
heavy-traffic  buffer  content  is  exponentially  distributed,  with 

P  [Q(oo)  >  x]  —  exp  ^ ^  ^  x^j  ,  x  >  0. 
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Finally,  in  the  stable  case,  we  need  to  introduce  the  class  of  Mittag-Leffler 
functions  [17,  p.  206]:  For  each  v  >  0  the  Mittag-Leffler  function  Ev  :  IR,  — >  1R,  is 
given  by 

OO 

K{x)  :=X)r(^  +  1)»  :rGR-  (3-19) 

Theorem  3.3.3  can  be  combined  with  established  facts  on  the  reflection  functional 
of  a  Levy  process  [6]  to  yield  the  following  conclusions. 

Theorem  3.4.3  Under  the  assumptions  of  Theorem  3.3.3,  the  distribution  of  the 
resulting  stationary  heavy-traffic  buffer  content  is  given  by 

P  [<2(oo)  >  x}  =  Ea_ i  a;a~1^  i  x  >  0,  (3.20) 

and  the  associated  heavy-traffic  buffer  asymptotics  are  hyperbolic,  with 

cK 

P  [Q(oo)  >  x]  ~  x 1~a  (x  — >  cx)).  (3.21) 


Proof.  Combining  Proposition  5a  of  [6,  p.  725]  (or  Theorem  in  [26,  p.  417])  with 
(3.18)  and  Theorem  3.3.3,  we  obtain 


E  [e-sQ(oo)'  = 


7 


— i -  ,  s  >  0. 

7—1  7  — 


7  +  CKsa 

Application  of  Fubini’s  theorem  in  (3.22)  yields 

P  [Q(oo)  >x\  dx  =  -  (1  —  E  [e-^(°°)]) 

Cr 


s 

1 


s  CK  +  ys 


and  a  simple  change  of  variable  gives 


f°°  T 

■  .  .  X- 

/  e-1  P 

Q( oo)  >  - 

/  0 

L  s  J 

dx  = 


C 


K 


Ck  +  7s 


1 —a 


s  >  0. 


(3.22) 


—  ,  a  >  0  (3.23) 


(3.24) 
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Letting  s  — *  oo  we  note  that  the  mapping  x  — >  P  Q( oo)  > 
increases  to  the  constant  mapping  x  — *  P  [Q(oo)  >  0],  whence 

/» rv* i 

x 


X 


monotonically 


P  [Q(oo)  >  0]  =  /  e  3:  lim  P  Q(oo)  >  — 

Jo  s^°°  L  s 

poo 


=  lim  /  e  x  P 

s — yoo  /n 


Q( oo)  > 


a; 


dx 


dx 


C, 


K 


—  lim  ■  . 

s— >oo  (_yK  'ys1  a 

=  1 

by  the  monotone  convergence  theorem,  so  that  Q(oo)  has  no  point  mass  at  0. 
For  |s|  >  (Ck/ the  right-hand  side  in  (3.23)  can  be  represented  by  an 
absolutely  convergent  geometric  series,  so  that 


poo  -I  00 

/  e~sx  P  [Q(oo)  >  x\  dx  =  - 

S  n=0 


7 


C 


K 


s(1"a)n,  |s|  >  (CK/ 7)1/(a-1}. 


Therefore,  by  Theorem  35.2  in  [13,  p.  192],  the  transform  can  be  inverted  term  by 
term  to  yield 


P  [Q(oo)  >  x]  = 


7 


X 


(a—l)n 


x  >  0 


n=0  v  CK)  r((a-l)n+l)  ’ 

and  (3.20)  readily  follows  from  the  definition  (3.19).  The  asymptotics  (3.21)  are 


verified  by  observing  that 

1  -  E  [e"'Q(oo)] 


s«-i 

7 


(■ 8  0+) 


and  by  making  use  of  a  standard  Tauberian  result  [7,  Corollary  8.1.7]. 


3.4.2  On  selecting  the  heavy  traffic  scaling 

As  the  appropriate  scaling  sequence  {([r,  r  =  1,  2, . . .}  is  revealing  of  the  nature  of 
the  limiting  heavy  traffic  process  {S(t):  t  >  0},  we  brieffy  discuss  here  its  selection. 
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In  Section  2.1.2  we  mentioned  that,  under  (3.9)  with  1  <  a  <  2,  the  M\G\oo 
busy  server  process  possesses  the  second  order  asymptotic  self-similarity  property, 
with  parameter  (3— oe)/2,  i.e.,  by  aggregating  the  original  process  {bn,  n  —  0, 1, . .  .} 
in  blocks  of  size  m  and  dividing  by  the  block  size,  we  obtain  in  the  limit  (as  m 
goes  to  infinity)  the  same  correlation  function  as  that  of  a  fractional  Gaussian 
noise  process.  Such  convergence  of  the  correlation  function  tempts  one  to  think 
that  the  appropriate  scaling  ensuring  (3.8)  might  be  the  one  that  balances  the 
rate  of  growth  of  the  partial  sums  variance,  so  that  convergence  (3.8)  occurs  to  a 
limiting  process  with  finite  variance.  We  now  explore  this  point  in  some  detail: 

By  standard  calculations,  we  find  the  variance  of  the  partial  sums  to  be 


var  [5r(t)] 


ArE  [a 

C2 

Sr 


/  M  \ 

( + 2  ~  ^)p  p  > ) » 


t  >  0 


for  all  r  =  1,2, . . ..  It  can  be  shown  that  when  the  tail  of  a  satisfies  (3.9)  with 
1  <  a  <  2,  the  candidate  scaling  {()r,  r  =  1,  2, . . .}  given  by 

r 

C  ■■=  r^p  [a  >  k]  ,  r  =  1,2,...  (3.25) 

k= 1 

indeed  results  in  a  finite  limiting  variance,  i.e.,  lim  var  [<S'r(f)]  exists  and  is  hnite 

r— >oo 

for  all  t  >  0.  In  addition,  invoking  (3.94)  we  see  that  the  scaling  (3.25)  has  the 
asymptotic  form 


(2  —  a)(a  —  1)E  [a 


[r  — »  oo 


(3.26) 


and  is  therefore  regularly  varying  of  order  (3  —  a)/ 2. 

On  the  other  hand,  from  Theorem  1.3.1  we  already  know  that  convergence  of 
a  normalized  partial  sum  process,  such  as  {Sr(t),  t  >  0},  can  only  be  to  a  self¬ 
similar  process,  and  that  the  corresponding  Hurst  parameter  H  may  be  determined 
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through  the  regularly  varying  scaling  { .  r  =  1,  2, . . .}  by 

i-  S[rx]  H  >  n 

Inn  -J—L  =  x  ,  x  >  0. 

r— >-oo  Qr 

Thus,  the  candidate  scaling  (3.25),  which  balances  the  growth  of  the  variance,  sug¬ 
gests  possible  convergence  to  a  fractional  Brownian  motion  with  Hurst  parameter 
(3  —  a)/2. 

In  the  present  heavy  traffic  setup  however,  convergence  of  the  rescaled  M\G\oo 
process  to  a  fractional  Brownian  motion  does  not  take  place.  The  candidate  scaling 
(3.25)  is  not  the  appropriate  scaling;  it  is  too  strong  and  yields  convergence  to 
a  degenerate  limit  -  the  identically  zero  process.  Theorem  3.3.3,  in  conjuction 
with  Proposition  3.3.1,  clearly  shows  that  the  correct  scaling  does  not  contain 
any  r^3-0^2  factor,  but  instead  contains  the  weaker  rl^a  factor  associated  with 
the  stable  law  to  which  the  service  rv  a  is  attracted.  As  a  result,  the  limiting 
heavy  traffic  process  turns  out  to  be  not  a  fractional  Brownian  motion  but  an 
a  stable  1/a-  self  similar  Levy  motion,  the  stable  analog  of  standard  Brownian 
motion,  which  has  independent  increments  with  infinite  variance.  In  heavy  traffic, 
the  corresponding  queue  length  asymptotics  are  not  Weibullian,  but  hyperbolic 
with  power  1  —  a.  Thus,  M\G\oo  processes  demonstrate  that,  within  long-range 
dependence,  fractional  Brownian  motion  does  not  assume  the  ubiquitous  role  that 
standard  Brownian  motion  plays  in  the  short-range  dependence  setup,  and  that 
modeling  possibilities  attracted  to  non-Gaussian  limits  are  not  so  hard  to  find. 
Clearly,  the  extent  to  which  such  non-Gaussian  processes  can  serve  as  useful  traffic 
models  deserves  further  consideration. 
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3.5  Outline  of  proof  and  preliminary  results 


In  this  section  we  organize  the  proof  of  Theorems  3.3. 1-3. 3. 3  into  a  series  of  steps 
which  we  formalize  as  Propositions;  their  proofs  are  given  in  Section  3.7. 

Look  at  the  rth  queueing  system  for  some  r  =  1,2,...,  and  fix  n  =  0, 1, . . .. 
We  note  the  decomposition  b rn  =  bfflr  +  bbf,r  where  the  rvs  b^r  and  b^)r  describe 
the  contributions  to  the  number  of  customers  in  the  system  at  the  beginning  of 
slot  [n,  n  +  1)  from  those  initially  present  (at  n  =  0)  and  from  the  new  arrivals, 
respectively.  It  is  easy  to  see  that 

br  n  PI 

bny  =  J21  >  aIld  ^  =  X  X  1  >n~k}- 

j= 1  k= 1  j=  1 

It  was  shown  in  [47,  Sec.  5]  that 

n  br 

*n)r  '■=  bT]  =  X  min(n’  _  X)  (3-27) 

j= 1  3= 1 

and 

n  n  Pk 

sn)r  '■=  X  bk)r  =  X  X  min(Wj,  n-k+1).  (3.28) 

k=l  k= 1  j= 1 

We  introduce  the  rescaled  versions 


and 


so  that 


S{0)r{t) 


1 

c 


t  >  o 


S^a)r(t) 


1 


t  >  0 


ST(t)  =  5,(0)r(t)+'5(o)r(t),  f  >  0. 


(3.29) 


Also,  for  each  T  >  0,  the  identically  zero  mapping  on  [0,  T]  is  the  element  of  D[ 0,  T] 
denoted  by  dT,  i.e.,  dT  :  [0,  T]  — >  R  with  9T(t)  =  0,  0  <  t  <  T. 
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We  first  show  that  the  initial  condition  plays  no  role  in  the  heavy  traffic  limit, 
as  should  be  expected.  This  reduction  step,  as  well  as  others  taken  in  this  section, 
is  accomplished  under  the  following  sufficient  condition. 


Condition  (B)  The  scaling  sequence  {Cr?  r-  =  1,2,...}  satisfies 


lim  7EP[?>^  =  0' 

r— > oo  (5,r  L ' 


3  = 1 


Condition  (B)  holds  under  each  set  of  assumptions  of  Theorems  3.3. 1-3. 3. 3;  this 
is  shown  in  Proposition  3.6.1  of  Section  3.6. 


Proposition  3.5.1  Under  Condition  (B),  for  each  T  >  0  it  holds  that 

{S,(0)r (t) ,  0  <  t  <  T}  r  dT  in  D[0,  T\. 

Thus,  in  order  to  get  (3.8)  it  suffices  to  consider  the  limiting  behavior  of  the 
rescaled  process  {S'^r(t),  t  >  0}.  To  that  end,  for  each  r  =  1,2,.. .,  we  introduce 
the  sequence  {wrn,  n  —  0, 1, . . .}  given  by 

n  hi 

w£:=0,  wn'-=^2^2°k,j,  n  =  1,2,...  (3.30) 

fc=i  j= i 

which  can  be  interpreted  as  the  sequence  of  partial  sums  associated  with  the  in¬ 
stantaneous  arrivals  (2.16).  The  corresponding  rescaled  process  {Wr(t),  t  >  0}  is 
given  by 

W'(t)  ■-  f  (w\H]  -  E  [mfrt|] ) .  t  >  0.  (3.31) 

Sr 

The  main  idea  driving  the  discussion  is  that  in  as  much  as  heavy  traffic  is  concerned, 
the  process  {Wr(t),  t  >  0}  acts  as  a  surrogate  for  {S^^t),  t  >  0}.  This  is  made 
precise  through  the  following  “convergence  together”  result. 
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Proposition  3.5.2  Under  Condition  (B),  for  each  T  >  0  it  holds  that 

{' Wr{t )  -  S{a)r(t),  0  <  t  <  T}  r  eT  in  D[ 0,  T\. 

Combining  Propositions  3.5.1  and  3.5.2,  we  immediately  get  the  following  con¬ 
clusion  from  the  decomposition  (3.29). 

Corollary  3.5.1  Under  Condition  (B),  for  each  T  >  0  it  holds  that 

{ Wr(t )  -  Sr(t),  0  <  t  <  T}  =¥> r  0T  in  D[ 0,  T], 

so  that  the  process  {Sr(t),  0  <  t  <  T}  is  S-weakly  convergent  if  and  only  if 
{Wr(t),  0  <  t  <  T}  is  S -weakly  convergent,  and  convergence  is  to  the  same  limit. 

Thus,  we  need  only  consider  the  convergence  of  the  process  {Wr(t),  t  >  0},  and 
characterize  the  limiting  process.  In  fact,  a  further  reduction  can  be  achieved  by 
noting  that  in  heavy  traffic  we  can  replace  {(3rk,  k  =  1,  2, . . .}  by  the  limiting  i.i.d. 
sequence  {/?*,,  k  =  1.2,...},  where  the  generic  rv  f3  is  a  Poisson  rv  with  parameter 
c/E  [cr] .  More  precisely,  consider  the  modified  workload  process  { vn ,  n  —  0, 1, . . .} 
given  by 

n  Pk 

Vo  =  0;  vn  =  V,  V,  crkj,  n  =  1,2,...  (3.32) 

k= 1  j=l 

under  the  assumption  that  the  rvs  {/?*,,  k  =  1,2,...}  are  independent  of  the  session 
duration  rvs  {<rnj-,  n, j  =  1,2,...}.  For  each  r  =  1,2,...,  the  corresponding 
rescaled  process  {Vr(t),  t  >  0}  is  defined  by 

Vr(t)  :=  y  (v[rt]  -  E  [v[rt]])  ,  t>  0.  (3.33) 

br 

Proposition  3.5.3  Under  Assumption  (B),  the  process  { Wr(t ),  0  <  t  <  T}  is 
S -weakly  convergent  if  and  only  if  {Vr(t),  0  <  t  <  T}  is  S -weakly  convergent, 
and  convergence  is  to  the  same  limit. 
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Corollary  3.5.1  and  Proposition  3.5.3  together  lead  to  the  following  conclusion: 


Corollary  3.5.2  Under  Assumption  (B)  and  Condition  (B),  the  process  {Sr(t), 
0  <  t  <  T}  is  S -weakly  convergent  if  and  only  if  {Vr(t),  0  <  t  <  T}  is  S -weakly 
convergent,  and  convergence  is  to  the  same  limit. 

3.6  Proofs  of  Theorems  3.3.1— 3.3.3 

First,  the  big  picture:  Corollary  3.5.2  and  Proposition  3.6.1  (given  below)  imply 
that  in  proving  Theorems  3.3. 1-3. 3. 3  we  need  only  investigate  the  convergence 
of  the  modified  workload  process  (3.33).  This  is  a  much  easier  task  as  we  now 
deal  with  the  (normalized)  partial  sums  process  associated  with  a  single  sequence 
of  i.i.d.  rvs,  of  finite  mean  but  possibly  infinite  variance,  an  extensively  studied 
situation  where  the  (functional  form  of  the)  classical  Central  Limit  Theorem  and 
its  generalization  to  i.i.d.  summands  with  infinite  variance,  are  expected  to  yield 
the  requested  convergence.  In  fact,  as  we  shall  see  shortly,  the  convergence  of 
the  finite  dimensional  distributions  of  {Vr(t),  t  >  0}  turns  out  to  be  an  easy 
by-product  of  classical  results  concerning  stable  distributions  and  their  domains 
of  attraction  [20,  pp.  574-581].  Finally,  the  desired  S  weak  convergence  of  the 
process  {Vr(t),  t  >  0},  thus  of  {Sr(t),  t  >  0},  will  be  validated  through  functional 
convergence  results  due  to  Skorokhod  [57].  This  approach  clearly  explains  the  form 
of  the  results  obtained  in  this  chapter,  providing  insights  as  to  when  the  process 
{Vr(t),  t  >  0}  is  expected  to  converge,  and  to  which  limit.  A  different,  analytic 
approach  using  characteristic  functions  was  pursued  in  the  technical  report  [60]. 

We  now  proceed  with  the  details:  In  Section  3.8  we  give  a  proof  that  the 
technical  Condition  (B)  required  to  establish  the  “convergence  together”  argument, 
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indeed  holds  under  the  assumptions  of  Theorems  3.3. 1-3. 3. 3. 

Proposition  3.6.1  Condition  (B)  holds  true  for  each  of  the  scaling  sequences 
{Cr,  t  =  1,  2, . . .}  in  Theorems  3. 3. 1-3.3. 3. 


Next,  we  consider  the  generic  compound  rv  Y  given  by 

P 

Y'=Y,r  I3'34) 

3= 1 

where  the  rv  (3  is  a  Poisson  rv  with  parameter  c/E  a\  and  independent  of  the  i.i.d. 
rvs  {(Jj,  .7  =  1,2,...}  which  are  distributed  according  to  a.  Fixing  t  >  0,  we  note 


that 


M 

^ (t)  =st  -  Y {Y*  -  E  ™  >  r  =  1, 2, . . .  (3.35) 

k= i 


where  the  i.i.d.  rvs  {Y*.,  k  =  1,2,...}  are  distributed  according  to  the  generic  rv 


Y. 


For  easy  reference,  we  restate  some  useful  facts  concerning  stable  distributions 
and  their  domains  of  attraction;  the  reader  is  refered  to  [20,  pp.  574-581]  for 
additional  material:  Let  L  be  a  rv  with  distribution  not  concentrated  at  one  point, 
and  let  {Xr,  r  =  1,  2, . . .}  be  a  sequence  of  i.i.d.  rvs,  with  generic  rv  X.  We  say 
that  X  belongs  to  the  domain  of  attraction  of  the  rv  L  if  there  exist  normalizing 
constants  (r  >  0  and  cr,r  =  1,2,...,  such  that 

X  x  -(-•••  T  Xr  —  r  cr 


Cr 


L. 


(3.36) 


By  Theorem  1  of  [20,  p.  576]  only  stable  rvs  possess  a  domain  of  attraction.  By 
Theorem  2  in  [20,  p.  577],  in  order  for  X  to  belong  to  the  domain  of  attraction  of 
a  stable  law  with  exponent  a,  0  <  a  <  2,  it  is  necessary  that  its  truncated  second 
moment  be  regularly  varying  with  exponent  2  —  a,  i.e., 

E[l[X  <  r]  X2]  ~  r2~ag(r)  (r  ->  oo),  (3.37) 
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for  some  slowly  varying  function  g  :  1R+  — )•  H+.  The  associated  scaling  sequence 
{£r,  r  =  1,  2, . . .}  in  (3.36)  must  then  satisfy 

lim  [1  [X  <  Crl  X2]  =  M  (3.38) 

r^ocCZ  L  J 

Sr 

for  some  constant  M  >  0  [20,  p.  579].  Moreover,  if  E  [X]  is  finite,  then  by  Theorem 
3 (ii)  of  [20,  p.  581]  we  can  take  cr  =  E  [X],  r  =  1,  2, . . .. 

We  are  now  ready  to  discuss  Theorems  3.3. 1-3. 3. 3  which  are  all  proven  in 
the  same  manner,  although  for  clarity  of  presentation,  we  shall  consider  each  of 
them  separately.  As  E  [Y]  is  finite  under  the  enforced  assumptions,  we  conclude 
from  (3.35)  and  (3.36)  that  for  each  t  >  0,  the  convergence  question  concerning 
{Vr(t),  r  =  1,  2, . . .}  is  equivalent  to  determining  whether  the  rv  Y  is  attracted 
to  a  stable  law,  and  to  which  one.  In  asserting  this  equivalence  we  rely  on  the 
fact  that  the  scaling  sequence  {Cr,  r  =  1,  2, . . .}  so  selected  is  regularly  varying,  as 
turns  out  to  be  the  case  by  Proposition  3.3.1,  so  that 

lim  =  t1/a,  t  >  0.  (3.39) 

r— »■  oo  Qr 

In  each  case,  we  show  that  both  the  necessary  condition  (3.37)  and  the  accompa¬ 
nying  sufficient  condition  stated  in  [20,  p.  577]  are  satisfied.  This  occurs  simply 
because  the  generic  rv  Y  inherits  the  tail  behavior  of  the  generic  service  time  a 
under  each  set  of  assumptions  of  Theorems  3. 3. 1-3. 3. 3. 

A  proof  of  Theorem  3.3.1. 

Since  E  [cr2]  <  oo,  the  variance  of  Y  is  also  finite,  and  is  given  by 

var  [Y]  =  var  \/3\  E  [a]2  +  E  \/3\  var  [cr]  =  L  (3.40) 

E  p] 

Hence,  the  truncated  second  moment  of  Y  varies  slowly,  i.e.,  (3.37)  holds  with 
a  =  2  and  as  Y  is  never  degenerate  at  one  point,  it  follows  from  Corollary  1  to 
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Theorem  2  in  [20,  p.  578]  that  Y  is  attracted  to  the  normal  distribution.  Obviously, 
the  scaling  (r  =  y/r,  r  =  1,2,...  satisfies  (3.38),  with  M  =  cE  [a2]  /E[a].  In 
fact,  by  a  well-known  result  of  Donsker  [5,  Thm.  16.1,  p.  137],  selecting  = 
y/r,  r  =  1,2, . . .  ensures  that  the  process  {Vr(t),  t  >  0}  is  S'-weakly  convergent 
to  a  Brownian  motion,  with 

{Ur(f),  t  >  0}  =^r  {VM  B(t),  t  >  0}.  (3.41) 


Combining  (3.41)  with  Proposition  3.6.1  and  Corollary  3.5.2  immediately  concludes 
the  proof.  ■ 


Under  the  assumptions  of  Theorems  3.3.2  and  3.3.3,  E  [a2]  is  infinite,  and  the 
compound  Poisson  rv  Y  now  has  infinite  variance.  Also,  if  a  satisfies  the  tail 
condition  (3.9),  so  does  Y  with 


P  [Y  >r 


P 


-  p 


~  E  \/3\  r~ah{r) 


( r  — *  cx)). 


(3.42) 


Li=i  J 

The  asymptotic  equality  in  (3.42)  is  stated  as  an  exercise  in  [20,  Ex.  31,  p.  288], 
where  the  reader  will  find  hints  for  its  proof  (see  also  [16]).  Next,  we  check  that 
the  truncated  second  moment  of  Y  is  given  by 


r  r—1 

E  [l  [Y  <  r }  y2]  =  2  nP  [Y  >  n\  —  r(r  +  2)P  [Y  >  r]  +  £p[^>n] 

n= 1  n— 0 


for  each  r  =  1,  2, . . ..  Using  (3.42)  in  this  last  expression,  we  find  that 

E  [1  [Y  <  r }  y2]  ~  E  \/3\  ^2  ^  n1"“h(n)  -  r2~ah(r)  j  (r  ->  oo)  (3.43) 

because  E  [y]  is  hnite  and  E  (y2]  infinite.  We  close  these  preliminary  remarks  by 
noting  that  the  truncated  second  moments  of  a  and  Y  are  obviously  related  to 
each  other  by 


E  1  [y  <  r]  y2  ~  E  \/3\  E  1  [a  <  r]  a2  (r  — »  oo).  (3.44) 
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A  proof  of  Theorem  3.3.2. 

Inserting  a  =  2  in  (3.43)  and  using  the  definition  (3.10)  (with  a  —  2),  we  get 
E  [1  [Y  <  r]  T2]  ~  E  \/3\  ( l2(r )  -  h(r))  (r  ->  oo). 


By  Proposition  3.8.1  (ii) ,  Z2  :  K+  — >  H+  is  slowly  varying  with 


lim 

r— >•  00 


h(r) 

Z2(r) 


0, 


so  that 

E  [1  [Y  <  r]  T2]  ~  E  \J3\  l2(r )  (r  ->  00). 


(3.45) 


This  time,  by  Corollary  1  in  [20,  p.  578],  the  slow  variation  of  the  truncated  second 
moment  is  a  necessary  and  sufficient  condition  for  Y  to  be  attracted  to  the  nor¬ 
mal  distribution,  with  normalizing  coefficients  selected  according  to  (3.38)  (despite 
the  fact  that  the  variance  of  Y  is  now  infinite).  Since  the  marginals  of  the  process 
{Vr(t),  t  >  0},  which  has  stationary,  independent  increments,  converge  to  a  Gaus¬ 
sian  distribution,  it  follows  by  [57,  Theorem  2.7]  without  any  additional  conditions 
that  (3.41)  takes  place.  Because  of  (3.45),  selecting  the  scaling  {£r,  r  =  1,2, . . .} 
according  to  (3.38),  with  M  =  cK/Y  [<x],  is  equivalent  to  (3.12).  Combining  (3.41) 
with  Proposition  3.6.1  and  Corollary  3.5.2  completes  the  proof.  ■ 


A  proof  of  Theorem  3.3.3. 

When  1  <  a  <  2  in  (3.9)  the  rvs  a  and  Y  have  infinite  second  moment,  and 
Proposition  3.8. l(i)  implies 


lim 

r—¥  00 


1 

r2~ah{r ) 


n 1  ah{n) 

n= 1 


1 

(2 -a)' 


(3.46) 
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Using  this  asymptotic  in  (3.43)  we  get 

Q, 

E  [1  [Y  <  r }  U2]  ~  - - E  \J3\  r2~ah(r)  (r  ->  oo).  (3.47) 

Zj  cy 

Invoking  Corollary  2  of  [20,  p.  578],  we  see  that  (3.47)  and  the  tail  condition 
(3.42)  are  sufficient  to  ensure  membership  of  Y  in  the  domain  of  attraction  of  a 
non-normal  stable  distribution  with  exponent  1  <  a  <  2.  The  associated  scaling 
sequence  {([r,  r  —  1,2,.. .},  selected  according  to  (3.38),  yields  convergence  of  the 
marginal  distribution  of  V ( 1 ) ,  as  r  goes  to  infinity,  to  that  of  an  a-stable  rv,  i.e., 

limE  [exp(i6»Ur(l))]  =  E  exp  (id  **)  cos^2  La(l)  ] 

r^oo  \  \  a{a  —  1)  2  /  / 

for  all  9  in  H.  The  exact  value  of  the  constant  given  above  can  be  easily  veri¬ 
fied,  by  recalling  the  expression  (3.13)  for  the  characteristic  function  of  La{  1)  and 
comparing  it  with  Eq.  (3.18)  of  [20,  p.  730]  (note  the  unfortunate  error  in  the  ± 
sign).  Next,  appealing  to  [57,  Theorem  2.7]  again,  we  conclude  that  convergence 
of  the  marginals  also  implies  S'- weak  convergence  of  the  process  |l/r(f),  t  >  ()}, 
which  has  stationary,  independent  increments,  to  an  a-  stable  Levy  motion.  More 
precisely,  it  holds  that 

t>  0}  {  ( M((3  ~  cos(x~ 1A>  L„(t),  t  >  0}.  (3.48) 

Using  (3.47)  in  (3.38)  with  M  =  cKa/( 2  —  cr)E  [a]  we  obtain  the  scaling  sequence 
{Cr,  r  =  1,2, . . .}  given  in  (3.14).  Finally,  combining  (3.48)  with  Proposition  3.6.1 
and  Corollary  3.5.2  shows  that  (3.15)  holds  true.  ■ 
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3.7  Proofs  of  Propositions  3.5.1,  3.5.2  and  3.5.3 


We  start  by  remarking  that  if  the  sequence  { (r .  r  =  1,  2, . . .}  is  regularly  varying 
(as  stated  in  Proposition  3.3.1),  then  Condition  (B)  also  implies 

1  [rf] 

lim  —  P  [<t  >  j]  =  0,  t  >  0.  (3.49) 

r— >oo  Cr  ^ ^ 

Sr  3= 1 

All  three  proofs  given  in  this  section  follow  the  same  pattern,  and  are  based  on 
the  following  simple  idea:  Consider  a  sequence  of  1R- valued  processes  {Xr(t),  t  > 
0},  r  =  1,2,...,  with  sample  paths  in  D[ 0,  oo).  Fix  T  >  0.  According  to  Theorem 
4.1  of  [5,  p.  25],  the  17-weak  convergence 

{Xr(f),  0  <  t  <  T}  r  dT  in  D[ 0, T], 
follows  from  the  convergence  in  probability 

sup  \Xr(t)\  >r  0.  (3.50) 

o  <t<T 


A  proof  of  Proposition  3.5.1. 

Fix  r  =  1,  2, . . .,  and  note  from  (3.27)  that 

br 

sup  |5'^r(i)|  <  'S''  min (dj  —  1,  [ rT ])  +  —  E  [min  (A  —  1,  [ rT ])]  . 

0 <t<T  C r  (r 

Hence,  for  every  e  >  0,  it  is  plain  that 


P 


sup  \s{0)r(t)\>  £ 

.0  <t<T 


< 


< 


P 


y  jrminfo-  -l,[rT])  + 
l>  j=l 


ArE  [a 

Cr 


E  [min  ( a 


2  Ar  E  <7 

E-Cr 


E  [min  (A 


1.  [rT])] 


1,  [rT])]  >  £ 

(3.51) 
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where  the  last  step  follows  by  Chebyshev’s  inequality.  It  is  also  the  case  that 

[rT]- 1 

E  [min(<r  —  1,  [rT])}  =  ^  P  [min(<j  —  1,  [rT])  >  n] 

n= 0 

[rT\~  1 

=  X]  p  p  ~  1  > 

n= 0 
[rT] 

=  X  p  p  > n]  • 

n=l 

Appealing  to  Condition  (B)  and  (3.49),  we  get 

lim  E  —  min  (a  —  1,  [rT])  =0  (3.52) 

r— >oo 

and  the  conclusion 


lim  P  sup 


\S{0)r{t)\>  £  =0 


immediately  obtains  from  (3.2)  upon  letting  r  go  to  infinity  in  (3.51). 


A  proof  of  Proposition  3.5.2. 

Fix  r  =  1,  2, _  From  (3.28)  and  (3.30)  we  note  that 

n  PI 

k= 1  j= 1 

n  Pn-k+l 

=  X  X  K-fe+u  -  k)+ 

k=  1  J=1 

n  0J 

XXK,-^  n  =  l,2,...  (3.53) 

k= 1  j=l 

where  the  last  step  made  use  of  the  mutual  independence  of  the  families  of  i.i.d. 
rvs  {/3[,  k  =  1,2,...}  and  {<Tk,j,  k,j  =  1,2,...}.  It  is  now  straightforward  to 


check  that 


,  [rT]  PI 

sup 


,  VT] 

^ r  k=  i 
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By  Chebyshev’s  inequality,  for  every  e  >  0  we  obtain 


P 


sup  \Wr{t)  ~S{a)r(t)\  >£ 

.0  <t<T 


< 


2Ar 


[rT] 


E  E  [(»■  -  fc)+] 


2ArE  <T 
E-Cr 


[rT] 


p  ^  ’ 

fc=i 


and  the  desired  convergence 


(3.54) 


sup  \VT(t)  -  S'(a>r(t)|  -A,  0 

0  <t<T 

follows  upon  letting  r  go  to  infinity  in  the  upper  bound  (3.54),  and  making  use  of 
(3.2),  Condition  (B)  and  (3.49).  ■ 


The  proof  of  Proposition  3.5.3  requires  estimates  that  derive  from  various  mar¬ 
tingales  inequalities;  we  now  state  them  in  Lemmas  3.7.1  and  3.7.2  for  easy  ref¬ 
erence:  Consider  a  collection  of  integrable  rvs  {Xj,  i  =  l,...,n}  adapted  with 
respect  to  the  filtration  {jTj,  i  =  1, . . . ,  n},  i.e.,  for  each  i  =  1, . . . ,  n,  the  rv  Xt  is 
measurable.  We  also  write 

Si  =  A  i  +  . . .  +  Xi ,  i  —  1 , . . . ,  n. 

Kolmogorov’s  maximal  inequality  [25,  Corollary  2.1,  p.  14]  is  given  first. 

Lemma  3.7.1  Assume  {(Si,  A*),  i  —  1, . .  .,n}  to  form  a  martingale.  Then,  for 
each  p  >  l,  it  holds  that 

<  A-pE[|Sn|p],  A  >  0. 

The  von  Bahr-Esseen  inequality  [64]  is  next. 


max  \Si 


>  A 
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Lemma  3.7.2  Assume  {(X* ,Ti),  i  —  1, . . .  ,n}  to  form  a  martingale  difference. 
7/E  [|Xj|p]  <  oo  for  all  i  —  1, . . . ,  n,  then 

n 

E  [\Sn\p]  <  2^E  [|Xi|p] ,  l<p<2. 

i= 1 

In  what  follows  Lemmas  3.7.1  and  3.7.2  are  applied  to  the  special  case  when 
the  rvs  {Xj,  %  —  1, . . . ,  n}  are  zero-mean  i.i.d.  rvs. 

A  proof  of  Proposition  3.5.3. 

Recall  that  the  rvs  {/3k,  k  =  1,2,...}  are  i.i.d.  Poisson  rvs  with  parameter  c/E  [cr] , 
which  are  independent  of  the  sequence  of  i.i.d.  session  duration  rvs  k,j  = 

1,2,...}. 

Fix  r  =  1,2,....  On  the  same  probability  triple  (0,X,  P)  where  the  previ¬ 
ously  mentioned  rvs  are  defined,  we  introduce  a  family  of  i.i.d.  {0,  l}-valued  rvs 
{Urk)],  k,j  =  1,2,...},  i.e., 

p  [ur  =  i]  =  Ar^  M  =  i  -  p  [ur  =  o] 

where  Ur  denotes  the  generic  rv  for  this  i.i.d.  sequence.  The  rvs  [Ufj,  k,j  = 

1,2,.. .}  are  assumed  independent  of  the  collections  of  rvs  mentioned  so  far.  Next, 
we  define  the  rvs  {(3rk,  k  =  1,  2, . . .}  by 

dk 

k  =  1,2,... 

3— 1 

We  also  define  the  workload  process  {wrn,  n  —  0, 1, . . .}  corresponding  to  {/ 3rk ,  k  = 

1,2,...}  by 

n  hi 

^o:=0>  W'i  :=  ^2^2  a  ko,  n  =  1,2,... 

fc=i  j= i 
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and  its  rescaled  version  {Wr(t),  t  >  0}  by 

W'(f)  :=  i  (^  -  E  [fi? [rt|]),  t>  0. 

Sr 

Under  the  enforced  independence  assumptions,  it  is  easy  to  check  that  {/3£,  k  = 
1,2,...}  =st  {/?£,  k  =  1,  2, . . .},  and  that 

0  <  t  <  T}  =st  { Wr[t ),  0  <  t  <  T}. 

Moreover,  these  rvs  are  all  defined  on  the  same  probability  triple  as  the  rescaled 
process  (Vr(t),  t  >  0}.  Thus,  the  result  will  be  established  if  it  holds  that 

{Vr{t)  -  Wr(t),  0  <t<T}  =^r  0T  in  D[0,  T], 


or  equivalently,  if  we  can  show  that 


sup  | Vr(t)  -  Wr{t)\  -Ar  0.  (3.55) 

0  <t<T 

To  that  end,  for  each  r  =  1,  2, . . .,  we  note  from  the  definitions  that 

n  /3k 

=  n  =  1,2,... 
k=  1  j= 1 


so  that 

n  0k 

Vn-™h  = 

k=  1  j= 1 

The  rvs  k  —  1, 2, . . .}  defined  by 


Ukj)°kj, 


n  =  1,2, . . . 


(3.56) 


0k 

:=  1  -  *  =  1,2,...  (3.57) 

j=l 

are  i.i.d.,  and  denote  by  Zr  the  corresponding  generic  rv  associated  with  this 
collection  of  rvs.  It  is  plain  from  (3.56)  and  (3.57)  that 


sup  \Vr(t)  -Wr[t)\ 

0  <t<T 


1 

—  sup 

Sr  l<n< [rT] 


-  E  Kl) 

fc=l 
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Fix  e  >  0.  Invoking  the  maximal  inequality  for  martingale  sequences  [Lemma 


3.7.1],  we  get 


T  _  l  i  [rTl 

P  sup  \\T(t)-W'(t)\>e  <——  E  £(Z£-E[Z£]) 


(3.58) 


with  p  selected  such  that  1  <  p  <  a  <  2.  This  selection  of  p  ensures  E  [ap]  < 
oo  both  under  short-range  dependence  and  under  the  assumption  of  regularly 
varying  tail  (3.9).  The  von  Bahr  -  Esseen  inequality  [Lemma  3.7.2]  for  martingale 


differences  can  now  be  applied  to  the  right-hand  side  of  (3.58)  to  yield 


(e(r)p 


L-4  2  IrTl 

E  £(ZJ-  E[Z3)  <JmE[|Z'-E[Z 


(3.59) 


By  the  convexity  of  x  — >  xp  (p  >  1)  on  1R+,  we  find 


E  [\Zr  -  E  [Zr]  |p]  <  2P_1(E  [| Zr\p ]  +  E  [| Zr\]p )  <  2PE  [\Zr 


(3.60) 


with  the  last  step  validated  by  Jensen’s  inequality.  Next,  using  the  definition  of 

Z' ,  we  obtain  by  the  same  convexity  argument  that 

-  0 

E[\Zr\p\(3}<  (3P-1  E  J^(l  -  U])pap  \(3  =  /3PE  [<jp]  E  [(1  —  Ur)p]  a.s.  (3.61) 

-3= 1 

under  the  enforced  independence  assumptions  (and  with  an  obvious  notation). 
Injecting  the  bounds  (3.60)  and  (3.61)  into  (3.59),  we  conclude  from  (3.58)  that 

P  sup  I  V(t)  -  Wr(t)  I  >  s  <  T  — ^E  H  E  m  -E  [(1  -  U’Y\  (3.62) 

_0<t<T  J  £pQr  r  Sr 

As  the  heavy  traffic  Assumption  (B)  implies 

lim^E[(l  -  Ur)p]  =  -,  (3.63) 

r— >oo  Qr  C 

the  desired  conclusion  (3.55)  now  follows  by  letting  r  go  to  infinity  in  (3.62)  and 
noting  that  lim  1/C P_1  =  0  for  p  >  1.  ■ 

r— >•  oo 
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3.8  A  proof  of  Proposition  3.6.1 


In  the  proof  of  Proposition  3.6.1  and  elsewhere,  we  make  use  of  the  following  fact. 
Lemma  3.8.1  For  any  slowly  varying  function  u  :  R+  — *  R+,  it  holds  that 

lim  xpu{x)  =  cx),  p  >  0,  (3.64) 

x—>oo 

while 

lim  xpu(x)  =  0,  p  <  0.  (3.65) 

x—>oo 


Proof.  By  the  Representation  Theorem  for  slowly  varying  functions  [7,  Theorem 
1.3.1,  p.  12],  we  can  write 


rx  £(t ) 

u(x)  ~  cexp  (  /  — -  dt  )  (x  —f  cx)) 

I A  t 


(3.66) 


with  constants  A  >  0  and  c  >  0,  and  Borel  mapping  £  :  R+  — >  R  such  that 


lime(t)  =  0.  Thus, 

t — yoo 


x  n 

xpu(x)  ~  cAp  exp  (  j  - — - — —  dt  )  (x  — >  oo). 


For  every  8  >  0  there  exists  t$  >  A  such  that  \e(t)\  <  8  for  t  >  t$,  whence 


±±E  <  £W±r  <  i±£  t>ts 

t  -  t  -  t  ’ 


so  that 


K  +  (-<5  +  p)  In  (  ^  +  -  dt  <  K  +  (8  +  p)  In  ( )  ,  x  >  ts 

JS  J  J  A  *  \*( 5. 


with 


if  :=  [“  £M±p  it. 
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The  conclusion  (3.64)  (resp.  (3.65))  follows  from  these  inequalities  when  selecting 
8  >  0  such  that  8  <  p  (resp.  8  <  — p )  -  such  a  selection  is  always  possible  when 
p  >  0  (resp.  p  <  0).  ■ 


The  limit  (3.64)  is  useful  in  the  proof  of  the  following  discrete  analogue  to  the 
direct  half  of  Karamata’s  Theorem  [7,  p.  26]. 


Proposition  3.8.1  Let  u  :  R+  — *  R+  be  a  slowly  varying  function.  Then  the 
following  statements  hold: 

(i)  For  any  p  >  —1,  we  have  the  asymptotics 

r  rp+ 1 


E 

n= 1 


muyn) 


p+ i 


u(r )  (r  — »  cx)); 


(ii)  For  any  p  <  —1,  we  have  the  asymptotics 

rP+ 1 


^^npu(n)  ~ 


p  +  1 


u(r )  (r  — »  cx)); 


(in)  The  mapping  u  :  R+  — »  R+  defined  by 


u{x):=  x  >  0 

n=  1 


is  a  slowly  varying  function  which  satisfies 


,  u  (x) 

lim  —7  r  =  0. 


x^-oo  u(x) 


(3.67) 


(3.68) 


(3.69) 


(3.70) 


Proof,  (i)  Under  the  condition  p+1  >  0,  it  follows  from  (3.64)  that  (3.67)  is 
equivalent  to 


lim  y 

i — s-oo  ^  r  \  r  I 


r-+oo  * — '  r\rJ  u(r)  p+1 

n=r*+ 1  v  ;  ^ 


(3,71) 


for  any  finite  r*  (which  is  now  fixed  for  the  remainder  of  the  proof). 
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To  prove  (3.71)  pick  5  in  (0, 1  +  p)  and  A  >  0.  By  Potter’s  bound  [7,  p.  25], 
there  exists  an  integer  r*  :=  r(S,  A)  such  that 


u(n)  .  fn\~& 

.  !  <  A  (  —  )  ,  r  <  n  <  r. 

u(r)  \r  J 


(3.72) 


Now,  pick  r  >  r*  and  note  that 

1  fn\P  u(n ) 


r  -  /  \  />! 

)  )  (  =  /  UJx)  dx 

r  \r  /  u(r)  J0 


E  -  - 

f  ^  T*  \  T* 


(3.73) 


n=r*+ 1 

where  the  function  Ur  :  R+  — >  R  is  defined  by 


Ur(x)  :=  1 


x  >  — 
r 


Tr(xY 


u(rTr{x )) 


x  >  0 


u(r) 


with 


m  .  ,  rx]  +  1 
Tr(x)  :=  - ,  x  >  0. 


(3.74) 


(3,75) 


For  every  x  in  (0,1),  we  have  x  <  Tr(x)  <  1,  and  the  Uniform  Convergence 
Theorem  for  slowly  varying  functions  [7,  Theorem  1.2.1,  p.  6]  thus  implies 


The  pointwise  convergence 


>  =  1. 

r— >oo  U\r ) 


lim  Ur(x)  =  xp 

r—>oo 


(3.76) 


(3.77) 


is  now  an  immediate  consequence  of  (3.74)-(3.76).  Moreover,  making  use  once 
more  of  Potter’s  bound  (3.72)  we  find  that 


0  <  Ur(x)  <  ATr(x)p  s  <  4max(l,xp  5) 


with  hnite  integral 


max(l,  xp  s) 


dx  <  1  + 


1 

p-5  +  1 


(3.78) 
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under  the  choice  of  8.  These  remarks  together  with  (3.76)  lead  to 


lim  /  Ur{x)dx  =  /  xpdx  = 


1 


p+  1 


(3.79) 


by  dominated  convergence,  and  the  desired  limit  (3.71)  follows  by  going  to  the 
limit  in  (3.73). 

(ii)  The  proof  is  similar  to  that  for  Part  (i).  We  note  that  under  the  condition 
p  +  1  <  0  it  follows  from  (3.65)  that  (3.68)  is  equivalent  to 

1  fn\p  u{n )  1 


lim  V-  (- 

r-Von  *  r  \  r* 


r— >oc  z — '  T  \  T  / 
n—r 


(3.80) 


■  7  u(r)  p  +  1 

Thus,  pick  d  in  (0,  —1  —p),  and  A  >  0.  By  a  Potter  bound  [7,  p.  25]  similar  to 
(3.72)  there  exists  an  integer  r *  :=  r(8,  A)  such  that 

<5 


u{n) 
u{r ) 


.  fn\° 

<  A  —  ,  r  <  r  <  n, 

\r  ) 


(3,81) 


and  we  conclude  that  npu(n )  is  finite  for  each  r  =  1,2, . . .,  under  the  current 

n=r 

choice  of  d.  Next,  we  write 

1  /"n\P  u{n) 


E 

n=r-\- 1 


r  \r 


*/  u(r) 


Ur(x)  dx 


(3.82) 


where  the  function  Ur  :  R+  — >  R  is  now  defined  by 

Ur{x)  :=  Tr(x)p  u(rTrW)  x  >  o 
u(r ) 


(3.83) 


with  Tr(x)  as  in  (3.75).  Since,  for  every  x  in  [1,  oo),  we  have  x  <  Tr(x)  <  x  4*  1, 
the  Uniform  Convergence  Theorem  for  slowly  varying  functions  [7,  Theorem  1.2.1, 
p.  6]  implies  (3.76)  and  the  pointwise  convergence  (3.77)  again  follows.  Making 
use  of  Potter’s  bound  (3.81)  we  find  that 


0  <  Ur(x)  <  ATr(x)p+s  <  Axp+S,  r>r * 
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with  finite  integral 


1 


/  xp+s  dx  = - - - 

h  p+S+1 


(3.84) 


under  the  current  choice  of  S.  These  remarks  together  with  (3.77)  lead  to 

1 


lim 

r— too 


Ur{x)dx  =  /  xpdx  = 


(3.85) 


'1  J 1  P  +  1 

by  dominated  convergence,  and  the  desired  limit  (3.80)  follows  by  going  to  the 
limit  in  (3.82). 

(iii)  We  begin  by  noting  that  the  standard  asymptotics 


M  / , 

lim  V  -  =  In  (  - 

a;->oo  n  V  a 

n—  [ax] +1 


0  <  a  <  b 


imply 


[bx] 


lim  — 

t. — Vnn  •  ^  T) 


1  U(nl  =lnf- 


x->oo  * — <  nu(x)  \a 

n=[ax]+ 1 


0  <  a  <  b. 


(3.86) 


Indeed,  for  some  x**  in  R+,  it  holds  that 


a  <  sup  :  n  =  [ax]  +  1, . . . ,  [6.x]  j  <  b,  x  >  x* 


(3.87) 


By  the  Uniform  Convergence  Theorem  [7,  Theorem  1.2.1,  p.  6],  for  each  5  >  0 
there  exists  x*  :=  xa^{8)  >  x **  such  that 

u{tx) 


sup 


U(X) 


1 1  :  t  E  [a,  b]  \  <  8,  x  >  x*. 


(3.88) 


Combining  (3.87)  and  (3.88)  we  readily  get  (3.86)  in  the  form 


[bx] 


lim 

rr — ^ 


1  u(*x) 


^  (  \  =  ln  ~~ 
x^-oo  z '  n  u(x)  \  a 

n=[ax\-\- 1 

Now,  pick  e  in  (0, 1).  In  view  of  (3.86)  we  have 


,  r  u(x)  ,  v--v  1  u(n) 

lim  inf  — — —  >  lim  > — - 

x-j-oo  u(x)  x^oo  -f— '  nu(x ) 

'  '  n=[ex]+l  ' 


=  —  111  £. 


(3.89) 
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Since  e  can  be  chosen  arbitrarily  small,  it  follows  that 


u{x 

hm  — =  cx) 


x^oouiyX) 


(3.90) 


or  equivalently,  (3.70).  To  prove  that  u  is  slowly  varying,  pick  y  >  1  and  note  that 
for  every  x  >  0,  we  have 

[vx]  -t  (  \ 

\  \  1  u\n) 

u(yx)  =  u(x)  +  u(x)  > - , 

v  '  w  nu(x) 

n=[x\+ 1  v  ' 

so  that 

u(yx)  u(x)  1  u(n) 

u(x)  u(x)  ^  nu(x) 

v  '  v  '  n=[x]+ 1  V  ' 

It  is  now  straightforward  from  (3.70)  and  (3.86)  to  obtain 


lim  =  1. 


*->•00  U[X) 


The  case  y  <  1  is  handled  in  a  similar  way,  and  the  slow  variation  of  u  follows. 


A  proof  of  Proposition  3.6.1. 

From  (2.3)  it  always  holds  that 


E  [cr]  =  ^  P  [<r  >  n] 

n= 0 


E  [cr2]  1 

2E  [cr]  +  2' 


(3.91) 


We  consider  each  of  the  scalings  {Cr,  r  —  1,2,...}  associated  with  Theorems  3.3.1 
-  3.3.3,  separately: 

[Theorem  3.3.1]  Under  short-range  dependence,  we  have  E  [cr2]  <  oo,  and  it  is 
immediate  from  (3.91)  that  Condition  (B)  holds  for  the  choice  Cr  =  \/r,  r  =  1,  2, . . . 
(in  fact  for  any  choice  such  that  limCr  =  oo).  ■ 

r— too 


We  next  turn  to  Theorems  3.3.2  and  3.3.3.  Upon  substituting  (3.9)  (with  1  <  a  < 
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2)  into  (2.3),  we  readily  get  from  Proposition  3.8.1  (ii)  that 


P  [<r  >  n 


1 

E  [a 


^2p[a>  j] 

j=n 


[a 


1 

1)E  [a 


n1~ah{n) 


(n  — >  cx)), 


whence 


^  P  [a  >  n]  ~ 

n=  1 


(a 


1 

1)E  [a ' 


n1~ah(n) 

n=  1 


(r  — *  cx)) 


(3.92) 


provided  E  [<r]  is  infinite. 

[Theorem  3.3.2]  When  a  =  2  in  (3.9),  the  condition  E  [a2]  =  oo  implies  that 
E  [<t]  is  infinite  by  (3.91).  Thus,  (3.92)  holds  in  the  form 


Epp>»]~^E 


h(n) 


n= 1 


E  \a\  ^  n 

n—  1 


r  — >•  oc ) 


which,  from  Proposition  3.8.1  (iii)  is  seen  to  be  slowly  varying.  By  Proposition 
3.3.1,  the  scaling  {£r,  r  =  1,  2, . . .}  is  1/2-regularly  varying,  so  that 


1 

C 


^  P  [a  >  n] 

n=  1 


1 

— 

E  [a 


1 

h(r) 


E 

n=  1 


h[n ) 
n 


(r  — *  cx)) 


(3.93) 


for  some  slowly  varying  function  h  :  1R+  — »  1R+.  The  ratio  of  slowly  varying 
functions  being  itself  slowly  varying,  we  readily  conclude  from  Lemma  3.8.1  and 
(3.93)  that  Condition  (B)  holds.  ■ 


[Theorem  3.3.3]  On  the  range  1  <  a  <  2,  E  [cr2]  is  infinite,  and  so  is  E  [a]  by  virtue 
of  (3.91).  Proposition  3.8.1  (i)  applied  to  the  right-hand  side  of  the  asymptotic 
equivalence  (3.92)  yields 

E  P  l?  >  n)  ~  (2  —  "a)(a~  —  1)E  [<j]  (r  (3J4) 

By  Proposition  3.3.1  the  scaling  {£r,  r  =  1,2,...}  is  l/a-regularly  varying,  so 
that 


1 

Cr 


^  P  [a  >  n] 

n= 1 


rsj 


1  h(r)  r 2  “ 

(2  -  a) (a  -  1)E  [a]  h(r)  r« 


(r  — >  cx)) 


(3.95) 
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for  some  slowly  varying  function  h  :  1R+  — >  R+.  The  ratio  of  slowly  varying  func¬ 
tions  is  itself  slowly  varying,  and  Condition  (B)  is  now  a  direct  consequence  of 
Lemma  3.8.1  once  we  note  that  2  —  a  —  a-1  <  0  for  a  >  0.  ■ 
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Chapter  4 


Light  traffic  limits 

4.1  Introduction 

We  seek  to  characterize  the  light  traffic  limiting  behavior  of  the  queueing  system 
with  M\G\oo  inputs.  Our  main  tool  for  accomplishing  this  task  is  a  methodology 
presented  in  a  series  of  papers  by  Reiman  and  Simon  [49,  50,  51].  Their  approach 
provides  a  general  framework  for  deriving  asymptotic  results  in  systems  where  the 
quantity  of  interest  can  be  expressed  as  a  function  of  a  Poisson-like  process.  If 
A  >  0  denotes  the  intensity  of  the  Poisson  process  driving  the  system,  then  the  light 
traffic  information  furnished  by  the  Reiman-Simon  technique  consists  of  derivatives 
of  the  quantity  of  interest,  with  respect  to  A,  evaluated  at  A  =  0+.  We  devote 
most  of  our  efforts  to  the  case  c  =  1,  because  then  additional  expressions  become 
available  by  relating  the  system  with  M\G\oo  inputs  to  that  with  instantaneous 
inputs.  The  results  quantify  the  effect  of  the  session  duration  distribution  G  and 
reveal  the  differences  between  the  gradual  M\G\oo  inputs  and  the  instantaneous 
arrivals  of  a  classical  GI\GI\1  queue.  Moreover,  in  the  special  case  c  —  1,  they 
suggest  a  classification  of  the  light  traffic  behavior  of  the  multiplexer  buffer  content 
in  terms  of  the  short-  vs  long-  range  dependent  property  of  the  M\G\oo  arrival 
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process. 


4.2  Reiman— Simon  theory 

4.2.1  Preliminaries 

To  place  our  problem  in  the  context  of  the  Reiman-Simon  methodology,  we  start 
by  introducing  bi-infinite  counterparts  to  the  sequences  of  IN  valued  rvs  repre¬ 
senting  session  arrivals  and  their  activity  durations.  That  is,  we  consider  the 
bi-infinite  sequences  of  W-valued  rvs  {/?n,  n  =  0,±1,±2, ...}  and  { anj ,  n  = 
0,  ±1,  ±2, . . . ;  j  =  1,2,...}  which  are  constructed  on  a  common  probability  space 
(Q,JT, PA),  are  mutually  independent  and  satisfy  Assumptions  (A)(i)  and  (ii),  re¬ 
spectively.  We  also  introduce  the  bi-infinite  sequences  of  M\G\oo  inputs  {bn,  n  = 
0,  ±1,  ±2, . . .}  and  queue  lengths  {qn,  n  =  0,  ±1,  ±2, . . .},  where  the  former  are 
given  by 

n  Pk 

bn-=  y:  yi  [<Tk,j  >  n  —  k] ,  n  =  0,±1, ..., 

k=— oo  j= 1 

so  that  bn  denotes  the  number  of  active  sessions  at  the  beginning  of  slot  [n,  n  +  1). 
In  this  setup,  instead  of  starting  the  Lindley  recursion  at  time  n  =  0  we  take  the 
viewpoint  that  the  system  has  been  operating  from  time  n  =  —  oo,  i.e., 

qn+i  =  [qn  +  bn+1  -  c}  + ,  n  =  0,±1,....  (4.1) 

Under  the  stability  condition  AE  [er]  <  c  convergence  to  the  stationary  ]R+-valued 
rv  has  already  taken  place  by  time  n  =  0,  that  is, 

qo  st  qoo •  (4.2) 

Application  of  the  Reiman-Simon  method  entails  conditioning  on  the  number 
of  arriving  sessions  and  their  corresponding  activity  durations.  We  introduce  the 
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necessary  notation.  First,  we  denote  by  the  empty  set  0  the  event  that,  from  time 
t  =  — oo  to  t  —  +oo,  there  are  no  session  arrivals  at  all.  Next,  £x  n  =  1,2, . . ., 
and  consider  the  following  event  u:  There  are  exactly  n  sessions  that  ever  become 
active.  For  each  i  =  1,  2, . . . ,  n,  fixed  ti  =  0,  ±1, . . .,  and  ki  =  1,2,.. .,  the  ith 
session  becomes  active  during  time  slot  time  [tt  —  l,  t,)  and  its  activity  period  lasts 
ki  time  slots.  We  denote  such  an  event  uj  by  {ti, . , . ,  tn;  ki, . . . ,  kn}.  In  other  words, 
{ti, . . . ,  tn\  ki, . . . ,  kn}  corresponds  to  an  event  where  n  sessions  arrive  to  the  system 
over  all  time,  and  these  arrivals  occur  during  time  slots  [ti  —  1,  ti), . . . ,  [tn  —  1,  tn ) 
with  respective  activity  durations  ki,...,kn.  We  drop  the  duration  indices  to 
denote  unions  of  events,  i.e., 

OO  OO 

{ti,...,tn]}  :=  |J  ...  1J  (4.3) 

k]_=  1  kn~  1 

is  the  event  that  exactly  n  sessions  arrive  to  the  system  over  all  time,  and  these 
arrivals  occur  during  time  slots  [ti  —  1,  ti), . . . ,  [tn  —  1,  tn). 

4.2.2  Light  traffic  derivatives 

Let  the  generic  system  performance  metric  0(A)  be  expressed  as 

0(A)  =  J  0  dP a  (4.4) 

for  a  suitably  chosen  rv  0  :  fl  — >  H.  For  example,  'tp(oj)  can  be  chosen  as  the  queue 
length  at  time  t  =  0,  corresponding  to  a  sample  path  u  in  12,  in  which  case,  from 
(4.2)  and  (4.4),  the  performance  metric  0(A)  is  the  expected  value  EA  [q^]  of  the 
stationary  queue  length  q^. 

Following  the  Reiman-Simon  method,  we  decompose  the  expectation  in  (4.4) 
according  to  occurrences  of  session  arrivals/activity  durations  events  of  the  form 
{ti, . . . ,  tn]  ki, . . . ,  kn}.  To  do  this,  for  each  n  =  1,  2, ...  we  associate  with  0  several 
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auxiliary  functions.  First,  the  expected  value  -0  of  0,  conditional  on  the  session 
arrivals  event  {ti, ... ,  tn:  }  is  given  by 

0({ti,  •  ■ , ,  t„;  })  :=  E  [0  |  {tu  . . . ,  tn]  }] ;  (4.5) 

this  does  not  depend  on  A.  Next,  we  dehne  the  function  0  :  {1,  2,  R  by 

+oo  +oo 

0(/ci, . . . ,  kn)  :=  Y  ...  Y  ip({ti,...tn;ki,...,kn}).  (4.6) 

t\ —  OO  tn —  OO 

where  0({fi, . . .  tn\  ki, . . . ,  kn})  is  the  value  of  -0  when  n  sessions  arrive  to  the  sys¬ 
tem  over  all  time,  these  session  arrivals  occur  during  time  slots  [fi  —  1,  ti), . . . ,  [tn  — 
1,  tn),  and  their  respective  activity  durations  are  given  by  ki,...,kn.  Further¬ 
more,  let  nnJ  denote  the  set  of  unordered  j- tuples  chosen  from  {1,2,...,  n)  (with 
repetitions  allowed),  where  for  7 r  =  {0,  i2, . . . ,  ij}  in  II nj,  we  use  the  notation 

f  7T  •  {^ll  )  02  )  •  •  •  J  )  }■ 

Now,  for  any  given  arrival  event  {tl5 . . . ,  tn\  },  we  dehne 

n 

t„;»  Y,  AU).  (4.7) 

j= o  Tren  nj 

For  instance,  we  have 

^({0})  =  0({0})  -0(0),  (4-8) 

^({0,0;})  =  0({0,0;})  -0({0;})  ~$({h;})  +  0(0)  (4.9) 

and  so  on. 

The  formulas  for  the  light  traffic  derivatives  can  be  obtained  by  considering  a 
system  where  only  arrivals  in  an  interval  of  the  form  [— T,  T),  for  T  =  1,2,...,  are 
ever  allowed  to  enter;  let  0T(A)  be  the  corresponding  performance  metric.  The  idea 
is  to  calculate  first  the  derivatives  of  0T( A)  with  respect  to  A,  at  A  =  0+,  and  then 
let  T  go  to  infinity.  Clearly,  it  is  necessary  to  justify  that  this  interchange  of  limits 
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ill  A  and  T  leads  to  the  correct  answer.  To  that  end  we  enforce  an  assumption  on 
the  finiteness  of  the  exponential  moment  of  a: 

Assumption  (C)  There  exists  9*  >  0  such  that  E  [ee><T]  <  oo  for  9  <  9* . 

In  [51]  it  is  shown  that  under  Assumption  (C)  the  interchange  of  limits  is  indeed 
valid;  here  we  simply  restate  this  conclusion  as 

Proposition  4.2.1  Under  Assumption  (C)  it  holds  that 

lim  0T(O+)  =  0(0+),  n  —  0, 1, - 

T^ood\n  v  ’  d\n  v  ’ 

The  following  result  is  essentially  a  discrete-time  version  of  Theorem  2  in  [51,  p. 
30],  and  enables  us  to  calculate  the  nth  order  derivative  of  0( A)  at  A  =  0+  by 
considering  scenarios  where  at  most  n  sessions  ever  arrive  to  the  system. 


Proposition  4.2.2  If  Assumption  (C)  is  satisfied,  then 


lim  0(A)  =  0(0), 
a^o+  v  '  v  ' 


(4.10) 


and  for  each  n  =  1,2,...,  it  holds  that 

+oo  +oo 


dn 

d\r‘ 


0(o+)—  ...  T({fi, . . . , t„; }). 


(4.11) 


t~\_ -  OO  1 77. —  OO 


Proof.  For  each  T  =  1,2,...  and  j  =  0, 1, . . .  let 

pj  (a)  := 


denote  the  probability  that  j  discrete-time  Poisson  session  arrivals  occur  during 
the  interval  [— T,  T).  For  each  n  =  0.1,...,  the  nlh  derivative  of  P-  (A)  with  respect 
to  A  is  given  by 


dn 

d\n 


min  (n,j) 

Pj  (A)  =  (2T)n  (-1) 

i= 0 
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so  that 


rln 

w  p>(0+)  = 


(2 T)"(  "  )  (-l)”--’  if  n>j 

3 


(4.12) 

0  if  n  <  j. 

Given  that  j  Poisson  arrivals  have  occured  in  -T.  T).  they  are  uniformly  dis¬ 
tributed  over  the  2 T  time  slots.  Thus, 

oo  T  T  1 

=  Y  •••  5Z  ••>**;}) 

j=0  t1=-T+l  t,=-T+l  '  ' 

and  using  (4.12)  we  get 


dn 

dX1 


(0+)  =  £>T)”  (  "  )  £  ...  J2 


3=0 

n 


ti=-T+l  tj  =  -T+ 1 

T 


E<  2T)"-1(-1)"-1  V  V  £(t„) 

J— 0  7rElInj  t7r=— T+l 

n  T  T 

E)-1)"^  E  E  ■  ■  ■  E 

j — o  ti —  r+i  tn — — t*+i 

J]  ...  ]T 

^  l — — r+l  tn=  T+l 


(4.13) 


Letting  T  go  to  infinity  in  (4.13)  and  invoking  Proposition  4.2.1  we  conclude  that 
(4.11)  holds  true.  ■ 

We  rely  on  Proposition  4.2.2  to  calculate  light  traffic  derivatives  of  system  quan¬ 
tities  in  the  sequel. 


4.2.3  Case  c  >  1 

We  now  consider  a  Lindley  recursion  (4.1)  with  release  rate  c  >  1.  Fix  some  integer 
p  =  1,2,....  Take  0  :=  Qop,  where  go  is  the  queue  length  at  time  n  =  0,  so  that 
the  performance  measure  of  interest  is  0(A)  =  [g^],  the  pth  moment  of  the 
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stationary  queue  size.  To  determine  its  light  traffic  derivatives  we  need  to  evaluate 
the  quantities  appearing  in  Proposition  4.2.2. 

Clearly,  if  at  most  [cj  sessions  are  ever  active,  then,  since  each  active  session 
generates  one  arrival  per  time  slot,  their  inputs  are  flushed  out  of  the  queue  by  the 
end  of  the  time  slot  and  the  queue  remains  empty.  So,  for  each  m  =  1,  2, . . . ,  |_cj 
and  ti, . . . ,  tm  =  0,  ±1,  ±2, ...  it  is  immediate  from  definition  (4.5)  that 

•0(0)  =  0  and  fi({tu . . . ,  tm;  })  =  0.  (4.14) 

Next,  fix  some  b  >  0  and  take  ip  \=  1  [q0  >  b\,  in  which  case  the  performance 
metric  of  interest  is  the  tail  probability,  0(A)  =  Ea  [ip]  =  Pa  [g<x>  >  b\.  If  at  most 
|_cj  sessions  become  active,  then  the  same  simple  considerations  as  before  apply, 
because  whenever  g0  =  0  we  also  have  1  [g0  >  b]  =  0.  Thus  relations  (4.14)  still 
hold  true  for  the  current  choice  ip  :=  1  [go  >  b]  as  well.  Consequently,  by  combining 
(4.14)  with  Proposition  4.2.2  for  each  of  the  functions  ip  :=  qop  and  ip  1  [go  >  b\ 
we  arrive  at 


Proposition  4.2.3  Consider  a  Lindley  recursion  (4-1)  with  release  rate  c  >  1.  If 
Assumption  (C)  is  satisfied  then  for  each  m  —  1, 2, . . . ,  |_cj,  it  holds  that: 

(a)  For  p  =  1,2,..., 


lim  Ea  [go/]  =  0 

A— >0+ 


=  0  and 


d\r 


Ea  [gQ 


=  0. 


A=0+ 


(4,15) 


(b)  For  b  >  0, 


dm 

lim  PA  [goo  >  b]  =  0  and  ——  PA  [g«>  >  b] 
a->o+  a\ni 


=  0. 


A=0+ 


(4-16) 


Before  proceeding  we  make  a  few  comments.  Note  that  when  c  =  1  the  server 
can  only  process  one  session  at  a  time,  so  that  the  system  can  be  viewed  as  a  single 
server  queue  operating  in  discrete  time.  Then  from  (4.16)  with  c  =  1  it  is  already 
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apparent  that  the  light  traffic  limits  of  a  queueing  system  with  M\G\oo  inputs 
differ  from  those  of  a  standard  single  server  GI\GI\1  queue.  For  the  system  first 
derivative  is  here  zero,  while  in  the  system  (2.17)  with  instantaneous  inputs  the 
first  light  traffic  derivative  is  positive.  This  is  a  manifestation  of  the  fact  that  work 
that  joins  the  system  gradually,  as  is  the  case  with  M\G\oo  inputs,  generates  less 
queueing  that  work  arriving  instantaneously.  Also,  relation  (4.16)  reflects  (though 
in  a  rough  manner)  the  statistical  multiplexing  gain:  All  powers  of  A  up  to  and 
including  A^  offer  no  contribution  to  the  tail  probability.  Thus  (4.16)  implies  that, 
in  light  traffic,  increasing  the  multiplexer  release  rate  c  while  maintaining  the  same 
system  utilization  AE  [cr]  /c  results  in  a  decreasing  tail  probability  Pa  [goo  >  b]. 

In  view  of  (4.15)  and  (4.16)  the  focus  shifts  to  the  calculation  of  the  derivative 
of  order  |_cj  +1.  This  is  the  first  non- zero  derivative  and  it  is  clearly  more  infor¬ 
mative  than  the  |_cj  +  1  lower  order  derivatives,  for  it  provides  the  leading  term  in 
expansions  of  system  quantities  around  A  =  0. 

4.2.4  Case  c  =  1 

In  this  section  we  consider  a  Lindley  recursion  (4.1)  with  release  rate  c  =  1.  This 
corresponds  to  the  situation  where  each  active  session  in  the  M\G\oo  input  gener¬ 
ates  arrivals  at  rate  equal  to  the  multiplexer  service  rate.  In  this  case  we  are  able 
to  obtain  explicit  expressions  for  the  second  order  light  traffic  derivatives  of  sys¬ 
tem  quantities  by  carrying  out  in  full  the  calculations  associated  with  Proposition 
4.2.2. 

We  take  ip  to  be  either  ip  :=  q0p  for  some  p  =  1,  2, . . .,  or  ip  1  [q0  >  b]  for 
b  >  0,  yielding  <p( A)  =  EA  [ q^ ]  and  <p{ A)  =  PA  [goo  >  b\,  respectively.  From  (4.11) 
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it  follows  that 


j2  +°°  +°° 

_0(o+)=  22  J2  ?(Oi>*2;}),  (4-17) 

t1=  — oo  t2  =  —  00 

where  we  have  also  taken  into  account  (4.14).  We  thus  need  to  evaluate  h ;  })• 

This  calculation  requires  consideration  of  a  system  where  only  two  sessions  are  ever 
active.  In  particular,  since 

OO  OO 

${{ti,t2;})  =^2^2  ^({ti,t2]kl,k2})P  [a  =  h]P  [a  =  k2]  ,  (4.18) 

k±=l  &2=1 

we  need  only  examine  the  queue  length  process  induced  by  events  of  the  form 
1 ,  f2;  k\ ,  k2 } .  It  is  convenient  to  interchange  the  order  of  summations  appearing 
in  (4.17)  and  (4.18)  -  this  can  be  done  without  qualms  for  the  summands  are  all 
non- negative.  Recalling  definition  (4.6),  we  calculate  the  sums  over  the  arrival 
times  first,  say 

+OO  +OO 

2{kX)k2)=  22  '•P({ti,t2;k1,k2}),  (4.19) 

tl  =  —  OO  t2  =  —  OC 

and  then  rewrite  (4.17)  as 

A^0(O+)  =  E  ^(<71,02)  ,  (4.20) 

where  <j\  and  n2  are  i.i.d.  copies  of  the  generic  activity  duration  rv  er.  Thus,  the 
calculation  of  the  second  derivative  is  reduced  to  the  evaluation  of  -?/>.  To  deter¬ 
mine  '0  consider  a  session  arrival/activity  duration  event  of  the  form  {t\,  t2 ;  ki,  k2}. 
Observe  that  if 

min(f2  +  k2,ti  +  ki)  >  max(f1;  t2)  (4-21) 

the  two  sessions  that  arrive  in  slots  [ti  —  l,ti)  and  [t2  —  1  ,t2)  are  simultaneously 
active  from  time  max(t1;  t2)  until  min(f2  +  k2,  t\  +  ki).  In  that  case  the  queue  size 
evolves  as  shown  in  Figure  4.1;  otherwise,  i.e. ,  if  min(t2  +  k2,  ti  +  ki )  <  max(f1;  t2), 
it  is  identically  zero. 
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max(t1  ,t2)  min^+k^+k,,)  max^+k^+k,,)  min^.y+k^ 

Time 


Figure  4.1:  Queue  length  evolution  under  the  event  {t\,  t2;  k\ ■  k2} 


By  inspection,  under  (4.21),  the  queue  length  at  time  t  =  0  is  given  by 


-max(£i,£2), 

min(£i  +  £q,  t2  +  k2 ) 

qodti,^^,^})  =  <  —  max(£i,£2), 

min(£i,£2)  +  ki  +  k2: 

0, 


if  max(£i,  t2)  <  0  and 
0  <  min(£i  +  £q,  t2  +  k2 ) 
if  min(£1  +  kXlt2  +  k2)  <0 
and  0  <  max(£i  +  ki,  t2  +  k2 ) 
if  max(h  +  ki,  t2  +  k2)  <0 
and  0  <  min(£1,  t2)  +  ki  +  k2 
otherwise. 


To  facilitate  the  evaluation  of  0  we  display  qo({t\,  t2,  ki,  k2})  in  the  M2-plane;  the 
values  for  k\  >  k2  and  k±  <  k2  are  shown  in  Figures  4.2  and  4.3,  respectively. 


We  carry  out  the  detailed  calculation  of  the  function  -0  corresponding  to  each 
of  the  choices  0  :=  q0,  0  :=  q0 2  and  0  :=  1  [q0  >  b],  b  =  0, 1, . . .,  in  Section  4.4. 
The  results  are  summarized  in  Lemma  4.2.1  below. 
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-<Vk2>  V  -k2  0 

Time  t 

Figure  4.2:  Values  of  qo({ti,t2,  ki,  k2})  when  k\  >  k2 


Figure  4.3:  Values  of  qo({ti,t2,  ki,  k2})  when  k\  <  h 


Lemma  4.2.1  Assume  c=l  in  the  Lindley  recursion  (4-1  )■  The  following  state¬ 
ments  hold: 

(a)  If  if  :=  q o  then 

^  k1k2(k1  +  k2),  h,k2  =  1,2, ....  (4.22) 

(b)  If  if  :=  q0 2  then 

k2)  =  ^  kik2(l  +  kik2),  ku  k2  =  1,  2, . . . .  (4.23) 

(c)  If  if  \=  1  [g0  >  6]  /or  some  6  =  0,1,...,  then 

*f{ki,  k2)  =  ^  {(/ci  -  b)2  +  (k2  -  b)2J  -  ^  (fci  +  /c2  -  26) 

+2(Aq  —  b){k2  —  6)  +  1  (4.24) 

for  ki,  k2  >  6  and  if(ki,  k2)  =  0  otherwise. 

In  principle  it  is  possible  to  evaluate  the  function  if  corresponding  to  if  :=  q0p  for 
p  >  2  in  a  similar  manner.  However,  such  calculations  become  increasingly  tedious 
for  large  p. 

We  now  obtain  explicit  expressions  for  the  second  order  light  traffic  derivatives 
of  system  quantities  by  combining  each  one  of  (4.22),  (4.23)  and  (4.24)  with  (4.20). 
This  leads  to  the  following 


Proposition  4.2.4  Let  c  =  1  in  the  Lindley  recursion  (4-1).  Linder  Assumption 
( C )  it  holds  that: 

(a)  The  moments  EA  [q^]  and  EA  [goo2]  satisfy 


d 2 

d\2 


Ea  loo 


A=0+ 


E  [cr]  E[cr2 


(4.25) 


and 


d 2 

d\2 


Ea  [goo2_ 


A=0+ 


\  (E  W2  +  EM2)  . 


(4.26) 
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(b)  For  each  b  =  0, 1, . . . 


d 2 

dX 2 


Pa  [goo  >  b\ 


A=0+ 


E 


{a  -  b)+2 1  P  [a  >  b]  +  2  E  [(a  -  6)+]: 


-3  E  [(<r  -  &)+]  P  [a  >  b]  +  P  [a  >  bf  .  (4.27) 


Proposition  4.2.4  delineates  a  light  traffic  behavior  for  the  queue  with  M\G\oo 
arrivals  that  is  certainly  different  from  that  of  a  classical  GI\GI\1  queue.  Indeed, 
when  considering  the  first  two  terms  in  a  light  traffic  expansion  of  PA  [gi0O  >  b] 
around  A  =  0  (for  c  =  1)  we  see  that  the  first  derivative  (4.16)  is  zero  so  that  the 
second  derivative  is  the  most  informative.  This  is  given  by  (4.27)  which  highlights 
the  role  of  the  activity  duration  rv  a,  through  both  its  distribution  and  its  first  two 
moments.  Thus,  a  light  traffic  expansion  of  the  tail  probability  Pa  [goo  >  b]  induced 
by  M\G |oo  arrivals  is  completely  different  from  the  corresponding  expansion  for 
the  system  with  instantaneous  inputs:  This  can  be  obtained  from  (2.38)  (or  via  the 
Reiman-Simon  method)  and  is  given  by  Pa  [q^  >  b\  rsj  AE  [a]  P  [a  >  b  +  1]  (A  -> 
0),  thus  starting  with  a  non-zero  first  order  term  A.  Notice  also  that  here,  even  if 
Assumption  (C)  were  to  be  relaxed,  (4.27)  shows  that  for  P\  [goo  >  b]  to  decay  like 
A2  for  small  A  it  is  necessary  that  E  [<r2]  be  finite.  If  E  [a2]  =  oo,  as  is  the  case 
for  long-range  dependent  M\G\oo  arrivals,  expression  (4.27)  yields  infinity  and  A2 
is  no  longer  the  correct  order  of  decay.  A  different,  smaller  exponent  should  be 
sought  in  the  long-range  dependent  case. 


4.2.5  A  heavy— light  traffic  relationship 

In  Sections  4.2.3  and  4.2.4  we  have,  at  least  partially,  mapped  out  the  light  traffic 
behavior  of  the  multiplexer  with  M\G\oo  inputs.  Of  course  this  partial  information 
is  augmented  whenever  expressions  for  the  next  higher  order  derivatives  become 
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available.  Unfortunately,  the  calculations  soon  become  intractable  and,  typically, 
explicit  expressions  are  available  only  for  the  first  non-zero  derivative.  In  principle 
however,  if  light  traffic  derivatives  of  every  order  were  known,  then  system  quan¬ 
tities  would  be  completely  determined  away  from  the  light  traffic  regime  by  their 
Taylor  series  expansion  (under  analyticity  assumptions).  In  particular,  knowledge 
of  all  light  traffic  derivatives  would  also  imply  full  knowledge  of  system  quantities 
in  heavy  traffic.  This  observation  raises  the  question  as  to  how  these  successive 
derivatives  at  A  =  0+  are  related  to  the  respective  limiting  behavior  of  system 
quantities  in  the  heavy  traffic  regime.  The  answer  is  given  in  [56],  where  a  sim¬ 
ple  yet  rather  unexpected  relationship  is  established.  This  relationship  suggests 
a  method  for  constructing  certain  approximations  of  system  quantities;  these  are 
precisely  the  interpolation  approximations  discussed  in  Chapter  5. 

The  link  between  heavy  traffic  limits  and  light  traffic  derivatives  is  provided  by 
the  following  proposition,  which  is  a  special  case  of  Simon’s  results  [56]. 


Proposition  4.2.5  Let  the  function  0  :  [0, 1)  — »  IR+  be  analytic  on  [0, 1).  Suppose 
that  the  limit 

H  :=  lim  (1  —  x)4>{x)  (4.28) 

x — >1 — 

exists  and  is  positive  and  finite.  Then 


1  dn 

lim  —  - —  0(0+)  =  H. 

n->oon!  dxn 


(4.29) 


We  apply  Proposition  4.2.5  in  the  context  of  the  Lindley  recursion  (4.1)  with  c  =  1, 
describing  a  system  with  M\G\oo  arrivals,  by  setting  x  :=  AE  [cr],  so  that  x  lies  in 
[0, 1)  whenever  the  system  is  stable.  Let  us  assume  0(A)  :=  EA  [g^]  is  analytic  in 
[0, 1).  From  Assumption  (B)  and  Theorems  3.3.1  and  3.4.1  we  infer  that  the  limit 

H  :=  ,  >,+  A  -  AE  H>  E»  fed  =  IB  (4-3°) 

A— >1/E[(t]  zrLi  [G 
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is  positive  and  finite  whenever  E  [<v2]  <  oo.  In  that  case  Proposition  4.2.5  reads 

1  1  dn  E  [a2] 

n^oon!  E  [a]n  d\n  A  ^°°  a=o+  2E  [<t] 

Recalling  Proposition  4.2.4(a),  which  provides  the  second  light  traffic  derivative, 
we  observe  that  (4.31)  already  holds  as  an  equality  for  n  =  2,  and  not  just  in  the 
limit  as  n  — »  oo.  This  hints  to  the  possibility  that  (4.31)  holds  with  equality  for 
all  n  =  2,3...,  namely 

^EV9<j|A_o+  =  YEM”~‘EK]'  n  —  2,3 ... .  (4.32) 

If  (4.32)  were  to  hold,  or  were  assumed  to  hold,  then  additional  results  would  be 
within  reach.  In  fact,  it  follows  from  (2.33)  that  (4.32)  indeed  holds  for  the  expected 
queue  size  Ea  [g^]  in  the  system  with  instantaneous  inputs.  These  observations 
motivate  a  closer  examination  of  the  relation  between  the  discrete-time  queue 
driven  by  M\G\oo  inputs  and  the  corresponding  system  with  i.i.d.  instantaneous 
inputs. 

4.3  Gradual  inputs 

In  this  section  we  revisit  the  system  with  gradual  input  process.  We  aim  at  obtain¬ 
ing  expressions  for  the  expected  queue  length  and  the  probability  that  the  queue 
length  is  zero.  We  adapt  the  arguments  of  [55]  and  suppose  that  the  system  with 
gradual  inputs  operates  in  parallel  with  a  second  system,  with  instantaneous  in¬ 
puts.  That  is,  we  consider  the  Lindley  recursion  (2.17)  with  instantaneous  inputs 
(2.16),  and  the  sequence  {q„\  n  —  0, 1, . . .}  evolving  according  to 

=  0;  Qnh  =  [Qn]  +  t>nh  ~  c]  +  ,  71  =  0,  1,  ...  ,  (4.33) 
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(4.34) 


where  the  sequence  of  gradual  inputs  {b^\  n  =  1.2,...}  is  given  by 

n  Pk 

bn]  =  1  [^4  >  n  ~  k\  >  n  =  0, 1, . . . , 

k= 1  j= 1 

with  the  convention  that  empty  sums  are  zero.  We  couple  the  two  input  sequences 
(4.34)  and  (2.16)  by  constructing  them  both  from  the  same  i.i.d.  rvs  {/3n+ 1,  n  = 
0, 1, . . .}  and  {oVi+M,  n  —  0, 1, ,  i  =  1,  2, . . .}  of  Section  2.1.1. 

4.3.1  Stationary  version 

The  gradual  input  sequence  {b^\  n  =  0, 1, . . .}  of  (4.34)  driving  recursion  (4.33) 
is  not  the  stationary  version  of  the  M\G\oo  busy  server  process.  This  stationary 
version  {bn,  n  =  0,1,...}  was  given  by  (2.1)  in  Section  2.1.1  and  contains  an 
additional  term  due  to  the  servers  initially  busy.  It  is  thus  necessary  to  establish 
a  relation  between  the  stationary  regimes  for  (4.33)  (with  non-stationary  inputs) 
and  for  (2.11)  (with  stationary  inputs).  This  is  done  in  the  following  Lemma. 

Lemma  4.3.1  Under  the  stability  condition  AE  [a]  <  c  the  sequences  { qn ,  n  = 
0, 1, . . .}  and  {q„\  n  =  0, 1, . . .}  associated  with  the  Lindley  recursions  (2.11)  and 
(4-33),  respectively,  converge  weakly  to  the  same  stationary  rv  q^. 

Proof.  Define  M  :=  max{cr0j-,  j  =  1, . . . ,  b},  where  b  and  {croj,  3  =  1, . . . ,  b} 
are  given  in  Assumption  (A),  so  that  M  is  the  (finite)  time  that  elapses  until  the 
activity  periods  of  all  initially  active  sessions  expire.  From  (2.1)  and  (2.2)  we  note 
that 

bn+M+l  —  ^n+M+1 1  U  =  0,  1,  .  .  .  .  (4.35) 

Thus,  on  every  sample  path,  the  queue  size  sequences  { qn+M ,  n  =  0, 1, . . .}  and 
{in+M,  n  =  0, 1)  •  •  •}  are  driven  by  the  same  input  sequence  (4.35)  with  possibly 
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different  (in  general)  initial  conditions  qM  and  ■  Because 

1  n 

lim  —  bf'  =  AE  [a]  a.s. 

n— >-+oo  77,  '  ^  J 
3= 1 

we  can  invoke  Lemma  6.1.4  in  [54,  p.  134],  to  conclude  that  under  the  stability 
condition  AE  [er]  <  c  the  sequences  {qn+M,  n  =  0, 1, . . .}  and  {g^+M,  n  —  0, 1, . . .} 
strongly  couple  on  every  sample  path,  hence,  due  to  the  already  established  fact 
that  qn  =^n  goo,  weakly  converge  to  the  same  stationary  rv  g^.  ■ 

4.3.2  Case  c  <  1:  A  stochastic  comparison 

In  this  section  we  assume  that  the  multiplexer  release  rate  is  c  <  1,  and  establish 
a  strong  stochastic  comparison  between  q ^  and  qj£.  For  each  n  =  1,2,...  let 
dn  denote  the  total  amount  of  work  processed  by  the  server  in  the  system  (2.17) 
with  instantaneous  inputs  (2.16)  during  the  interval  (0,  n\.  The  sequence  {dn,  n  = 
1,2,.. .}  satisfies 

n 

dn  =  ^2uk-  Q  n\  n  =  1,2, ... .  (4.36) 

k= i 

Next,  look  at  the  sibling  system  (4.33)  with  gradual  inputs.  Under  the  condition 
c  <  1,  a  single  active  session  suffices  to  make  full  use  of  the  server  capacity  c, 
despite  the  gradual  nature  of  its  input.  This  is  no  longer  true  if  c  >  1;  in  that  case 
if  there  is  only  one  active  session  in  the  system  it  is  served  at  unit  rate  and  the 
portion  c  —  1  of  the  capacity  remains  unused.  Thus,  under  the  assumption  c  <  1, 
in  the  system  (4.33)  (with  gradual  inputs)  the  server  processes,  in  every  sample 
path,  exactly  the  same  amount  of  work  as  the  server  in  the  coupled  system  (2.16) 
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(4.37) 


(with  instantaneous  inputs)  .  This  crucial  observation  enables  us  to  write 

n 

dn  =  ^b[a)  ~  Qn\  n  —  1,2,  ...  . 
k= 1 

for  the  same  completed  work  sequence  {dn,  n  =  1,2, . . .}  as  that  of  (4.36).  We 
can  thus  obtain  a  relation  between  the  queue  sizes  in  the  two  systems:  For  each 
n  =  1,2,...,  set 

n  0n-k+ 1 

Vn  :=  E  E  (^n-fc+ij  -  &)  +  (4.38) 

fc=l  j= 1 

and  note  from  relations  (4.36)  and  (4.37),  and  definitions  (4.34)  and  (2.16)  that 

n  0k 
fc=i  i=i 

n  0n-k+ 1 

=  E  E  K-fe+ij  -  A;)+ 

fc=i  j=i 

=  u„.  (4.39) 

Observe  here  that  (4.39)  implies  the  sample  path  inequality 

q(n]  <  q(n\  n  =  1,2,....  (4.40) 


Take  n  going  to  infinity  in  (4.40)  and  recall  that  the  stochastic  ordering  <st  is 
stable  under  weak  convergence  (Proposition  B.3  of  Appendix  B).  It  is  now  plain 
from  Lemma  4.3.1  that 

<„  «£’  (4.41) 

with  q ^  =st  goo,  and  the  following  stochastic  comparison  is  obtained: 

Proposition  4.3.1  Consider  the  Lindley  recursions  (2.11)  and  (2.17)  with  inputs 
characterized  by  a  common  pair  (A,  a)  and  release  rate  c  such  that  AE  [a\  <  c  <  1. 
The  respective  stationary  rvs  q ^  and  q£)  satisfy 

Qoc  <st  q^-  (4.42) 


87 


4.3.3  Case  c  <  1:  Expected  queue  size  E  [q^] 

In  addition  to  the  stability  condition  AE  [a]  <  c  we  assume  that  E  [ex2]  <  oo, 
which  ensures  E  \q^]  <  oo,  and  proceed  again  via  (4.39)  and  Lemma  4.3.1  to 
express  E  [q^]  in  terms  of  E  [g^] .  For  each  n  =  1,  2, ....  it  holds  that 

E  [,<“>]  =  E  [«£•>]  -  E  [„„] .  (4.43) 

From  (4.38)  by  Wald’s  identity  we  have 

n 

E[vn]  =  Aj]E  [(v-k)+] 

k= 1 

n 

=  AE  [a]  ^  P  [a  >k] 

k= 1 

whence 


lim  E  [vn\ 

n— >•  oo 


AE  [a]  (E  [a]  —  1) 


AE  [a]  ^ 
^E  [a(a 


E  [cr(cr  +  1)] 
2E  [a] 

-1)] 


(4.44) 


where  the  second  equality  follows  from  (2.4).  Turning  to  the  sequence  {q^\  n  = 
0,1,.. .}  we  note  that  (2.17)  is  driven  by  the  i.i.d.  rvs  {un,  n  =  1,  2, . . .}  of  (2.16) 
and  satisfies  q ^  =  0  <  q[u\  so  that  the  monotonicity  result  of  Proposition  B.7 
(Appendix  B)  applies.  We  get 


Au)  < 


Sstq£\ 


n  =  1,2,... 


(4.45) 


and  because  of  (4.40)  we  also  obtain 


/y(a)  /y(W)  71  =.  1  2 

S/n  ioo  ?  ^  J-j  . 


(4.46) 


Since  E  [,£>]  <  oo  whenever  E  [cr2]  <  oo,  inequalities  (4.45)  and  (4.46)  imply 
that 

E  [<&)]  =  lim  E  [g£*>]  and  E  [g£>]  =  lim  E  [g£*>]  (4.47) 

L  J  n— >-oo  L  J  J  n— too 

by  the  dominated  convergence  theorem.  Therefore,  letting  n  go  to  infinity  in  (4.43) 
and  using  (4.44),  (4.47)  and  Lemma  4.3.1  we  collect 


Proposition  4.3.2  Consider  the  Lindley  recursions  (2.11)  and  (2.17)  with  inputs 
characterized  by  a  common  pair  (A,  a)  and  release  rate  c  such  that  AE  [cr]  <  c  <  1. 
If  E[  cr2]  <  oo,  then  the  respective  stationary  rvs  q <*,  and  q(f)  satisfy 

E  boo]  =  E  [,«]  -  jE  H«t  -  1)] .  (4.48) 


In  the  case  where  the  multiplexer  release  rate  is  c  =  1/m,  for  some  integer 
m  =  1,2,...,  we  can  combine  (4.48)  with  (2.51).  Doing  so  yields  an  explicit 
expression  for  the  expected  queue  size  induced  by  M |  G |  oc  inputs,  namely 

AE  [a2]  (AmE  [a]  +  m  -  1) 


E  [goo.  = 


2(1  —  AmE  [a]) 


m  =  1,  2, . . . 


(4.49) 


provided  that  AE  [cr]  <  c  =  1/m  and  E  [cr2]  <  oo.  Of  course  formula  (4.49)  should 
be  consistent  with  the  light  traffic  derivatives  of  E  [g^],  calculated  in  Section  4.2.4 
through  the  Reiman-Simon  method.  Indeed,  for  m  =  1,  in  which  case  c  =  1,  it  is 
now  easy  to  differentiate  (4.49)  twice  and  verify  that  (4.25)  holds  true. 


4.3.4  Case  c  =  1:  Determination  of  P  [q^  =  0] 

Let  the  multiplexer  rate  be  c  =  1.  We  now  obtain  an  exact  expression  for 
P  [goo  =  0],  i.e.,  the  probability  that  the  queue  size  is  zero  in  the  discrete-time 
queue  with  M |  G |  oc  input  process.  This  is  accomplished  by  combining  the  corre¬ 
sponding  result  on  the  instantaneous  inputs  queue  with  a  sample  path  argument 
for  the  coupled  system  with  gradual  inputs.  The  details  are  as  follows: 
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Consider  the  systems  (2.17)  and  (4.33),  with  coupled  instantaneous  and  grad¬ 
ual  inputs  respectively.  In  Section  2.3.4  we  introduced  a  decomposition  of  the 
queue  length  evolution  in  system  (2.17)  in  terms  of  “idle”  and  “busy”  periods 
{(/<“’, 4”1),  n  =  l,2,...}  in  the  sense  defined  by  (2.40)  and  (2.41).  In  the  same 
manner  we  introduce  { ( I n ,  Bn),  n  =  1,2,.. .},  the  idle  and  busy  period  lengths  as¬ 
sociated  with  (4.33).  These  rv  pairs  are  also  i.i.d.;  let  (/,£?)  be  the  corresponding 
generic  pair.  We  first  note  that  the  coupled  systems  both  start  empty  and  have 
identical  cycle  lengths, 

In  +  Bn  =  I^  +  B^\  n  —  1,2, ... ,  (4.50) 

so  that  it  suffices  to  focus  on  one  such  common  regenerative  cycle,  say  the  first 
one. 

Next,  look  at  system  (2.17),  with  instantaneous  inputs.  Starting  from  q ^  =  0, 
the  queue  finds  itself  in  an  idle  period.  In  the  time  slot  preceding  the  first  busy 
period  there  must  be  either  a  single  session  arrival  whose  workload  exceeds  one, 
or  more  than  one  session  arrivals.  Now,  consider  the  coupled  system  (4.33),  with 
gradual  inputs.  Here,  on  every  sample  path,  the  queue  is  empty  whenever  it  is 
empty  in  the  sibling  instantaneous  input  system.  Clearly,  if  a  busy  period  in 
(2.17)  is  triggered  by  at  least  two  arrivals  in  the  preceding  slot,  then  the  queue  size 
in  (4.33)  also  grows  positive,  the  two  systems  become  simultaneously  busy,  and 
the  idle  period  I\  for  (4.33)  is  equal  to  the  respective  idle  period  /{u)  for  (2.17). 
However,  if  a  busy  period  in  (2.17)  is  initiated  by  a  single  session  arrival,  the  queue 
size  in  (4.33)  remains  zero,  as  long  as  only  one  session  is  active.  This  is  due  to  the 
assumption  that,  in  system  (4.33),  an  active  session  generates  input  at  rate  equal 
to  the  multiplexer  release  rate  c  =  1.  In  that  case  the  idle  period  I\  for  (4.33)  is 
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longer  than  the  respective  idle  period  /[u)  for  (2.17).  Let 

J:=h-  l[u)  (4.51) 

denote  their  difference.  Since  J  =  0  if  the  generic  busy  period  / ^  of  (2.17)  starts 
with  more  than  one  session  arrivals,  it  only  remains  to  condition  on  the  event  that 
it  starts  with  a  single  arrival. 

To  calculate  the  probability  that  the  busy  period  of  (2.17)  is  initiated  by  a  lone 
arriving  session  we  write 

P  [/(«>  =  k,  (3iM  =  1]  =  P  [c q lu)  =  •  •  •  =  =  0,  (3k  =  1,  <rM  >  1 

=  r)k-1P\J3  =  l]P[<j>i\,  k  —  1,2, ... , 
with  rj  given  by  (2.43),  so  that 

OO 

P  [&(«)=!]  =  ^  p[/(u)  =  fc>  &(»>  =  !] 

k= 1 

=  T— P  [8  =  1]  P  [(7  >  1] 

=  E  [/(u)j  P  \J3  =  1]  P  [(7  >  1]  ,  (4.52) 

where  the  last  equality  in  (4.52)  follows  from  (2.45).  Next,  note  that  on  the  event 
{/3j(u)  =  1}  the  queue  size  associated  with  (4.33)  becomes  positive  only  after  an 
arrival  of  a  second  session  which  takes  place  before  the  activity  period  of  the  first 
one  expires.  Let  the  rv  X  count  the  number  of  successive  time  slots  that  elapse 
until  a  second  session  arrives,  inclusive  of  the  slot  where  this  arrival  occurs.  The 
rv  X  is  geometrically  distributed  with 

P  [X  >  k\  =  P  \j3  =  Of  ,  k  =  0, 1, ... .  (4.53) 

From  the  discussion  above  it  follows  that  the  generic  rv  J  satisfies 

(  min([(7  —  l|cr  >  1],X)  if  /3lW  =  1 
J  =st  <  (4.54) 

0  if  (3j(u)  ^  1. 
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We  can  now  calculate  the  expectation  E  [J]  without  difficulty.  Using  (4.53)  we 


write 


E  [min ([<7  —  l|cr  >  1],  X)]  =  P  [min ([<7  —  1| o  >  1],  X)  >  k\ 

k= 0 
oo 

=  [a-  1  >  k\a  >  1]P  \J3  =  Of 


k= 0 


l 


P  [cr  >  1] 


P  [<r  >  k  +  1]  P  \/3  =  0]' 


k= 0 


P  [(7  >  1]  P  \J3  =  0] 


E  [min(<j,  X)  —  1]  . 


Therefore,  (4.54)  via  (4.52)  implies 


E  [  J]  =  p  || ~  (,j  E  [/<“>]  E  [min(a,  V)  -  1] . 


(4.55) 


By  application  of  the  Renewal-Reward  theorem,  in  conjuction  with  Lemma  4.3.1, 
we  determine  P  [q^  =  0]  as 

E  [I]  E  [/(“)]  +  E  [J] 


p  [qoc  =  o]  = 


(4.56) 


E  [B]  +  E  [I]  E  [£(“)]  +  E  [/(«)]  ’ 
where  the  second  equality  in  (4.56)  follows  from  (4.50)  and  (4.51).  Thus,  making 
use  of  (2.42)  and  (4.55)  in  (4.56)  above  we  collect 


Lemma  4.3.2  Consider  the  Lindley  recursions  (2.11)  and  (2.17)  with  inputs  char¬ 
acterized  by  a  common  pair  (A,  a),  with  generic  IN -valued  session  arrival  rv  /3  such 
that  X  =  E  \(3\  <  oo  and  release  rate  c  =  1.  If  AE  [cr]  <  1,  then 

p  [fc  =  0]  =  P  [<&)  =  0]  (l  +  p h  =  11  E  [min(<r,  X)  -  1])  ,  (4.57) 

where  the  rv  X  follows  the  geometric  distribution  (f.  53). 

We  stress  that  (4.57)  is  valid  for  any  LV-valued  sequence  of  i.i.d.  rvs  {/3n+i,  n  = 
0,1,.. .}  with  E  \/3\  <  oo,  as  the  arguments  leading  to  Lemma  4.3.2  do  not  require 
that  fi  be  Poisson  distributed. 
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By  specializing  (4.57)  to  M\G\oo  input  processes  and  invoking  (2.31)  we  obtain 
an  explicit  expression  for  P  [q^  =  0] ,  in  the  case  where  the  multiplexer  release  rate 
is  c  =  1. 


Proposition  4.3.3  Consider  the  Lindley  recursion  (2.11)  with  c  —  1.  Under  the 
stability  condition  AE  [a]  <  1,  it  holds  that 


P  [<?oo 


0]  =  (1  — AE  [a))ex 


1  +  a  e~Xkp  la  >  k\ 


k= 1 


(4.58) 


4.3.5  Case  c  =  1:  Short—  vs  long-range  dependence 


Let  the  multiplexer  release  rate  be  c  =  1.  In  this  section  we  seek  to  develop 
some  understanding  as  to  what  kind  of  light  traffic  results  should  be  expected 
when  Assumption  (C)  fails,  as  would  be  the  case  if  a  follows  some  subexponential 
distribution.  To  that  end  we  discuss  expansions  based  on  the  closed  form  expression 
(4.58)  for  P  [q^  =  0].  For  short-range  dependent  inputs  we  verify  that  (4.58)  is 
in  agreement  with  the  partial  light  traffic  information  for  P  [q^  =  0]  obtained  in 
Section  4.2.4.  More  interestingly,  we  show  that  (4.58)  can  be  exploited  to  provide 
a  light  traffic  limit  for  long-range  dependent  inputs;  such  a  result  could  not  have 
been  obtained  via  the  Reiman-Simon  theory  (at  least  in  its  present  form). 

We  start  by  considering  the  case  where  the  M\G\oo  input  process  is  short-range 
dependent.  From  (4.58)  it  follows  that 


P  [<Zoo  >  0]  =  1  —  (1  —  AE  [a])ex  1  +  A  ^(e“Afc  -  1)P  [a  >  k\  +  A(E  [a]  -  1) 

V  k= i  , 

(4.59) 

We  now  show  that  the  leading  term  in  the  light  traffic  expansion  of  P  [g^  >  0]  is 

of  order  A2  and  evaluate  this  term  explicitly.  To  do  this  note  that  the  mapping 

k  — »  —  (1  —  e~Xk)  monotonically  increases  to  the  mapping  k  — >  k  as  A  — »  0+. 
A 
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Therefore,  since  E  [a2]  <  oo  under  short-range  dependence,  we  conclude  that 

OO  OO  1 

.woi  -  e_Afc) p  w  >  k\  =  k  p  ^  =  2 E  _  x)]  (4-6°) 

k= 1  fc=l 

by  the  monotone  convergence  theorem,  while 

-  eA(l  -  AE  [cr])(l  -  A  +  AE  [a]))  =  E  [a]  (E  [a]  -  1)  +  (4. 

Zj 

Thus,  combining  (4.61)  and  (4.60)  with  (4.59)  yields 

Corollary  4.3.1  (Short-range  dependence)  In  the  setup  of  Proposition  4-3.3, 
with  E  [<r2]  <  oo,  it  holds  that 

lim  4  P  [9oo  >  0]  =  i  (E  [a(a  -  1)]  +  2E  [<r]  (E  [<r]  -  1)  +  1)  .  (4.62) 

A— >o+  2 

This  is  precisely  the  result  implied  by  (4.16)  and  (4.27)  for  b  =  0.  In  addition, 
we  observe  that  Assumption  (C)  is  superfluous  in  this  case,  for  what  is  needed  to 
obtain  (4.62)  is  that  E  [a2]  be  finite.  This  indicates  that  the  conclusions  of  Section 
4.2.4  may  still  be  valid  when  a  is  subexponentially  distributed  with  E  [cr2]  <  oo. 

Next,  we  turn  our  attention  to  the  situation  where  the  tail  of  a  is  regularly 
varying  with  index  —a  (1  <  a  <  2),  i.e.,  of  the  form 

P  [a  >  n\  =  n~ah(n),  n  =  1,2,...  (4.63) 

for  some  slowly  varying  function  h  :  H+  — >  H+.  In  this  case  E  [cr2]  =  oo,  the 
associated  M\G\oo  process  is  long-range  dependent  and  the  arguments  used  above 
do  not  apply.  First  write 

OO  -j  /  OO  \ 

e~XkP  [<T  >  k]  =  — — -  e~Xkp  W  =  k]~  e~x  (4.64) 

k= i  e  -  i  yk=i  J 

and  manipulate  (4.58)  to  obtain 

P[goo>0]  =  l-eA(l-AE[a])(^l  — A— ^Efajj 

(1  -  AE  [<t])  (  £  e-AtP  [<r  =  k]  -  1  +  AE  [cr])  .  (4.65) 


lim  —  (l 

a^o+A2  v 
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Expression  (4.65)  has  the  advantage  of  explicitly  displaying  the  Laplace-  Stieltjes 
transform  of  the  distribution  of  a.  This  puts  us  in  position  to  invoke  a  Tauberian 
result  on  the  asymptotic  behavior  of  Laplace-Stieltjes  transforms  at  the  origin.  In 
particular,  Theorem  8.1.6  in  [7,  p.  333]  provides  the  asymptotics  of  the  second 
term  in  (4.65)  as 

Y  e~XkV  [a  =  k]-l  +  AE  [a]  ~  Xah(l/X)  f  (A  -t  0+).  (4.66) 

fc=i  a 

For  the  first  term  in  (4.65),  making  use  of  the  fact 

e~x  —  1  T  A 

lim  - - - =  0,  1  <  a  <  2 

A— >0  A" 

we  find  after  straightforward  calculations  that 

&S  F  0  -  eA(1  -  AE  M)  (l  -  A  -  m)  )  =  0.  (4.67) 

Consequently,  (4.65)  via  (4.66)  and  (4.67)  leads  to 


Corollary  4.3.2  (Long-range  dependence)  In  the  setup  of  Proposition  4-3.3, 
with  the  tail  of  a  given  by  (f.63),  it  holds  that 


X~ 


lim  ,  .  . , . 
\^o+h(l/X) 


P  [^oo  >  0] 


T(2-q) 

a  —  1 


(4.68) 


The  effect  of  the  distribution  of  a,  and  in  particular  of  its  second  moment, 
on  the  light  traffic  asymptotics  of  P  [goo  >  0]  is  now  more  apparent:  As  A  — *  0+, 
the  “busy”  queue  probability  P  [q^  >  0]  exhibits  a  A2  decay  under  short-range 
dependence,  and  the  slower  non- integer  power  decay  A",  when  a  satisfies  (4.63),  in 
which  case  the  M\G\oo  process  is  long-range  dependent.  The  limit  (4.68)  prompts 
us  to  conjecture  that,  when  c  =  1,  under  long-range  dependence,  such  a  A"  decay 
in  A  of  the  tail  probability  P  [q^  >  b]  holds  more  generally  for  all  b  >  0.  Yet  this 
problem  remains  as  of  now  unresolved. 
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Figure  4.4:  Calculation  of  -0(Ay,  ^2);  example  for  Ay  —  2,  Ay  —  5. 

4.4  Proof  of  Lemma  4.2.1 


We  consider  parts  (a),  (b)  and  (c)  separately. 

Proof  of  (a).  By  symmetry  it  suffices  to  determine  i/>(Ay,  Ay )  for  k\  <  Ay.  To  do 
this,  we  refer  to  Figure  4.3,  and  split  the  area  where  Ay,  Ay})  >  0  in  ten 

numbered  regions,  adopting  the  convention  that  edges  between  region  i  and  any 
other  higher  numbered  region  all  belong  to  region  i.  The  splitting  we  choose  is 
shown  in  Figure  4.4,  accompanied  by  a  simple  example  for  the  case  Ay  =  2,  Ay  =  5. 

The  double  sum  (4.19)  is  broken  up  into  sums  over  each  region,  which  are 
calculated  below. 


Region  1: 


Region  2: 


— k2  — fc2 

E  E  +  b 2 )  = 

tl=~ (fel+fe2)  t2=tl 
1  — fc2 


Ay(Ay  +  1)  (A^i  +  2) 


E  (Ai  —  (1  —  Ay)  +  Ay)(Ay  —  k 


i  = 


ki{k2 


ti=i— (fci+fc2) 

— fci  *i+fci 

Region  3:  E  E  <*.-*+*) 

2/ci+l  t2=— 


ki(k\  +  1)  (Aii  —  1) 
6 


&i)(^i  +  1) 

~2 


o  o 

Region  4:  E  E(-« 

tl=l  —  ^1  t2=t\ 


Aii  (Aii  +  l)(^i  —  1) 
6 
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ki(ki  —  1)  {k\  +  1) 


h(k2  -  h)(h  +  1) 


ki(ki  —  1)  (k\  +  1) 


o  h-i 

Region  5:  E  E  (-n) 

ti=l—ki  t2=—k\ 

0  -fci-1 

Region  6:  (-*i)  = 

t\=—k\  t2=—k2 
1  -k2-l 

Region  7:  E  E  {t2  —  t\  +  k2 ) 

tl=— fcl  t2=tl—k2 
-k2 

Region  8:  -/ci  +  (t2  +  +  fc2)(/c2  -  fci) 

*2  =  —  (^1+^2) 

— k2  ti  — 1  — k2  —  1 

Region  9:  E  E  (t2  +  +  k2 )  —  ^  )  (t2  +  +  A;2) 

tl=—  (fcl+fc2)+l  t2  =  —  (fcl+fe)  t2  =  —  (kl+k2) 


6 


&i((^i  +  1)(^2  —  ^1)  —  2) 


—  /Cl  ti—  1 

Region  10:  E  E  h  —  (2  k2  —  2k\  —  l)ki 

t\=—k2+l  t2=—k2 


-  l)(fci  -  2) 
6 


k\{k2  —  k\  —  1  ){k2  -  k\  -  2) 


Adding  the  results  from  regions  1-10  shows  that  the  claim  holds  true. 


Proof  of  (b).  It  suffices  to  evaluate  i/j  for  ki  <  k2.  We  refer  to  the  proof  of 
Part  (a),  split  the  t\t2  plane  in  the  same  manner  and  carry  out  the  algebra  for  the 
double  sum  (4.19)  corresponding  to  ^  :=  g„.  The  summands  are  listed  below. 

v— ;  v-^ ,  s2  fci(fci  +  l)2(/ci  +  2) 

Region  1:  zJ  ^  +  k2)~  = - 

ti=—{ki+k2)  t2=ti 


12 


1 — k2 


Region  2:  (k2  -  ki)  ^  (ii  -  (1  -  k2)  +  k2)‘ 

tx=l—(k\+k2) 


k\{k\  +  l)(2fc1  +  1  ){k2  —  fci) 
6 
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—fcl  tl+fcl 

Region  3:  ^  ^  (ti  -  t2  +  hf 

t\= — 2fci+l  t2— — fciH-1 
0  0 


*?(*?  - 1) 


12 


Region  4:  E  E<-*>! 

tl=l  —  fcl  t2=t\ 

0  tl-1 

Region  5:  E  E(-^ 


t\=l—ki  t2=—ki 
0  — fci  — 1 


Region  6:  E  E(-M 


tl=— fcl  t2  =  —  k2 
1  -k2- 1 


-  1) 

12 

~  1) 

12 

k\{k\  +  1)(2A;1  +  1)(&2  — 

'  6 


Region  7:  E  E  (#2  —  ^1  +  ^2^ 


*f(*f  *  1) 


ti=—ki  t2=ti—k2 


12 


~k2 


Region  8:  (k2  -  k±)  ^  (t2  +  fa  +  k2)2  -  k\ 

t2=—{ki+k2) 

fa({fa  +  l)[2ki  +  1  )(k2  -  ki)  -  6fa) 


6 


— &2  ti  — 1  —  k2  —  1 

\2 


Region  9:  E  E  [1 2  +  fa  +  k2)2  —  ^  ^  (^2  +  +  k2)~ 

t]_=  —  (fci+fc2)  +  l  ^2  —  —  (fcl_l_fc2)  ^2  = —  (fcl_l_fc2) 

/c1(/c1-1)2(A:1-2) 


12 


—  fcl  ti—  1 

Region  10:  E  E  kf  —  (2  k2  —  2/ci  —  l)/c2 

tl  =— A.2+1  t2  =  ~k2 


k\{k\  —  k2  +  l)(ki  —  &2  +  2) 


Adding  the  summands  from  regions  1-10  verifies  the  asserted  expression. 


Proof  of  (c).  By  symmetry  we  need  only  calculate  ip  for  fa  <  k2.  We  go  back  to 
Figure  4.3,  where  ip  corresponding  to  ip  :=  q0  was  displayed.  Clearly,  in  order  for 
q0  to  be  greater  than  b  it  is  necessary  that  fa  >  b.  If  this  is  not  true,  i.e.  if  fa  <  b, 
then  1  [g0  >  b]  is  zero,  hence  ip  corresponding  to  ip  :=  1  [go  >  b\  is  also  zero  and 
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-(k,+k2)  b+l-^+kp  -k(  -b-1  0 

Time  t, 


Figure  4.5:  Values  of  1  [g0  >  b]  ({t1; t2 ;  ki,  k2})  for  ki  <  k2. 

the  second  part  of  the  claim  holds  true.  Now,  if  fci  >  6  then  the  tit2-plane  does 

indeed  contain  a  region  where  ip  —  1.  We  show  ip({ti,  t2\  ki,  k2})  corresponding  to 

ip  :  =  1  [q  >  b]  and  identify  this  region  in  Figure  4.5.  Using  this  graph,  calculation 

of  ip(ki,k2)  from  the  double  sum  (4.6)  is  a  matter  of  algebra.  From  the  points 

along  the  line  t\  =  —ki,  which  splits  the  ip  =  1  region  in  two  trapezoids,  we  collect 

-6-1 

Y  l  =  k1  +  k2-2b-l.  (4.69) 

i=-(fcl+fc2)+b+l 

On  the  right  trapezoid  the  double  sum  can  be  calculated  as 

k2—b—2  ^ 

(j  +  ki  -  b)  =  ~(2ki  +  k2  -  3b  -  2)(fc2  -  b  -  1).  (4.70) 

3=0 

On  the  left  trapezoid  we  find 

ki—b—2  ^ 

Y  {k1  +  k2-2b-2-  j )  =  i(2  k2  +  kx  -  3b  -  2)(h  -  b  -  1).  (4.71) 

3=0 

Adding  (4.69),  (4.70),  (4.71)  and  manipulating  validates  the  claim  for  ki  >  b.  ■ 
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Chapter  5 


Interpolation  approximations 

5.1  Introduction 

In  this  chapter  we  address  the  problem  of  evaluating  the  buffer  content  distribution 
at  a  multiplexer  fed  by  an  M\G\oo  arrival  process,  across  the  entire  range  of  sys¬ 
tem  utilizations  between  zero  and  one.  For  arbitrary  session  duration  distribution 
G  the  system  lacks  the  desired  Markovian  structure  and  obtaining  exact  analyti¬ 
cal  results  appears  very  difficult,  if  not  impossible.  To  circumvent  the  difficulties 
of  an  exact  analysis  one  typically  seeks  to  devise  approximations,  by  relying  on 
information  gleaned  from  asymptotic  regimes.  A  promising  approach  consists  of 
deriving  approximations  from  the  analysis  of  large  buffer  asymptotics  [29,  39,  47]; 
these  estimates  are  exact  in  the  limit  as  the  buffer  level  goes  to  infinity.  Here,  we 
propose  an  alternative  class  of  simple  approximations,  justified  by  the  character¬ 
izations  of  the  buffer  content  distribution  in  the  heavy  and  light  traffic  regimes, 
obtained  in  Chapters  3  and  4,  respectively.  These  approximations  are  termed  inter¬ 
polation  approximations  [49],  because  they  arise  by  suitably  interpolating  between 
the  heavy  and  light  traffic  limits  of  the  quantity  of  interest.  Such  approximations 
are  asymptotically  exact  in  the  limits  as  the  system  utilization  goes  to  zero  and 
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one.  We  provide  examples  of  interpolation  approximations  to  the  buffer  content 
distribution,  for  several  commonly  chosen  G,  covering  both  short-  and  long-range 
dependent  M\G\oo  inputs.  Comparisons  with  simulation  estimates  suggest  that 
the  approximants  capture  accurately  the  queue  size  distribution  at  small  buffer 
sizes,  for  which  approximations  based  on  large  buffer  asymptotics  are  often  ill  fit¬ 
ted.  On  the  other  hand,  when  G  has  finite  exponential  moment,  we  do  not  expect 
the  heavy-light  traffic  interpolation  to  be  accurate  for  buffer  sizes  much  larger 
than  the  maximum  burst  length:  The  approximation  exhibits  the  appropriate  ex¬ 
ponential  decay  in  the  buffer  size,  yet  the  rate  is  only  asymptotically  exact  as  the 
system  utilization  tends  to  one,  i.e.,  in  the  heavy  traffic  limit.  This  drawback  is 
often  absent  under  long-range  dependence,  as  there  are  cases  where  the  queue  size 
distribution  has  hyperbolic  (in  the  buffer  size)  asymptotics  with  the  same  expo¬ 
nent  for  all  traffic  intensities!  Moreover,  an  approximation  is  more  valuable  in  the 
presence  of  heavy  tails,  when  considering  that  alternative  estimates  by  means  of 
simulation  take  an  unreasonably  long  time  to  obtain. 

5.2  Summary  of  asymptotics 

In  the  context  of  the  Lindley  recursion  (2.11),  describing  the  evolution  of  the  buffer 
content  at  a  multiplexer  with  M\G\oo  inputs,  we  are  interested  in  approximating 
the  probability  that  the  stationary  buffer  content  exceeds  b  when  the  system  uti¬ 
lization  is  p  :=  AE  [cr]  /c, 

P(b,  p)  :=  PA  [goo  >  b]  ,  b>  0,  0  <  p  <  1.  (5.1) 

Let  us  point  out  that  the  system  dynamics  depend  on  A,  G  and  c  jointly,  and  not 
simply  on  the  utilization  p.  However,  for  developing  approximate  expressions  it  is 
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convenient  to  fix  G  and  c  and  adopt  the  view  suggested  by  (5.1),  that  is  consider 
the  buffer  content  distribution  as  a  function  of  the  system  utilization  p.  In  a  similar 
manner  we  introduce  the  moments  of  the  stationary  buffer  content 

Qk(p)  :=  Ea  [qook]  ,  0<p<l,  k  =  1,2,...; 

these  should  not  be  confused  with  the  heavy  traffic  queue  length  process  {Q(t),  t  > 
0}  of  Chapter  3. 

The  interpolation  approximations  we  have  in  mind  hinge  on  the  availability  of 
explicit  expressions  for  limits  of  system  quantities  as  A  — >  c/E  [a]  (heavy  traffic 
limits),  and  derivatives  with  respect  to  A  as  A  — >■  0  (light  traffic  derivatives). 
For  notational  convenience  we  now  rephrase  the  required  light  and  heavy  traffic 
asymptotics  in  terms  of  the  utilization  p,  with  the  understanding  that  when  p  — >  0 
or  p  — >  1  it  is  actually  A  which  goes  to  the  corresponding  limit,  while  both  c  and 
cr  remain  fixed. 

We  start  from  the  results  in  the  light  traffic  regime,  obtained  by  the  Reiman- 
Simon  method: 


Proposition  5.2.1  Consider  the  Lindley  recursion  (4-1)  with  release  rate  c  >  1 
and  let  b  >  0.  If  Assumption  (C)  is  satisfied,  then  the  following  hold: 

(a)  For  each  n  —  0, 1, . . . ,  [c\,  we  have 


dn 

dp1 


-P(b,  0+)  —  0; 


(5,2) 


(b)  In  addition,  for  c  =  1 ,  we  have 


d2 


E[a]2—p(b,0+)  =  E  \(a-b)+~\  P  [a  >  b]  +  2  E  [(a  -  b)+] 


-3  E  [(a  -  b)+]  P  [a  >  b]  +  P  [a  >  bf  . 


(5.3) 
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Considered  next  is  the  behavior  of  the  queue  with  M\G |oo  arrivals  in  heavy 
traffic,  that  is,  as  the  packet  arrival  rate  AE  [cr]  tends  to  the  multiplexer  release 
rate  c  from  below.  We  tacitly  assume  that  the  heavy  traffic  limit  of  the  stationary 
distribution  coincides  with  the  stationary  distribution  of  the  heavy  traffic  limit. 
Using  Assumption  (B)  we  express  the  relevant  facts  from  Theorems  3.4.1  and 
3.4.3  in  terms  of  the  utilization  p  in  the  following 

Proposition  5.2.2  The  heavy  traffic  limits  of  the  stationary  queue  length  distri¬ 
bution  associated  with  (2.11)  can  be  classified  as  follows: 

(a)  If  E  [a2]  <  oo,  then 

hmP a  [(1  -  p)  <?oo  >  x\  =  exp  ( -  ~E  xV  x  >  0.  (5.4) 

p-> i  \  E  [az \  J 

(b)  If  P  [a  >  n\  =  n~a ,  n  =  1,2,...,  with  1  <  a  <  2,  then 

limP a  [(1  -  d)1/(Q_1)<7oo  >  x\=  Ea_ i  'P  xa~1')  ,  x  >  0,  (5.5) 

p->i  \  1  (2  —  a)  J 

where  for  v  >  0,  Ev  :  R  — >  R  is  the  Mittag-Leffler  special  function  defined  by 
(3.19). 

Part  (a)  of  Proposition  5.2.2  addresses  the  classical  short-range  dependent 
case,  for  which  the  heavy  traffic  normalizer  is  1  —  p  and  the  limiting  heavy  traffic 
distribution  is  exponential.  Part  (b)  deals  with  a  long-range  dependent  M\G\oo 
arrival  process,  in  which  case  the  heavy  traffic  queue  length  distribution  is  ex¬ 
pressed  through  a  Mittag-Leffler  function  with  hyperbolic  decay,  while  the  heavy 
traffic  normalizer  is  (1  —  p)1/^"-1)  ancl  has  power-law  behavior. 

The  results  under  the  fight  and  heavy  traffic  regimes  are  subsequently  combined 
into  approximations  for  all  values  of  p  in  the  interval  [0, 1). 


103 


5.3  Heavy— light  traffic  interpolations 


Whenever  Assumption  (C)  is  satisfied,  Pa  [Qoo  >  b\  is  infinitely  differentiable  with 
respect  to  p  at  p  =  0,  hence  it  can  be  approximated  by  bringing  together  heavy 
traffic  limits  and  light  traffic  derivatives  into  a  Taylor  series-like  expansion.  To  this 
end  we  enforce  Assumption  (C)  throughout  Sections  5.3.1  and  5.3.2  and  follow  the 
approach  proposed  in  [21].  In  passing,  we  also  discuss  approximations  for  Qk(p), 
k  —  1,2.  The  details  are  given  below: 

5.3.1  Tail  probability  approximations 

Consider  the  normalized  queue  length  rv  (1  —  p)  and  define 


F(x,p)  :=  PA  [(1  -  p)  <?<x>  >  x] 
x 


=  P 


1  ~P 


P 


0  <  p  <  1,  £  >  0 


and 


F(x,  1)  :=  limPA  [(1  -  p)  Qoo  >  x\  ,  x  >  0. 

p-S>l 


(5.6) 


(5.7) 


Assume  that  partial  derivatives  of  F(x,  p)  with  respect  to  p,  up  to  order  n,  at  p  — 
0T,  are  available.  Construct  Fn(x,p),  the  nth  order  interpolation  approximation 
to  F(x,  p),  by  means  of  the  polynomial 


Fn(x,  F(x,  0+)  +  (  F(x,  1)  -  ^  ^  F{x,  0+) 


1 


i= 0 


i\  dp 


i=0 


i\  dp1 


P 


n+1 


(5.8) 


Observe  that 


^  dl  ^  dl 

Fn{x,  1)  =  F(x,  1)  and  —  Fn(x,  0+)  =  —  F(x,  0+),  i  =  0,l,...,n, 

that  is,  Fn(x,  p)  is  precisely  that  unique  n+1  degree  polynomial  in  p  which  matches 
the  n  +  1  partial  derivatives  of  F(x,p)  at  p  =  0+  and  its  heavy  traffic  limit. 
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Now,  by  reversing  the  (1  —  p)  normalization  in  Fn(x,p )  we  generate  the  nth  order 
interpolation  approximation  to  Pa  [goo  >  b]  as 

Pa  [goo  >  b]  w  Fn((  1  -  p)  b ,  p).  (5.9) 

Note  that,  in  principle,  this  may  lie  outside  [0, 1],  in  which  case  it  is  obviously  a 
poor  approximation. 

To  calculate  the  quantities  associated  with  (5.9)  it  remains  to  express  the  partial 
derivatives  appearing  in  (5.8)  in  terms  of  the  light  traffic  derivatives  of  Pa  [goo  >b]. 
We  have 

F(x,  0+)  =  P(x,  0+)  +  x  J^P(x,  0+)  (5.10) 

and 

i?F(x' 0+)  =  b-p(x'0+)+2xid~F(x'0+) 

d  d 2 

+2x  — P(x ,  0+)  +  x2  — -P(x,  0+).  (5-11) 

ax  ax2 

In  case  additional  light  traffic  information  is  available,  repeated  application  of  the 
chain  rule  will  yield  higher  order  derivatives,  as  needed. 

We  are  now  ready  to  write  approximate  expressions  anchored  on  the  heavy  and 
light  traffic  results  of  Section  5.2.  Proposition  5.2.2(a)  provides  the  limit  (5.7)  that 
should  be  inserted  in  (5.8).  Proposition  5.2.1(a)  can  be  used  to  substitute  for  the 
partials  in  (5.10)  and  then  in  (5.8).  Thus,  if  the  multiplexer  release  rate  is  c  >  1, 
the  \c\th  order  interpolation  approximation  to  Pa  [goo  >  b]  is  simply 

Pa  [goo  >  b]  «  FLcJ  ((1  -  p)  b ,  p)  =  pLcJ+1  exp  i1  ~  p)  b^j  •  (5-12) 

More  can  be  accomplished  in  the  case  c  =  1,  since  Proposition  5.2.1(b)  affords  us  a 
promising  2nd  order  interpolation  approximation.  Insertion  of  (5.11)  in  (5.8)  yields 

F2 (b,  p)  =  ip2(l  -  p)  P(b,  0+)  +  p 3  exp  b^j  (5-13) 
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and  the  latter  leads  to  the  2nd  order  approximation 


Pa  [qoo  >  b]  «  F2((l  -  p)  b,p),  c=l,  (5.14) 

where  Proposition  5.2.1(b)  is  used  to  supply  the  second  partial  derivative  in  (5.13). 


5.3.2  Moment  approximations 


Next,  we  briefly  deal  with  moment  approximations.  We  restrict  attention  to  the 
case  c  =  1  and  consider  only  the  queue  length  first  and  second  moment.  The 
relevant  light  traffic  limits  are  given  by  (4.15),  (4.25)  and  (4.26).  In  heavy  traffic 
we  see  from  (5.4)  that 

lim(l  -  p)k  Qk(p)  =  k\  >  k  =  1,2, . . . . 


Moment  approximations  are  then  developed  by  interpolating  for  (1  —  p)Qi(p)  and 
(1  —  p)2  Qoip)-  in  very  much  the  same  manner  as  distribution  approximations.  We 
skip  the  details  of  the  derivation  and  list  the  resulting  final  expressions 

E  k2l  p2 


Qi(p)  = 


2E  [a]  1  -p 


c  =  1 


(5,15) 


and 


Q^ip) 


P ' 


4(1  -  Py 


E  \a 


2)2 


EM 


-1  + 


E\a 


21 2 


EH 


+  1 


c  =  1. 


(5.16) 


Note  that  formula  (5.15)  is  in  fact  exact,  as  it  coincides  with  result  (4.49)  (with 
m  —  1)  and  with  a  continuous  time  analog  established  for  a  fluid  model  in  [55,  p. 
23].  This  match  clearly  validates  the  interpolation  method.  On  the  contrary  we 
note  that  (5.16)  cannot  be  exact.  To  see  this,  consider  the  example  where  a  —  l 
deterministic.  This  corresponds  to  i.i.d.  Poisson  arrivals,  for  which  a  probability 
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generating  function  of  the  queue  length  rv  is  available.  Using  (2.32)  it  can  be 
shown  that  the  exact  expression  is 

p 2 

Q2(p)  =  6(1  _  (p2  -  p  +  3),  c=l,  <7  =  1  a.s. 

a  formula  that  clearly  cannot  be  recovered  using  only  two  light  traffic  derivatives. 
Still,  when  <7  =  1  approximation  (5.16)  is  within  9%  of  the  correct  value,  for  all  p 
in  [0,1). 

5.3.3  Long-range  dependence 

It  is  apparent  from  the  developments  of  Section  5.2  that  the  light  traffic  results,  as 
stated  in  Proposition  5.2.1,  do  not  cover  several  interesting  distributions  belonging 
to  the  subexponential  family.  Such  is  for  example  the  lognormal  distribution,  which 
violates  Assumption  (C)  despite  having  finite  kth  moment  for  all  k  —  0, 1, . . ..  In 
view  of  Corollary  4.3.1,  it  is  natural  to  expect  that  Assumption  (C)  can  be  relaxed 
to  require  that  E  [crfc]  be  finite,  for  appropriate  k  >  2,  in  order  for  Proposition 
(5.2.1)  to  go  through.  This  would  still  not  address  the  case  of  long-range  depen¬ 
dence,  characterized  by  E  [cr2]  =  oo.  We  are  however  able  to  construct  a  sharp 
approximation,  based  on  heuristic  arguments  presented  below. 

Consider  the  Pareto  distribution  P  [cr  >  n\  =  n~a,  n  =  1,2,...,  with  1  <  a  <  2; 
in  this  case  not  only  Assumption  (C)  fails,  but  as  seen  from  Proposition  5.2.1,  (5.3) 
yields  infinity.  This  indicates  that  Pa  [qoo  >  b]  may  not  be  an  analytic  function  of  p 
under  long-range  dependence.  When  c  =  1  we  recall  that  Corollary  4.3.2  strongly 
suggests  that  limp~"PA  [<7oo  >  b]  is  the  sought  after  non-trivial  limit,  for  all  b  >  0. 

p— >0 

These  considerations  lead  us  to  postulate  that,  when  c  =  1, 

limp-“PA  [Qoo  >b}  =  K(b)  (5.17) 
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for  some  unknown  mapping  K  :  R+  — y  R+,  which,  by  Corollary  4.3.2,  satishes 
K( 0)  =  T(2  —  a) /{a.  —  1)E  [cr]Q.  On  the  other  hand,  the  heavy  traffic  result  of 
Proposition  5.2.2(b)  hints  at  developing  an  approximation  around  the  normalized 
rv  (1  —  p)1/!"- b  q^.  Then,  taking  advantage  of  Proposition  5.2.2(b)  we  propose 
the  approximant 


P>  [fa  >  6|  «  [  -(°„7.1)E!g|  ^  6”- 


c—1. 


(5,18) 


T(2  —  a )  p° 

This  expression  is  in  agreement  with  the  heavy  traffic  limit  (5.5).  In  addition,  from 
the  Mittag-Leffler  function  asymptotics  given  in  [17,  p.  207],  we  have 


Ea_i(-x) 


1  1 


(X  — »  oo 


x  T(2  —  a) 

which  ensures  that,  as  p  — >  0,  approximation  (5.18)  conforms  with  the  conjectured 
fight  traffic  limit  (5.17). 


We  close  the  presentation  of  the  approximate  expresssions  with  a  comment. 
Recall  that  each  active  source  in  the  M\G\oo  arrival  process  generates  one  infor¬ 
mation  unit  per  time  slot.  So,  c  =  1  corresponds  to  the  case  where  the  amount  of 
service  in  one  slot  is  exactly  equal  to  the  amount  of  information  that  one  active 
source  generates  in  one  slot.  When  c  =  la  single  active  source  suffices  to  make 
full  use  of  the  server  capacity;  in  this  system  there  is  never  any  leftover  capacity 
to  simultaneously  serve  more  than  one  source.  On  the  contrary,  when  c  >  1,  the 
server  can  attend  to  more  than  one  source  during  one  time  slot,  so  that  there  is  a 
multiple  service  feature  to  the  system  behavior.  An  exact  or  approximate  analysis 
in  this  regime  is  clearly  more  challenging. 
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5.4  Numerical  results 


To  gauge  the  accuracy  of  the  proposed  expressions  we  have  carried  out  simulation 
experiments  under  various  choices  for  the  distribution  of  the  session  duration  rv 
cr.  The  experimental  values  are  obtained  by  regenerative  simulation  and  relative 
widths  accompanying  them  correspond  to  95%  confidence  intervals.  We  almost 
exclusively  (with  one  exception)  confine  ourselves  to  the  simple  situation  where  the 
multiplexer  release  rate  is  c  =  1.  While  the  list  of  examples  below  is  not  exhaustive, 
it  does  serve  to  illustrate  the  ability  of  the  heavy-light  traffic  interpolation  to 
“ballpark”  the  true  tail  probabilities,  as  well  as  its  limitations. 

Deterministic  When  the  session  duration  is  deterministic,  a  =  D  a.s.,  for  some 
positive  integer  D ,  approximation  (5.14)  reads 

Pa  [fe>6]  «  dTld  (3[£,  _  (l  _  /&]+([£>  _  (l  _  rt6]+  _  l) 

+  1  \D  >  (1  -  p)b\)  +  /ex p  (-4(1  -  />)'>)  ,  6=0,1,.... 

We  let  the  session  duration  he  cr  =  3  and  obtain  simulation  estimates  for  the 
steady  state  probability  Pa  [qoo  >0].  Of  course,  in  this  case  the  exact  expression 
(4.58)  is  available.  In  Table  5.1  we  list  simulation  estimates  and  numerical  values 
from  (4.58)  and  from  the  light-heavy  traffic  interpolation.  A  comparison  of  the 
exact  formula  (4.58)  to  the  light-heavy  traffic  interpolation  shows  that,  in  this  case, 
the  agreement  is  excellent.  Since  we  expect  the  approximation  to  be  asymptotically 
exact  at  the  endpoints  p  =  0  and  p  =  1,  it  is  not  surprising  that  the  largest  errors 
occur  in  moderate  traffic. 

In  the  same  setup,  we  next  consider  the  tail  probability  Pa  [qoo  >4].  No  exact 
expressions  are  available  in  this  case.  From  Table  5.2  we  see  that  although  the 
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Tail  probability  Pa  [q 

oo  >  0] 

p 

Exact 

Simulation 

Approximation 

Error  (%) 

0.1 

1.0478e-02 

1.0469e-02±0.2% 

1.0500e-02 

-0.21 

0.2 

4.1622e-02 

4.1668e-02±0.3% 

4.1778e-02 

-0.38 

0.3 

9.3042e-02 

9.3020e-02±0.2% 

9.3500e-02 

-0.49 

0.4 

1.6441e-01 

1.6442e-01±0.2% 

1.6533e-01 

-0.56 

0.5 

2.5545e-01 

2.5533e-01±0.2% 

2.5694e-01 

-0.58 

0.6 

3.6594e-01 

3.6607e-01±0.1% 

3.6800e-01 

-0.56 

0.7 

4.9573e-01 

4.9601e-01±0.1% 

4.9817e-01 

-0.49 

0.8 

6.4470e-01 

6.4488e-01±0.1% 

6.4711e-01 

-0.37 

0.9 

8.1279e-01 

8.1264e-01±0.1% 

8.1450e-01 

-0.21 

Table  5.1:  Pa  [goo  >  0]  for  deterministic  session  duration  a  =  3. 

approximation  yields  estimates  in  the  correct  order  of  magnitude,  the  errors  are 
substantial  when  not  in  the  moderate-to-heavy  traffic  regime.  This  can  be  ex¬ 
plained  as  follows:  When  a  =  3,  in  order  for  the  queue  to  build  up  to  4  at  least 
3  sources  should  be  simultaneously  active.  Note  that  the  light  traffic  component 
of  the  approximation  consists  of  the  second  derivative,  which  can  be  obtained  by 
considering  sample  paths  with  at  most  two  source  activations  in  the  system.  Thus, 
any  effects  due  to  the  activation  of  more  than  two  sources  are  not  adequately 
accounted  for  in  fight  traffic. 

Uniform  We  now  specialize  (5.14)  to  the  case  where  a  is  uniformly  distributed 
on  {1, . . . ,  M},  i.e.,  P  [a  =  n]  =  1/M,  n  —  1,  2, . . . ,  M.  This  yields 

Pa  [loo  >  6]  «  1  [M  >  (1  -  p)i>]  3jffM~+rt1)2  (1  +  5(M  -  (1  -  p) bf) 
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Tail  probability  Pa  [goo  >  4] 

p 

Simulation 

Approximation 

Error  (%) 

0.1 

1.1271e-04±1.7% 

9.0718e-05 

19.51 

0.2 

1.3444e-03±1.7% 

9.4753e-04 

29.52 

0.3 

6.1736e-03±0.9% 

5.9952e-03 

2.89 

0.4 

1.9246e-02±0.6% 

1.4415e-02 

25.10 

0.5 

4.7745e-02±0.5% 

3.9894e-02 

16.44 

0.6 

1.0292e-01±0.4% 

9.5777e-02 

6.94 

0.7 

2.0035e-01±0.3% 

1.9757e-01 

1.39 

0.8 

3.6093e-01±0.3% 

3.6379e-01 

-0.79 

0.9 

6.1407e-01±0.3% 

6.1902e-01 

-0.81 

Table  5.2:  Pa  [goo  >  4]  for  deterministic  session  duration  a  =  3. 

x  (M  —  (1  —  p)b)2  +  p3  exp  ^-6  2M+^l)  ’  ^  =  0, 1, . . . . 

For  M  =  5  we  compare  simulation  vs  approximation  in  Tables  5.3  and  5.4,  for 
system  utilizations  p  =  0.2  and  p  =  0.8  respectively.  Once  more,  the  approxima¬ 
tion  is  very  sharp  for  small  buffer  sizes.  As  the  buffer  size  increases  beyond  the 
maximum  burst  length  and  the  true  probabilities  become  smaller,  the  approxima¬ 
tion  lingers  on  in  the  correct  order  of  magnitude,  but  clearly  deteriorates  away 
from  heavy  traffic.  Eventually,  as  the  buffer  size  tends  to  infinity,  the  interpolation 
approximation  overestimates  the  actual  probabilities. 
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Tail  probability  Pa  [goo  >  b] 

Buffer  size 

Simulation 

Approximation 

Error  (%) 

0 

4.4907e-02±0.2% 

4.5333e-02 

-0.95 

2 

1.1747e-02±0.5% 

1.1399e-02 

2.97 

5 

1.4543e-03±1.4% 

9.7380e-04 

33.04 

8 

1.9456e-04±3.7% 

2.4379e-04 

-25.30 

10 

5.1620e-05±7.0% 

1.0186e-04 

-97.31 

Table  5.3:  Utilization  p  =  0.2;  a  ~  uniform(l,  5). 


Tail  probability  Pa  [^oo  >  b] 

Buffer  size 

Simulation 

Approximation 

Error  (%) 

0 

6.5489e-01±0.1% 

6.6133e-01 

-0.98 

10 

2.0086e-01±0.4% 

1.9161e-01 

4.60 

20 

6.3303e-02±0.9% 

5.8057e-02 

8.29 

30 

1.9964e-02±1.7% 

1.9406e-02 

2.79 

40 

6.2827e-03±3.1% 

6.5188e-03 

-3.75 

50 

1.9703e-03±5.6% 

2.1897e-03 

-11.13 

Table  5.4:  Utilization  p  =  0.8;  a  ~  nniform(l,  5). 
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Figure  5.1:  Geometric  7  =  0.8  session  duration. 

Geometric  Taking  a  to  follow  the  geometric  distribution.  P  [er  >  n]  =  7",  n  = 
0,1,.. with  0  <  7  <  1  we  obtain  from  (5.14)  that 

Pa  [goo  >  b]  «  ^(l-p)(l+7)272(1'p)6  +  P3exp  (-2  ^(1  -  p)bj  ,  b  =  0, 1, . . . 

As  an  example  we  set  7  =  0.8  and  plot  simulated  and  approximate  values  in 
Figure  5.1,  for  system  utilizations  p  =  0.1,  0.4,  0.6  and  0.9.  In  all  cases  confidence 
interval  widths  were  within  10%  of  the  mean.  The  linear  decrease  of  the  simulated 
values  suggests  an  exponential  decay  of  the  queue  length  distribution,  in  agreement 
with  large  deviations  results.  Figure  5.1  clearly  indicates  that  the  heavy-light 
traffic  interpolation  is  sufficient  for  providing  rough  estimates  for  a  wide  range  of 
probabilities  and  buffer  sizes. 

An  example  for  c  >  1  If  the  multiplexer  release  rate  is  c  >  1,  expression  (5.14) 
is  not  applicable,  the  first  non-zero  light  traffic  derivative  is  not  available,  and  the 
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Figure  5.2:  Geometric  7  =  2/3  session  duration;  release  rate  c  =  4. 

only  available  heavy-light  traffic  approximant  is  (5.12).  To  illustrate  its  behavior 
we  provide  an  example  in  Figure  5.2,  where  we  have  picked  c  =  4  and  a  geometric 
with  parameter  7  =  2/3.  The  results  correspond  to  system  utilizations  p  =  0.2, 
p  =  0.6  and  p  =  0.8.  The  approximation  fares  very  well  in  moderate  to  heavy 
loads,  yet  it  obviously  yields  inaccurate  results  for  p  =  0.2  due  to  insufficient  light 
traffic  information. 


Pareto  Let  the  session  duration  rv  follow  the  Pareto  distribution  P  [cr  >  n\  = 
n~a,  n  =  1,2,...,  with  1  <  a  <  2,  which  case  the  M\G\oo  process  is  long- 
range  dependent  and  the  approximate  expression  (5.18)  is  in  effect.  Assessing  the 
performance  of  (5.18)  requires  numerical  evaluation  of  the  Mittag-Leffler  function. 
In  general,  a  calculation  based  on  the  series  expansion  (3.19)  is  not  recommended. 
Instead,  the  Laplace  transform  of  the  Mittag-Leffler  law  can  be  inverted  by  contour 
integration  along  a  suitably  chosen  path  in  the  complex  plane;  details  are  deferred 
to  Section  5.5.  We  finally  arrive  at  the  alternative  expression 


K(-x)  = 


sin(i/7r) 


p-it/2 

Jo  1 


,-{x  tan  0)1/" 


vi r  J o  1  +  sin(20)  cos(^7r) 

which  is  evaluated  by  numerical  integration. 


d6, 


x  >  0,  0<v<l,  (5.19) 
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Figure  5.3:  Pareto  a  =  1.5  session  duration. 
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Figure  5.4:  Pareto  a  =  1.7  session  duration. 


115 


We  have  tested  approximation  (5.18)  for  system  utilizations  p  =  0.2,  0.5  and 
0.8.  Under  long-range  dependence  simulation  estimates  converge  very  slowly; 
moreover  confidence  intervals  based  on  the  regenerative  method  cannot  be  con¬ 
structed,  because  the  underlying  period  has  infinite  variance.  In  the  results  shown 
the  runs  were  109  time  slots  long,  and  by  that  time  the  estimates  had  stabilized. 
The  log-log  scale  plots  in  Figures  5.3  and  5.4  correspond  to  two  Pareto  distribu¬ 
tions  with  parameters  a  =  1.5  and  a  =  1.7.  Observe  that  the  heavier  a  =  1.5 
Pareto  tail  induces  larger  tail  probabilities  than  a  =  1.7,  at  the  same  system  uti¬ 
lizations.  In  both  Figures  5.3  and  5.4  we  see  that  simulated  and  approximate  values 
are  very  close,  suggesting  that  expression  (5.18)  provides  a  satisfactory  approxi¬ 
mation.  Note  also  the  almost  linear  shape  of  the  curves  in  log-log  scale,  reflecting 
the  power  law  asymptotics  of  the  queue  size  distribution  announced  in  [29,  39,  46]. 

Truncated  Pareto  When  a  follows  a  truncated  Pareto  distribution,  the  result¬ 
ing  M\G\oo  process  is  short-range  dependent.  Yet,  over  a  finite  range  of  time 
scales,  it  can  display  dependencies  similar  to  those  of  a  long-range  dependent 
M\G\oo  process.  Specifically,  for  1  <  a  <  2,  pick  some  N  =  2,  3, . . .  and  consider 
the  truncated  Pareto  distribution  on  {1,  2, . . . ,  N}  given  by 

(n~“-  N~°),  n  —  1,2, ....  JV.  (5.20) 

The  distribution  (5.20)  has  finite  support  and  clearly  satisfies  Assumption  (C). 
The  integer  parameter  N  provides  the  desired  flexibility  in  controlling  the  tail 
behavior  of  <7,  hence  the  dependencies  in  the  M\G\oo  arrivals.  When  N  —  2 
the  corresponding  rv  a  is  deterministic,  <7  =  2  a.s.  As  N  increases  the  second 
moment  of  the  session  duration  distribution  also  increases,  thus  leading  to  stronger 
dependencies  in  the  M\G\oo  process.  In  the  limit  as  N  goes  to  infinity  the  rv  a 
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Figure  5.5:  Truncated  Pareto  a  =  1.7,  N  =  50. 
converges  weakly  to  the  standard  Pareto  rv 

P  [a  >  n]  =  n~a,  n  =  1,2, , 

so  that  we  traverse  the  boundary  from  a  rv  a  with  finite  exponential  moment  to 
one  that  has  infinite  variance,  hence  to  a  corresponding  M\G\oo  process  that  is 
long-range  dependent. 

To  study  the  effect  of  the  truncation  level  N  on  the  queue  size  distribution  we 
carry  out  simulation  experiments  for  a  =  1.7  and  two  truncation  levels  N  =  50 
and  N  =  1000.  Results  are  collected  for  system  utilizations  p  =  0.2,  0.5,  0.8.  As 
approximation  (5.14)  does  not  assume  any  simple  closed  form  expression  in  this 
case,  we  calculate  the  various  quantities  entering  (5.14)  numerically. 

In  Figure  5.5  we  compare  simulation  results  for  truncation  level  N  =  50  with 
the  approximate  values  obtained  from  (5.14).  Confidence  intervals  widths  are 
not  shown,  since,  with  the  exception  of  the  three  points  at  the  bottom  of  the 
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plot,  they  were  well  within  10%  of  the  mean.  The  pairs  of  curves,  corresponding 
to  p  —  0.2,  0.5  and  0.8,  show  that  the  approximation  tracks  the  true  queue  size 
probabilities  satisfactorily,  especially  for  small  buffer  sizes.  However,  it  is  clear  that 
as  the  target  probabilities  become  smaller  and  the  buffer  size  of  interest  larger,  the 
quality  of  the  approximation  degrades. 

Further  increase  in  N  is  expected  to  be  even  more  revealing  of  the  limitations 
of  approximation  (5.14).  Note  that,  maintaining  a  =  1.7  and  increasing  N  from 
50  to  1000  results  in  a  small  increase  in  the  expectation,  from  E  [<r]  =  2.9  to 
E  [a]  =  3.035,  and  a  large  increase  in  the  second  moment,  from  E  [a2]  =  15.863  to 
E  [cr2]  =  42.49.  Figure  5.6  depicts  the  queue  size  probabilities  for  truncated  Pareto 
session  duration  with  a  =  1.7  and  N  =  1000,  at  system  utilizations  p  =  0.2,  0.5 
and  0.8.  Along  with  simulation  estimates  we  also  plot  the  values  obtained  from 
approximation  (5.14)  (labeled  SRD)  and  those  from  expression  (5.18)  (labeled 
LRD).  The  latter  is  appropriate  for  a  long-range  dependent  M\G\oo  process,  so 
it  does  not  strictly  apply  to  the  truncated  Pareto  setup.  However  (5.18)  becomes 
applicable  in  the  limit  as  the  truncation  level  N  goes  to  infinity.  Thus,  it  is 
representative  of  the  shape  to  which  the  simulation  curve  tends  as  N  grows  larger. 
From  Figure  5.6  it  becomes  clear  that  while  the  estimates  from  the  second  order 
interpolation  (5.14)  are  adequate  for  small  buffer  sizes,  they  fail  to  track  the  true 
probabilities  as  the  buffer  size  increases.  In  fact,  for  N  =  1000  and  within  the 
buffer  range  shown  in  the  figures,  it  is  the  curve  from  the  LRD  approximation 
(5.18)  that  is  closer  to  the  simulated  values. 

Obviously,  the  larger  the  variance  of  the  truncated  Pareto  rv,  the  closer  the 
simulation  curve  will  be  to  that  of  the  LRD  approximation  (5.18).  Moreover,  we 
see  that  for  a  fixed  variance  of  a  the  match  between  the  simulation  curve  and  that 
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Figure  5.6:  SRD  vs  LRD  approximation  for  truncated  Pareto  session  duration. 
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Figure  5.7:  Truncated  Pareto  a  =  1.7,  N  =  1000. 

of  the  LRD  approximation  is  better  at  smaller  buffer  sizes.  On  the  other  hand,  as 
the  buffer  size  grows  to  infinity,  the  LRD  approximation  (5.18)  will  also  eventually 
fail:  It  exhibits  a  hyperbolic  decay,  while,  as  follows  from  the  developments  in  [47], 
the  queue  size  distribution  induced  by  truncated  Pareto  session  durations  decays 
exponentially  fast.  This  can  also  be  verified  visually,  from  the  linear  shape  of 
the  rightmost  part  of  the  simulation  curves  in  Figure  5.6.  This  contrast  between 
the  truncated  Pareto  and  standard  Pareto  session  durations  is  more  evident  when 
comparing  the  log-log  plots  of  simulated  values  in  Figures  5.4  and  5.7.  Although 
the  shape  of  the  left  part  of  the  curves  in  Figure  5.7  is  almost  linear,  suggestive  of 
a  hyperbolic  initial  segment,  the  righmost  part  decreases  rapidly  and  the  curves 
become  concave,  in  accordance  with  the  anticipated  exponential  decay. 
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5.5  On  the  Mittag— Leffler  distribution 


In  Chapter  3  we  saw  that  the  heavy  traffic  buffer  content  distributions  for  long- 
range  dependent  M\G\oo  arrivals  are  given  through  Mittag-Leffler  special  func¬ 
tions.  Here,  we  leave  the  series  expansion  (3.19)  for  the  Mittag-Leffler  function 
aside,  and  discuss  instead  an  alternative  representation.  This  offers  an  interest¬ 
ing  interpretation  for  this  class  of  distributions  and  is  often  more  amenable  to 
numerical  calculations,  such  as  the  ones  performed  in  Section  5.4. 

It  is  known  [17,  p.  207]  that  for  0  <  v  <  1,  the  Mittag-Leffler  function  Eu(—x) 
is  completely  monotone  on  [0,  oo),  i.e.,  for  all  n  =  0, 1, . . . 

(— l)n  ~—Eu(—x)  >0,  x  >  0,  0  <  v  <  1. 

The  class  of  completely  monotone  functions  is  characterized  by  the  following  the¬ 
orem  [20,  p.  439],  due  to  Bernstein: 

Theorem  5.5.1  The  function  (j>  :  [0,  oo)  — >  R  is  completely  monotone  if  and  only 
if  it  is  of  the  form 

P+  oo 

(J)( A)  =  /  e~Xx  F{dx)  (5.21) 

Jo 

where  F  is  a  measure,  not  necessarily  finite,  on  [0,  oo). 


It  is  then  clear  that  for  0  <  v  <  1,  the  Mittag-Leffler  function  x  — >  Eu(—x), 
x  >  0,  admits  the  representation  (5.21)  and  hence  can  be  viewed  as  a  mixture  of 
exponential  distributions.  Its  corresponding  measure  F  in  (5.21)  is  determined  as 
follows: 

First,  we  recall  the  established  Laplace  transform  relations.  In  particular,  (3.20) 
and  (3.23)  provide  the  Laplace  transform  pair 

e~sxEJ-xu)  dx  =  -  - - - ,  s  >  0,  0  <  v  <  1;  (5.22) 

v  ’  s  IT  s~v  ’  v  ’ 
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this  formula  can  also  be  derived  from  [17,  Eq.  (18)  p.  209].  We  obtain  a  rep¬ 
resentation  of  the  form  (5.21)  for  Eu(—xu)  by  using  the  definition  of  the  Laplace 
inversion  integral 


E1.(—xv)  =  lim 


1 


vd+jY 


esx  - 


1 


_  _  ds,  x  >  0,  0  <  v  <  1,  (5.23) 

y->+oo2tt/  Jd_jY  s  1  +  s-"  -  V  7 

where  we  pick  d  to  be  any  strictly  positive  abscissa  on  the  real  axis.  The  denom- 


•  2m+l 


inator  1  +  s  v  has  roots  of  the  form  sm  =  e  3  v  n,  m  =  0,  ±1, . . .,  and  since 

1 2m  +  1| 


0  <  v  <  1  we  have 


>  1  for  every  m  =  0,  ±1, . . .,  so  that  there  is  no  root 


with  argument  in  [ — 7r,  7 r].  Thus,  the  integral  (5.23)  can  be  evaluated  along  the 
clockwise  path  from  H  to  A  shown  in  Figure  5.8,  namely 


Eu(—xu)  =  lim 


1 


1 


1 


R— >+oo  27 r  7 
r—>0 


ds, 


x  >  0,  0  <  v  <  1. 


(5,24) 


I H  5  1  +  5”  ' 

In  the  limit  as  the  radius  R  goes  to  infinity  the  contribution  of  the  arcs  CB  and 

GF  to  the  integral  above  is  zero,  by  Jordan’s 
lemma.  The  limiting  contribution  of  the  arcs 
HG  and  BA  is  also  zero,  because  their  length 
is  bounded  and  esx  is  also  bounded  along 
these  arcs.  Finally,  on  the  circular  arc  ED 
we  set  s  =  r  e 0  in  [ — 7r,  tt]  ,  and  see  that 
the  resulting  integrand  vanishes  as  r  — »  0. 
Thus,  in  (5.24)  only  the  integrals  along  the 
segments  FE  and  DC  remain  in  the  limit. 
Setting  s  =  ye~3n  and  s  =  yejn,  y  >  0,  for 
FE  and  DC,  respectively,  and  manipulating 
we  collect 


Figure  5.8:  Integration  path 


"+QO 


EJ-xA  = 


e  xy  fu(y)  dy,  x>0,  0  <  v  <  1, 


(5.25) 
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with 


fM  ■  = 


sin(z/7r) 

Tryl~u  (y2u  +  2 yv  cos(z/7r)  +  1)  ’ 


y  >  0,  0  <  v  <  1. 


(5.26) 


Expression  (5.19)  is  now  obtained  by  a  change  of  variable  tan#  =  yv  in  (5.25). 
Relations  (5.25)  and  (5.26)  provide  the  desired  interpretation  of  the  Mi t tag  -Leffler 
distribution  in  terms  of  an  infinite  mixture  of  exponentials.  The  specific  form 
of  (5.26)  shows  that  when  0  <  v  <  1,  the  density  function  fu(y)  weighting  an 
exponential  distribution  with  parameter  y  increases  to  infinity  as  y  goes  to  zero,  so 
that  the  Mi  flag  -Leffler  distribution  contains  no  single  dominant  exponential  with 
strictly  positive  parameter,  as  expected.  This  infinite  mixture  of  exponentials  is  to 
be  contrasted  with  the  classical  heavy  traffic  queue  size  distribution  under  short- 
range  dependence,  given  by  a  single  exponential;  the  one  to  which  expression  (5.25) 
collapses  as  v  — *  1. 

Representation  (5.25)  and  convergence  to  the  exponential  distribution  as  v  — »  1 
are  illustrated  in  Figure  5.9.  The  bottom  plots  show  the  density  function  fv(x) 
for  v  =  0.6,  0.8  and  0.99,  and  the  top  log-log  scale  plots  show  the  corresponding 
Mittag-Leffler  distributions  Eu(—xu).  The  dash-dotted  line  depicts  the  negative 
exponential  e~x .  It  is  clear  that  as  v  — »  1  the  density  fv{x)  tends  to  place  all  of  its 
mass  at  one,  i.e.,  at  a  single  exponential,  and  this  also  becomes  apparent  from  the 
top  plot  for  v  =  0.99.  In  addition,  we  observe  that,  even  for  v  =  0.99,  although 
e~x  and  Ely(—x1')  are  very  close  for  small  values  of  the  argument  x,  they  remain 
strikingly  different  as  x  grows  to  infinity. 
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v  =  0.6  v  =  0.8  v  =  0.99 


Figure  5.9:  Mittag-Leffler  distributions  Ev{—xl/)  and  spectral  densities  /„(x). 
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Chapter  6 


Conclusions 


Recent  measurements  studies  have  demonstrated  that  communication  networks 
carry  traffic  much  burstier  than  expected  (self-similar,  long-range,  or  subexpo- 
nentially  dependent).  These  findings  have  generated  interest  in  the  potential  im¬ 
plications  of  high  variability  and  dependence  on  network  performance. 

In  this  dissertation,  we  sought  to  understand  the  impact  of  (strong)  correla¬ 
tions  in  the  input  packet  stream  on  the  performance  of  a  single  network  multi¬ 
plexer.  This  was  modeled  as  a  discrete-time  queue  driven  by  a  family  of  M\G\oo 
correlated  arrival  processes.  Given  that  obtaining  exact  solutions  is,  in  general, 
extremely  difficult,  we  instead  focused  on  the  analysis  of  the  system  behavior  in 
two  asymptotic  regimes,  namely  light  and  heavy  traffic. 

In  heavy  traffic,  we  distinguished  between  M\G\oo  arrival  processes  with  short- 
and  long-range  dependence,  identifying  for  each  case  the  appropriate  heavy  traffic 
scaling  that  results  in  non-degenerate  limits.  The  resulting  limits  for  short-range 
dependent  inputs  involve  the  standard  Brownian  motion.  Of  particular  interest  are 
the  conclusions  for  the  long-range  dependent  case:  The  normalized  queue  length 
can  be  expressed  as  a  function  not  of  a  fractional  Brownian  motion,  but  of  an  ca¬ 
stable,  1/a  self-similar  independent  increments  Levy  process.  The  buffer  content 
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distribution  in  heavy  traffic  is  expressed  through  a  Mittag-Leffler  special  function 
and  displays  a  hyperbolic  decay,  of  power  1  —  a. 

Investigation  of  the  system  behavior  in  fight  traffic  reveals  the  effect  of  two 
aspects  of  the  M\G\oo  arrivals,  i.e.,  the  session  duration  distribution  G  and  the 
gradual  nature  of  the  inputs,  as  opposed  to  the  instantaneous  inputs  of  a  standard 
GI\GI\1  queue.  However,  the  arising  limits  cannot  be  fully  classified  in  terms  of 
short  vs  long-range  dependence  property  of  the  M\G\oo  process,  hence  demon¬ 
strating  that  the  latter  is  not  the  only  factor  that  impacts  performance. 

Exploiting  the  results  above,  we  proposed  a  family  of  heuristic  approximations 
for  a  multiplexer  with  M\G\oo  inputs.  These  interpolation  approximations  were 
developed  by  combining  the  asymptotic  characterizations  of  the  buffer  content 
distribution  under  heavy  and  fight  traffic  conditions  and  are  applicable  to  all  traffic 
intensities.  For  several  common  pmfs  G  the  approximants  assume  a  simple  final 
form,  and  are  capable  of  providing  quick  and  reliable  estimates  of  the  buffer  content 
distribution,  especially  for  small  buffer  sizes. 

In  closing,  we  mention  two  open  questions.  First,  Corollary  4.3.2  lead  us  to 
conjecture  that,  in  light  traffic,  if  c  =  1  and  the  session  duration  distribution 
has  a  regularly  varying  tail  of  order  —a,  with  1  <  a  <  2,  then  the  buffer  con¬ 
tent  distribution  also  exhibits  a  power  law  behavior  of  the  form  A“.  Second,  the 
stochastic  comparisons  discussed  here  could  be  complemented  by  stochastic  com¬ 
parisons  within  the  family  of  M\G\oo  inputs.  For  two  discrete-time  queues,  with 
c  =  1,  each  driven  by  instantaneous  inputs  (2.16),  representation  (2.38)  shows 
that  a  convex  comparison  between  the  respective  distributions  of  a  translates  into 
a  strong  stochastic  comparison  between  the  corresponding  queue  lengths.  This 
can  be  viewed  as  a  discrete-time  analog  of  the  well  known  result  for  the  wait- 
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ing  time  in  a  standard  M\G\1  queue  and  is  another  instance  of  a  folk  theorem 
of  queueing  “variability  increases  delays”.  We  conjecture  that  the  folk  theorem 
holds  for  the  gradual  M\G\oo  inputs  as  well  (when  c  =  1),  and  simulation  results 
indeed  suggest  that  a  convex  stochastic  ordering  between  session  durations  leads 
to  a  strong  stochastic  comparison  between  the  queue  lengths.  However,  as  the  ana¬ 
log  to  (2.38)  is  not  available,  the  conjecture  remains  unsettled,  and  may  require 
multidimensional  stochastic  comparisons  between  M\G\oo  arrival  vectors. 
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Appendix  A 


Asymptotic  invertibility  of  regularly  varying 
functions 


We  establish  the  asymptoptic  invertibility  property  given  in  Proposition  3.3.1  by 
relying  on  the  following  facts: 


Proposition  A.l  Consider  a  Borel  measurable  function  [3  :  R+  — »  R  such  that 
lim  f3{x)  =  0.  With  A  >  0  and  y  >  0,  the  equation 

x—>oo 

y  =  exp  (^JX  1  d?j  ,  x>A  (A.l) 

has  a  unique  solution  x  :=  x(y)  for  all  y  large  enough.  Moreover, 


lim 

y— >-oo 


y) 

x(y) 


=  7,  7  >  0, 


(A.2) 


or  equivalently,  the  mapping  y  — >  x(y)  is  regularly  varying  of  order  1,  i.e.,  x(y)  ~ 
yv(y)  (y  — >  oo)  for  some  slowly  varying  function  v  :  R+  — >  R+. 


Proof.  Set 

B{x)  :=  1  - 13 ^  dt ,  x  >  A 

J  A  t 
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and  pick  £  in  (0, 1).  Since  lim (3(x)  =  0,  there  exists  A *  =  A*(e)  >  A  such  that 

#—>•00 


1  <  1  -  P{t) 

t  t 


< 


1  +  £ 


t  >  A*. 


(A.3) 


It  is  straightforward  to  see  that  lim  B(x)  —  oo  and  B(A)  =  0,  and  by  continuity, 

#—>•00 

the  range  of  x  — >  exp (B(x))  contains  the  semi-  infinite  interval  [l,oo).  We  also 
conclude  from  (A.3)  that  x  — »  B(x)  is  strictly  monotone  increasing  on  the  interval 
[A*,oo),  and  the  existence  and  uniqueness  of  a  solution  to  (A.l)  follows  when¬ 
ever  y  >  y*  with  y*  :=  exp(B(A*)).  The  solution  mapping  y  — >  x(y)  is  strictly 
increasing  on  [y*,  oo). 

We  now  turn  to  proving  (A. 2).  There  is  nothing  to  prove  when  7  =  1.  With 
7  >  1,  whenever  y  >  y*,  we  get 

7  =  —  =  exp(£(x(qy))  -  B(x(y )))  =  exp  (  [  - — dt 

V  \Jx(y)  1 

and  the  use  of  the  inequalities  (A.3)  yields 


^  ■ 

1 - 

V 
C 

V 

U) 

1 

'x(nr  y)~ 

Jcy 

l+£ 


or  equivalently, 


i  x('yy)  i 

ql+J  <  <  -y  — 

x{y ) 


Letting  y  go  to  infinity  in  (A. 4),  we  conclude 


(A.4) 


A-  ^  r  ■  f  x(w)  ^  r  ®(7 V)  ^  W- 

qi+s  <  innmt  — — —  <  hmsup  — — —  <  q1^, 

v+ °°  x{y)  y->  oo  x(f/) 

and  (A. 2)  is  obtained  as  we  note  that  £  is  arbitrary  in  (0, 1).  The  case  q  <  1  is 
handled  in  a  similar  way;  details  are  omitted  in  the  interest  of  brevity.  ■ 
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Lemma  A. 2  Consider  slowly  varying  functions  u,  w  :  R+  — »  R+;  such  that 
u(x)  ~  w(x)  (x  — >  oo),  and  let  a  >  1.  For  any  sequences  {£r,  r  =  1,2,...} 
and  {rjr,  r  =  1,2,...}  with  lim  ( r  =  lim  r/r  =  oo  such  that 

r— >oo  r—>  oo 

lim  r(fau((r)  =  lim  rrjfaw(r)r )  =  K  (A. 5) 

r— too  r— too 

for  some  finite  constant  K  >  0,  it  holds  that  ( r  ~  r]r  (r  — >  oo). 


Proof.  We  first  look  at  the  special  case  when  u 
(A. 5)  implies 


CMC- 


lim 

r—yoo  r]r  aU[rjr 


=  1. 


w,  in  which  case  condition 


(A.6) 


We  refer  to  the  proof  of  Lemma  3.8.1,  where  we  introduced  the  asymptotically 
equivalent  representation  (3.66)  of  the  slowly  varying  function  u.  Substituting 
(3.66)  in  (A.6),  we  see  that 


lim  exp 

r—¥  oo 


a  —  eft) 
t 


dt  -|~ 


=  1, 


or,  equivalently, 


lim 

r—¥  oo 


rVr  a  —  eft) 


dt 


=  0. 


(A.7) 


Pick  S  in  (0,o;).  Because  lim  fr  =  lim  r/r  =  oo,  there  exists  r$  such  that  for  r  >  r$ 

r— too  r—t’oo 

we  have  |e(t)|  <  5  whenever  t  >  min((‘r,  rjr).  Thus, 


f a  —  5) 


In^ 

< 

r a  - e(t)  m 

hr  t 

r  >  r5, 


and  combining  this  last  inequality  with  (A.7)  we  obtain  the  desired  conclusion 
lim  r/r  /  (r  =  1. 

r—¥  oo 

In  general,  when  u  and  w  are  not  necessarily  equal,  we  note  the  easy  relation 

Cau((r)  =  r(fau((r)  _  W(r]r)  r  =  12 

rj~au(r]r)  rq~aw(r]r)  ufqr)  ’ 
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Condition  (A.  5)  and  the  asymptotic  equivalence  of  u  and  w  together  imply  that 
the  relation  (A. 6)  still  holds,  and  the  conclusion  (r  ~  yr  (r  — >  oo)  follows  from  the 
hrst  part  of  the  proof.  ■ 


Proposition  A. 3  Consider  a  slowly  varying  function  u  :  R+  — »  R+,  and  let 
a  >  1.  For  any  sequence  {Cr,  r  =  1,  2, . . .}  with  lim  (r  =  oo  such  that  (A. 5)  holds, 

r—>oo 

we  have 


Cr 


raw{r ) 


(r  — )•  cx)) 


for  some  slowly  varying  function  w  :  R+  — »  R+. 


Proof.  We  go  back  to  the  proof  of  Lemma  3.8.1,  where  we  introduced  the 
asymptotically  equivalent  representation  (3.66)  of  the  slowly  varying  function  u.  In 
view  of  Lemma  A. 2,  it  suffices  to  consider  a  sequence  {Cr,  r  =  1,  2, . . .}  determined 
by  the  relations 

rC~“  •  cexp  (/  dt'j  =  K,  r  >  r*  (A. 8) 

for  some  r*  large  enough,  with  constants  A  >  0  and  c  >  0,  and  Borel  mapping 
£  :  R+  — >  R  such  that  lim  e{t)  =  0.  We  can  write  (A. 8)  in  the  equivalent  form 

t—¥  OO 


Br « 


exp 


{jCr  i  -  m 


r  >  r* 


with 

B  :=  (tm)  a  and  ^  :=  f  -  °' 

Hence,  by  Proposition  A.l,  for  large  enough  r  we  see  that  fr  is  the  unique  solution 
x(y )  of  the  equation  (A.l)  with  y  =  Br* .  By  the  second  part  of  Proposition  A.l, 


131 


we  have 


( r  =  x{Bra)  ~  Brav{Bra)  (r  — *  cx)) 

for  some  slowly  varying  function  v  :  1R+  — »  H+.  The  desired  conclusion  is  now 
immediate  once  we  note  that  the  mapping  w  :  x  — >  Bv(Bx^)  is  slowly  varying 
whenever  v  is.  ■ 
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Appendix  B 


Stochastic  orderings 


We  collect  here  some  definitions  and  properties  concerning  stochastic  orderings. 
The  material  is  drawn  mostly  from  [59],  where  additional  information  is  available. 

Throughout  let  X  and  Y  denote  two  IR, -valued  rvs,  and  let  V  denote  the  set 
of  all  probability  distribution  functions  of  IR-valued  rvs. 

Definition  B.l  Let  X  and  Y  have  distribution  functions  F  and  G,  respectively. 
We  say  that  X  is  stochastically  smaller  than  Y ,  and  write  X  <st  Y ,  or,  equiva¬ 
lently,  F  <st  G,  if 

F(x)  >  G(x),  x  G  R . 

Proposition  B.2  It  holds  that  X  <st  Y  if  and  only  if 

E[#X)]<E[#y)]  (B.l) 

for  all  non- decreasing  functions  <f  :  IR  — *  R  for  which  the  expectations  in  (B.l) 
exist. 

Proposition  B.3  Let  {Fn,  n  =  1,2, . . .}  and  { Gn ,  77.  =  1,2,...}  be  two  subsets 
ofV  such  that  Fn  F  and  Gn  =^„,  G  for  limits  F  and  G  in  V,  respectively. 
If  Fn  <st  Gn  for  all  n  =  1,2,...  then  F  <st  G. 
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Definition  B.4  We  say  that  X  is  smaller  in  the  convex  stochastic  ordering  than 
Y ,  and  write  X  <cx  Y ,  if 

E[#X)]<E[#y)]  (B.2) 

for  all  convex  functions  <f  :  R  — *  1R  whenever  the  expectations  exist  in  (B.2). 

We  similarly  define  the  increasing  convex  stochastic  ordering. 

Definition  B.5  We  say  that  X  is  smaller  in  the  increasing  convex  stochastic  or¬ 
dering  than  Y,  and  write  X  <icx  Y ,  if 

EWX)]<EWF)]  (B.3) 

for  all  increasing  convex  functions  4>  :  R  — »  H  whenever  the  expectations  exist  in 
(B.3). 

The  <icx  ordering  admits  the  following  characterization  [59,  p.  8]. 

Proposition  B.6  It  holds  that  X  <iCX  Y  if  and  only  if  X  is  smaller  in  mean 
residual  life  than  Y ,  i.e., 

E  [( X  -  x)+]  <  E  [(y  —  x)+]  ,  X  e  R 

provided  the  expectations  above  are  finite. 

Consider  now  the  Lindley  recursions 

wi+i  =  iw{nX)  +  Xn+1}+  and  =  [n;(y)  +  Yn+1}+  n  =  0, 1, ,  (B.4) 

with  initial  conditions  Wq  ;  and  w(  ,  these  are  independent  of  the  driving  se¬ 
quences  of  i.i.d.  rvs  {Xn,  n  =  1,2,...}  and  { Yn ,  n  =  1,2,.. .},  with  generic  rvs  X 
and  Y,  respectively.  Let  -<  denote  either  <st  or  <iCX,  and  let  {wn,  n  =  0, 1, . . .} 
be  either  {w^\  n  —  0, 1, . . .}  or  \  n  =  0, 1, . . .}. 
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Proposition  B.7  [59,  p.  79]  If  wq  -<  w\  in  (B.f)  then  for  all  n  =  0,1,...,  we 
have 

wn  "" '  n  —  1  ■ 


Moreover,  if  the  stationary  rv  exists,  then 

wn -<  w^,  n  —  0,1, , 

where,  if  -<  denotes  <iCX,  it  is  further  assumed  that  Woo  has  finite  expectation. 
Proposition  B.8  [59,  p.  80]  If  X  -<Y  and  WqX>  -<  in  (B.f)  then 

(X)  y  (Y)  n  1 

Wn+l^Wn+V  n  =  0,1,.... 

Moreover,  if  the  corresponding  stationary  versions  w and  w exist,  then  X  -< 
Y  is  sufficient  to  ensure  that 


where,  if  -<  denotes  <iCX, 


W{X)  W(Y) 

^OO  ^  ^OO 


it  is  further  assumed  that  w has  finite  expectation. 
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